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This article is devoted to the study of the critical properties of classical XY and Heisenberg frus- 
trated magnets in three dimensions. We first analyze the experimental and numerical situations. 
We show that the unusual behaviors encountered in these systems, typically nonuniversal scaling, 
are hardly compatible with the hypothesis of a second order phase transition. Moreover, the fact 
that the scaling laws are significantly violated and that the anomalous dimension is negative in 
many cases provides strong indications that the transitions in frustrated magnets are most prob- 
ably of very weak first order. We then review the various perturbative and early nonperturbative 
approaches used to investigate these systems. We argue that none of them provides a completely 
satisfactory description of the three-dimensional critical behavior. We then recall the principles 
of the nonperturbative approach — the effective average action method — that we have used to 
investigate the physics of frustrated magnets. First, we recall the treatment of the unfrustrated 
— O(N) — case with this method. This allows to introduce its technical aspects. Then, we show 
how this method unables to clarify most of the problems encountered in the previous theoretical 
descriptions of frustrated magnets. Firstly, we get an explanation of the long-standing mismatch 
between different perturbative approaches which consists in a nonperturbative mechanism of anni- 
hilation of fixed points between two and three dimensions. Secondly, we get a coherent picture of 
the physics of frustrated magnets in qualitative and (semi-) quantitative agreement with the numer- 
ical and experimental results. The central feature that emerges from our approach is the existence 
of scaling behaviors without fixed or pseudo-fixed point and that relies on a slowing-down of the 
renormalization group flow in a whole region in the coupling constants space. This phenomenon 
allows to explain the occurence of generic weak first order behaviors and to understand the absence 
of universality in the critical behavior of frustrated magnets. 
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I. INTRODUCTION 



Understanding phase transitions and, specifically, crit- 
ical phenomena has been one of the central issues of sta- 
tistical mechanics during these last decades Q and the 
field theoretical renormalization group (RG) approach to 
these phenomena has been one of the great successes of 
theoretical physics. This is so true that it is generally 
believed that, apart from specific problems — disordered 
and glassy systems for instance — , an almost complete 
understanding of the physics occurring at a phase tran- 
sition has been reached. This is certainly due to the fact 
that it is indeed the case for all the systems belonging 
to the so-called Wilson-Fisher universality classes of d- 
dimensional systems whose symmetry breaking scheme 
is given by O(N) — > 0(N — 1). In fact, although they 
have become the archetype of systems displaying critical 
phenomena well described by perturbative field theoreti- 
cal approaches, these O(N) symmetric systems turn out 
to be exceptions rather than the rule. For most systems a 
quantitative and, even sometimes, a qualitative descrip- 
tion of the critical physics is either still lacking or very 
difficult to obtain by perturbative RG methods. This is 



the case, for instance, in the Potts model m mag- 

netic systems with disorder Qjin superconductors 
in Josephson junction arrays 0, in He3 H,@, in smec- 
tic liquid crystals jlCj . in electroweak phase transitions 
[inil2| and in frustrated magnets like helimagnets or ge- 
ometrically frustrated magnets (triangular for instance) 
which are our main purpose in this article (see jl3j | for a 
review) . 

Actually, it should not be surprising that a qualitative 
difference exists between the critical behaviors of systems 
belonging to the O(N) universality class and the oth- 
ers: among the systems where the order parameter has 
N real components, O(N) corresponds to the maximal 
symmetry and, thus, to the simplest situation. Think, 
for instance, at a unit-norm constraint imposed on the 
microscopic degrees of freedom (S 2 = 1): the maximal 
symmetry compatible with it is indeed O(N). From a 
perturbative point of view, this means that the Ginzburg- 
Landau- Wilson (GLW) Hamiltonian of an O(N) sym- 
metric model involves only one (marginal) interaction 
term — {4> 2 ) 2 — and, thus, only one coupling constant. 
As a consequence, the perturbative RG flow of the criti- 
cal theory takes place in a one-dimensional space of cou- 
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pling constant and is thus simple. In particular, only 
one nontrivial perturbative fixed point can exist [14} . On 
the contrary, the Hamiltonian of systems having a TV- 
dimensional order parameter and displaying a symmetry 
group G smaller than O(N) involves also terms that ex- 
plicitly break O(N). It thus contains several interaction 
terms and, therefore, several coupling constants. The 
RG flow then takes place in a multi-dimensional space 
and is thus far less simple: it can, in particular, involve 
fluctuation-induced first order transitions — runaway in 
a region of instability — and several fixed points with 
different symmetries. Universality itself is not guaran- 
teed in these systems since the basins of attraction of the 
fixed points can be highly nontrivial. 

Of particular interest for us, it is generically observed 
in these systems that, by varying N and/or d, the critical 
physics changes qualitatively: low dimensions (d — > 2) 
and large number of spin components (N — > oo) favor 
smooth fluctuations and, thus, second order phase tran- 
sitions, while larger dimensions (d — ■> 4) and small N 
(N ~ 1) f avor larger fluctuations and thus first order 
transitions |206| |. Therefore, in many systems — and 
notably in frustrated systems — the critical behavior 
changes qualitatively i) for the physically interesting val- 
ues of N — N ~ 1 — when the dimension varies between 
d = 2 and d = 4 ii) at fixed dimension when the number 
N of components varies between N — oo and N ~ 1. 
Thus, the different perturbative approaches are in the 
worst possible position: it is quite difficult to obtain def- 
inite conclusions in d = 3 and for N ~ 1 from extrap- 
olations of perturbative results even if they are valid in 
the domains where they have been established: d = 2 for 
the Nonlinear Sigma (NLtr) model, d — 4 for the GLW 
model and for large N in a 1/N expansion. This is one 
of the reasons why, after more than twenty five years 
of considerable efforts, the situation is still not clear for 
most three-dimensional systems that do not belong to 
the 0(N)/0(N — 1) universality classes. 

Let us now discuss two concrete problems encountered 
in the perturbative RG studies performed on the three- 
dimensional systems we are interested in. First, the com- 
putational difficulties encountered in perturbation theory 
are non negligible. Within the NLer model approach, 
the series are generally considered as useless due to the 
lack of Borel-summability (see however |Tq). Within 
the GLW approach, the perturbative computations al- 
most always call for resummation procedures. In gen- 
eral, these procedures are not as easy as they are in the 
O(N) fO(N — 1) models. The series are either not proven 
to be Borel summable or are even suspected to be non 
Borel summable. This is the root of a lot of difficulties en- 
countered in this approach (see [IB, H3, EH H3, ; 
for a review in the case of the diluted Ising model, see 
|23|. for the presence of non analyticities in perturbative 
series, see j24[ H^, |2(| |22) and, for a general review, see 
pfij). The second point is more conceptual: although it 
is generically possible to perform a perturbative expan- 
sion of the critical theory around d = 2 — within the 



NLcr model approach — and around d = 4 — within 
the GLW approach — it has not been possible to relate 
these two expansions within the usual field theoretical ap- 
proach (except for large enough N where the 1/N expan- 
sion allows to recover, at leading orders, the perturbative 
results obtained in the NLcr and GLW approaches). From 
this point of view, even high-order perturbative calcula- 
tions performed in the GLW model do not help since the 
perturbative expansion cannot be extrapolated down to 
d = 2 for N > 2. For instance, the critical exponent v 
diverges in d — 2 when it is calculated as a power series 
in e = 4 — d. This fact, which is not crucial for systems 
whose critical behavior does not change qualitatively be- 
tween d = 2 and d = 4 (e.g. the 0(N) fO(N - 1) models) 
forbids for the others to obtain a completely coherent 
picture of the physics between d = 2 and d = 4. Most 
of the time one of the perturbative approach — usually 
the NLcr one — is dismissed without real justifications 
and the other is blindly trusted. Since all the RG equa- 
tions are smooth in N and d, it is not clear if and also, 
why, this procedure is legitimate. It would, of course, be 
much more satisfactory to have a unified approach not 
linked to a particular value of d or TV and that allows to 
interpolate between both approaches. 

All these drawbacks of the usual perturbative RG 
methods call for a nonperturbative approach. Such an 
approach is, in fact, already known and its foundations go 
back to Kadanoff and Wilson with the idea of block spin 
and effective, scale-dependent theory 0, |2^| . It is some- 
times called the exact renormalization group method 
but we prefer to call it the nonperturbative renormal- 
ization group (NPRG) method (for contributions of dif- 
ferent authors to the early attempts to use NPRG, see 
pjol I31II32I l33ll . for an exhaustive bibliography of the sub- 
ject, see )• This idea has been turned into an efficient 
computational tool du ring th e last ten years, mainly by 
Ellwanger M JM pjjj Hall , Morris Hi EH and Wet- 
terich [I^, |43, |4J, [45j ■ It has allowed to determine the 
critical exponents of the O(N) models with high preci- 
sion without having recourse to resummation techniques 
EH El MJM \m EU. It has also allowed, for the 
first time |5lf~to relate, for any N, the results of the 
0(N)/0(N - 1) model obtained near d = 4 and d = 2, 
a fact of major importance for our purpose. Also impor- 
tant for the present purpose, it has allowed to tackle with 
genuinely nonperturbative situations. For instance, the 
Berezinskii-Kosterlitz-Thouless phase transition Ell 
has been recovered directly from a study of the GLW 
mo del, i.e. without introducing explicitly the vortices 
[54ll55| . To cite just a few other successes of this method, 
let us mention low-energy Quantum Chromodynamics 
[45) . the abelian Higgs model relevant for superconduc- 
tivity j5||, the study of the Gross-Neveu model in three 
dimensions [Ell E3j phase transitions in He3 the 
study of cubic anisotropy in all dimensions as well as the 
randomly diluted Ising model [60l |. the two-dimensional 
Ising multicritical points |6l|. etc. 

In this article, we study by means of NPRG methods 
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one of the most famous systems exhibiting the changes 
of critical behavior previously described: it is the system 
consisting of XY or Heisenberg spins on the triangular 
lattice (stacked triangular in d — 3) with antiferromag- 
netic nearest neighbor interaction (Sectional). This sys- 
tem is the archetype of frustrated spin systems and is 
supposed to be in the same universality class as another 
set of frustrated magnets: the helimagnets. Almost all 
these systems have been intensively studied both numeri- 
cally and experimentally these last twenty five years (see 
Section However, their behavior remains unclear 

and displays quite unconventional features. For instance, 
almost all experiments exhibit scaling laws around the 
transition temperature — which suggests a second order 
phase transition — but with critical exponents that de- 
pend on the particular material studied, on microscopic 
details, etc, which is incompatible with the standard phe- 
nomenology of a second order phase transition. In some 
experiments or numerical simulations, the scaling laws 
are sometimes significantly violated while the anomalous 
dimension r\ is found negative, a fact forbidden by first 
principles if the theory is <fi A GLW-like (see the follow- 
ing). The theoretical situation in these systems is also 
not clear from the perturbative point of view (Sections 
II VI and Iv )) : first, independently of the experimental con- 
text, the results obtained within the usual perturbative 
approaches — in dimensions d = 2 + e and d = 4 — e — 
conflict. Second, neither the low-temperature expansion 
around d = 2 nor high-order weak-coupling calculations 
performed around d = 4 or directly in d = 3 succeed in 
reproducing satisfactorily the phenomenology. We show, 
in this article, that the NPRG approach ( Section IVI|) to 
frustrated systems fSection IVTIf) clarifies almost entirely 
the situation. First, it allows to smoothly interpolate 
between d = 2 and d = 4 and to clarify the mismatch 
between these approaches. In particular, a mechanism of 
annihilation of fixed points, already identified for a long 
time around 4 — e dimensions for N ~ 21.8 is shown to 
operate around two dimensions for N ~ 3 nonpertuba- 
tively with respect to the low-temperature approach of 
the NLcr model jS^, H^- This explains the irrelevance 
in d = 3 of the 0(4) fixed point obtained within a low- 
temperature approach in d = 2 + e. Second, our approach 
provides a description of the physics in d = 3, in terms 
of weakly first order behaviors, compatible with the phe- 
nomenology fSections IVIIII and IIX|I . In this respect, an 
important feature of our work is that it explains the oc- 
curence of scalin g in frustrated magnets without fixed or 
pseudo- fixed |6J, |65, 66] point. This phenomenon relies 
on a slowing-down of the RG flow in a whole region in 
coupling constants space. This allows to explain one of 
the most puzzling aspect of the critical physics of these 
systems, i.e. the occurence of scaling without universal- 
ity. We discuss (Section[Xj possible experimental and nu- 
merical tests of our scenario. We then comment (Section 
IXI|) the consequences of our work for the perturbative ap- 
proaches that have been used to investigate the physics 
of frustrated magnets. Finally, we give our conclusions 



(Section EH- 

II. THE STA MODEL AND GENERALIZATION 

A. The lattice model, its continuum limit and 
symmetries 

We now describe the archetype of frustrated spin sys- 
tems, the Stacked Triangular Antiferromagnets (STA). 
This system is composed of two-dimensional triangular 
lattices which are piled-up in the third direction. At 
each lattice site, there is a magnetic ion whose spin is 
described by a classical vector. The interaction between 
the spins is given by the usual lattice Hamiltonian: 

// ^./,,.S,..S. (1) 

(H) 

where, depending on the anisotropics, the Si are two or 
three-component vectors and the sum runs on all pairs of 
nearest neighbor spins. The coupling constants equals 
J_l for a pair of sites inside a plane and Jy between planes. 

The interactions between nearest neighbor spins within 
a plane is antiferromagnetic, i.e. J± > 0. This in- 
duces frustration in the system and, in the ground state, 
gives rise to the famous 120° structure of the spins, see 
Fig. [Qi. This nontrivial magnetic structure is invariant 

• S V -% : ^ % 

S '% y % if- y 



>" ^ | ^ ^ 




b) 



FIG. 1: The ground state configurations a) of the spins on 
the triangular lattice and b) of the order parameter made of 
two orthonormal vectors. The plaquettes, which constitute 
the magnetic cell, are indexed by I and are shaded. 
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under translations of length y3 times the initial lattice 
spacing. The magnetic cell, indexed by /, which is repli- 
cated all over the system, is a plaquette of three spins 
S[, Si and Si, see Fig. UK- 

Note that the nearest-neighbor out-of-plane interaction 
Jy is, depending on the compounds, ferro- or antiferro- 
magnetic, but the two cases can be treated simultane- 
ously since no extra frustration appears through this 
interaction. Finally, interactions between more distant 
spins (next-to- nearest neighbors, etc) also exist but are 
neglected in the following since they are supposed to be 
irrelevant. 

There have been numerous derivations of the long dis- 
tance effective field theory supposed to describe the criti- 
cal physics of this system [67L I68L 1691 17Q| . We here sketch 
the derivation which is the most appropriate for our pur- 
pose. The Hamiltonian (|TJ has the usual rotational sym- 
metry acting on the spin components: 0(2) or 0(3) for 
XY or Heisenberg spins, respectively. To identify the 
order parameter, it is also necessary to consider the sym- 
metry of the magnetic cell. For the triangular lattice, 
this is the C?.v group that interchanges the spins inside a 
plaquette 207]. 

The identification of the order parameter is close in 
spirit to what is done in the nonfrustrated case, e.g. for 
the antiferromagnets on a square lattice. At zero temper- 
ature, the sum of the three spins for a given plaquette /: 

E 7 = S[ + Si + Si (2) 

is vanishing (Let us note that E 7 is analogous to the 
local magnetization of nonfrustrated antiferromagnets — 
E 7 = S[ + Si in this last case — that also vanishes in 
the ground state). In average, this is also the case at any 
finite temperature so that the thermal average of: 

E = ^E 7 , (3) 
i 

where the sum runs on all plaquettes, cannot be an order 
parameter: the associated modes are never critical. We 
therefore replace E 7 by its average value: 

E 7 -» (E 7 ) = (4) 

which is equivalent to freezing the fluctuations of the 
spins inside each plaquette. The constraint E 7 = is 
called the "local rigidity constraint" . Having eliminated 
E, we keep only two vectors per plaquette (0{, 4>i) which 
represent the local order parameter. For 02, we choose 
one of the spins of the plaquette, see Figs.QJi and b. For 
the other, 0{, we choose the linear combination of the 
spins which is orthogonal to <f>i and of unit norm, see 
Fig. The local order parameter thus obeys on each 
plaquette: 

l\4] = ki with ye{l,2}. (5) 



The dihedral (0 7 , 0/,) plays a role analogous to the stag- 
gered magnetization in the nonfrustrated case. 

As usual, once the model is reformulated in terms of 
its order parameter, the effective interaction — from 
plaquette to plaquette — becomes ferromagnetic, see 
Fig. By taking the dihedral (0 7 ,02) on the center 
of the plaquette /, we indeed find that it interacts fer- 
romagnetically with the dihedral (0f, 0J 7 ) defined on the 
center of the plaquette J — the plaquettes / and J being 
nearest neighbours — such that 0{ interacts only with 0f 
and 02 only with 0J 7 . A more detailed analysis shows that 
the two vectors 0{ and 02 play symmetric roles jg^- As 
a consequence, the effective Hamiltonian reads: 

H = -JJ2 ($[4( + $ I 24i) (6) 
(i, J) 

with the same coupling constant J > for the 0f's and 
for the 4>i 's. Moreover, since the anisotropies resulting 
from the stacked structure of the lattice are supposed 
to be irrelevant, we take the same coupling constant for 
the interactions inside a plane and between the planes. 
The continuum limit is now trivial and proceeds as in the 
usual ferromagnetic case. The effective Hamiltonian in 
the continuum thus writes, up to constants: 

H = -Jd d x ((50! (x)) 2 + (50 2 (x)) 2 ) (7) 

with the constraint that 0i and 02 are orthonormal. This 
model is called the Stiefel Vjv,2 model with N = 2 in the 
XY case and N — 3 in the Heisenberg case. In Vm,2 the 
index 2 means that we are considering two orthonormal 
vectors 0i and 02. 

It is convenient to gather the vectors 0i and 02 into a 
rectangular matrix: 

$ = (01,02) (8) 

and to rewrite H as: 

H = —J rf d xTr (9 *$(x).9<i>(x)) (9) 

where = 

In the following two sections, we consider successively 
the case of Heisenberg and XY spins. 

B. The Heisenberg case 

In this case, H is invariant under the usual left 0(3) 
rotation and inversion group acting on the spins: 

$' = R<f>, R £ 0(3) . (10) 

It is also invariant under a right 0(2): 

= $U, U e 0(2) . (II) 
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This last symmetry encodes the fact that <p\ and 4>2 play 
the same role which, itself, is reminiscent of the C% v 
symmetry of the triangular plaquette. The system is 
thus symmetric under G — 0(3) x 0(2). In the low- 
temperature phase, a typical ground state configuration 
is given by (see Fig. 

fl 0\ 

$o « 1 . (12) 

It is symmetric under the diagonal group — ^(2)di a g ~~ 
built from the right 0(2) and from a particular left 0(2) 
in 0(3): 

(e cos 9 — sin # \ / n ■ n \ 
esmO cos6> $ • a a ) ( 13 ) 
o oi/ \- smd cos0 J 

where e = ±1 encodes the Z2 part of 0(2)^^. Apart 
from the previous Z2 contained in the 0(2)^^, another 
Z2 is also left unbroken. It is the combination of a Z2 
included into the right 0(2) of G, Eq. fTT|) . and of a 
rotation of 7r around the x-axis contained in the rotation 
group 50(3) of G. Thus, G is spontaneously broken 
down to H = Z2 x 0(2)jj a g. As a consequence, the 
symmetry breaking scheme reads: 

G = 0(3) x 0(2) -» H = Z 2 x 0(2) diag (14) 

which is often referred to, once all the Z2 groups have 
been cancelled, as the 50(3) x 50(2)/50(2) model. 

Here appears the main feature of frustrated magnets: 
the 50(3) group is fully broken in the low-temperature 
phase whereas it is only broken down to 50(2) in non- 
frustrated magnets. This has two important conse- 
quences that are at the very origin of the nontrivial crit- 
ical behavior encountered in frustrated magnets. 

First, there are three Goldstone modes in the broken 
phase instead of two in the nonfrustrated case. This im- 
plies a physics of spin waves different from that of the 
0(3)/0(2) model. Second, the order parameter space 
50(3) having a nontrivial first homotopy group 

tt 1 (50(3)) = Z 2 (15) 

there exist stable nontrivial topological configurations 
called vortices. Because of the Z2 homotopy group, only 
one kind of vortex exists, contrarily to the well-known 
case of XY ferromagnets where there are infinitely many 
different kinds of vortices, each one being indexed by an 
integer, the winding number. 

It has been established firstly by Kawamura and 
Miyashita E2 that the existence of vortices is impor- 
tant at finite temperature in two dimensions. This has 
been largely confirmed by subsequent works studying 
the temperature dependence of thermodynamical quanti- 
ties such as the correlation length, the spin-stiffness, etc 
El El El El El El Actually, although this has not 



been directly established, they certainly also play an im- 
portant role for the critical physics of the STA in three 
dimensions. A simple argument allows to argue to that 
end: let us go back on the lattice and introduce, on each 
plaquette /, together with (f>[ and 2 , a third vector 0| 
defined by: 

$i = ${ A $i . (16) 
Let us then gather them into a 3x3 matrix: 

= • (17) 

Since (cf>[, <\>\ , (^3) are three orthonormal vectors, one has 
*$ / $ / = t and, therefore, $ J is a 50(3) matrix. This 
allows to rewrite the Hamihonian JfJ) on the lattice as: 

H = -^2Tr(V (18) 

(i, J) 

where V is a diagonal matrix of coupling constants that 
characterizes the interaction between the ^('s, between 
the <j>2 's and between the ^3's. One deduces from the mi- 
croscopic derivation that V = diag(J, J, 0), i.e. that the 
interaction is the same between the <^{'s and beween the 
^2's and that there is no interaction between the ^3's. 
However, for the present purpose, we consider, without 
loss of generality, the case where the interaction is nonva- 
nishing and identical between all vectors. One thus has 

V = Jll. Now, we use the decomposition of a rotation 
matrix $ 7 of 50(3) in terms of a four-component unit 

vector S 1 = (5^5f,5|,5|): 

^ = 2 (slS* - \8 k ^j +2 efeim 5 / 5i + 2 (V - ^ 4; . 

(19) 

In terms of the vector S 1 , the Hamihonian H18|) writes: 

H = -4J^ (S^.S^Y (20) 

(i, J) 

which is the Hamihonian for /owr-component nonfrus- 
trated spins with a particularity that each vector S 1 ap- 
pears quadratically. Therefore, the Hamihonian (|20|l is 
invariant under a global 0(4) group and under a local — 

gauge — Z 2 group that changes S 1 to — S 1 . It corre- 
sponds to the i?P 3 = 50(4)/(50(3) x Z 2 ) model. Note 
that, had we kept the microscopical coupling constants: 

V = diag(J, J, 0), the Hamihonian l|20|) would be supple- 
mented by terms breaking the 50(4) global symmetry 
and leaving untouched the Z2 local symmetry which is 
the important point for our purpose (see [78j for details). 
For three-component spins, an analogous Hamihonian — 
the RP 2 model — had been introduced by Maier and 
Saupe EH and by Lebwohl and Lasher [8(j to investigate 
the isotropic-nematic transition in liquid crystals. An 
extensive study of the RP 2 model, as well as a detailed 
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investigation of the role of vortices in this transition, has 
been performed by Lammert et al. j^, • These au- 
thors have shown, in particular, that these nontrivial 
topological configurations favor the first order character 
of the transition. In the case of four-component spins, no 
such detailed analysis has been performed. However, the 
RpN = SO{N)/(SO{N - 1) x Z 2 ) models that general- 
izes Hamiltonian l|2U[) to A^-component spins have been 
numerically studied in (8j| for 2 < N < 4. These systems 
have been shown to undergo a first order phase transi- 
tion. Since the only difference between the RP N and the 
0(N)/0(N - 1) — or, equivalently, SO(N)/SO(N - I) 

— models lies in their topological properties, one is nat- 
urally led to attribute the origin of the first order char- 
acter of the phase transition in the RP N models to the 
Z2 vortices. Finally, since the Hamiltonian l|18|) . rele- 
vant to STA, can be mapped onto the Hamiltonian i|20|) 

— up to the 0(4)-breaking terms — one can expect that 
the topological configurations also favor first order phase 
transitions in frustrated magnets in three dimensions. 

C. The XY case 

In the XY case, the Hamiltonian 10 is still invariant 
under a right 0(2) group, see Eq. (|II|I . while the left 
symmetry group becomes 0(2). In the low-temperature 
phase, the rotational symmetry is broken and, since the 
spins are constrained to be in a plane, the permutation 
symmetry between 4>\ and 4>2 is also broken. As a conse- 
quence, the symmetry breaking scheme is: 

G = 0(2) x 0(2) - H = 0(2) diag . (21) 

This symmetry-breaking scheme is usually referred to as 
50(2) x Z2 — > 1. The Z? d egre es of freedom are known 
as chirality variables |70l \72\ . |8J, |85| . 

In this case, there also exist topological defects since: 

tti(50(2)) =Z . (22) 

These defects are identical to those of the ferromagnetic 
XY model that drive the famous Berezinskii-Kosterlitz- 
Thouless transition in two dimensions |52ll53l |. However, 
in the frustrated case, they very likely interact non triv- 
ially with the chirality degrees of freedom which are crit- 
ical in d = 3 at the same temperature as the spin wave 
degrees of freedom. This is apparent from the fact that 
one observes a unique phase transition and not two dis- 
tinct Ising-like and XY-like transitions (8(|. As a conse- 
quence, one can expect, in the frustrated case, a physics 
different and probably more complicated than in the non- 
frustrated 0(2) model that undergoes a standard second 
order phase transition in three dimensions. 

D. Generalization 

For reasons that will become clear, we consider the gen- 
eralization of Hamiltonian Q to A-component spins. It 



is straightforward to extend the previous considerations 
to this case. One finds the symmetry breaking scheme: 

G = 0(N) x 0(2) -> H = 0(N - 2) x 0(2) diag . (23) 

In the following, we shall drop the "diag" index for 
simplicity. Note that the previous Heisenberg and XY 
cases are recovered trivially provided that we identify 
0(0) with the trivial group 1 and O(f) with Z2. 

We now give a review of the experimental and numeri- 
cal results for both the XY and Heisenberg systems. We 
will argue that a critical analysis of these results is cru- 
cial to understand that, up to now, the critical behavior 
of these systems is still unexplained. 

III. EXPERIMENTAL AND NUMERICAL 
SITUATIONS 

A. Preliminaries 

In this section, we analyze the experimental and nu- 
merical results relevant to the physics of frustrated mag- 
nets. Our aim is to show that these data are hardly com- 
patible with a second order phase transition since, in par- 
ticular, they show that frustrated magnets display scaling 
without universality. Moreover, we show that there are 
even some direct indications for weak first order behav- 
iors in these systems. We recall that a phase transition 
is said to be weakly of first order when, at the transition, 
the jump of the order parameter is small and the cor- 
relation length is large. Thus, scaling behaviors can be 
observed on a large range of temperatures so that these 
transitions look like second order phase transitions ex- 
cept very close to the critical temperature where scaling 
aborts. 

We emphasize that, by itself, the analysis of the ex- 
perimental and numerical results would not be sufficient 
to firmly conclude on the first order nature of the tran- 
sitions. It has to be seen as one of the pieces of the 
argumentation that, together with a theoretical analysis, 
will lead to a coherent picture of the critical physics of 
frustrated systems. 

To perform this analysis we need, in our discussion, 
to compare experimental results among themselves, as 
well as with numerical and theoretical calculations. Let 
us explain how we extract average values and error bars 
out of a set of experimental determinations of critical ex- 
ponents. In the experimental literature, only one error 
bar is quoted, which merges the systematic and statis- 
tical errors. Our first — minimal — hypothesis is that 
error bars have a purely statistical origin (no systematic 
error). Under this assumption, we can trivially compute 
the (weighted) average values of the exponents together 
with their error bars. This is the meaning of the numbers 
we give in the following when we deal with average values 
of the critical exponents. It is clear that this hypothesis 
is too simple to be realistic since the experimental sys- 
tematic errors cannot be neglected. Thus, the values we 
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compute, especially the error bars, should be taken with 
caution. We however show in the course of this article 
that our conclusions are robust to a possible underesti- 
mation of the error bars in our calculations, see Section 
IVB71 

Let us also notice that a possible source of error in the 
estimation of the critical exponents themselves could be 
the existence of corrections to scaling that could bias all 
the results. As we now argue, we can however expect 
that these effects are not dramatic. Let us consider the 
well-documented case of the ferromagnetic Ising model 
in d = 3. Most of the time corrections to scaling are not 
considered in the determination of the critical exponents 
and the associated error bars. When they are taken into 
account, they induce a tiny change in the critical expo- 
nents, i.e. at most of few percents (see for instance (S3 
and |S3 f° r a review) . It is therefore reasonable to think 
that neglecting corrections to scaling induce an error of 
few percents on the critical exponents while this proba- 
bly leads to largely underestimated error bars when those 
are announced to be of the order of 1% [208J. 

In the case of frustrated magnets, if we make the as- 
sumption that the corrections to scaling are comparable 
with those found in the ferromagnetic Ising model and 
bear in mind that the error bars quoted in the literature 
are of the order of 5 - 10% (see Tables III ITT1 and llVll . 
we are led to the conclusion that corrections to scaling 
are significant neither for the exponents nor for the error 
bars. 



B. The XY systems 

Let us first discuss the XY case since the experimen- 
tal situation is richer than in the Heisenberg case. Also, 
the symptoms of the existence of a problem in the in- 
terpretation of the results are clearer than in this latter 
case for reasons that shall be explained in this article and 
particularly in Section llXI 



1. The experimental situation 

Two classes of materials are supposed to be described 
by the Hamiltonian ©. The first one is made of ABX3 
hexagonal perovskites — where A is an alkali metal, B 
a transition metal and X a halogen atom — which are 
physical realizations of XY STA. The most studied ones 
are CsMnBr 3 , CsCuCl 3 , CsNiCl 3 and CsMnI 3 (see Hj 
for a review and [9(j for RbMnBr 3 . We have excluded 
this material since the measurement of its specific heat 
presents a shoulderlike anomaly near T c which renders 
the determination of a and (3 doubtful) . The second one 
is made of rare earth helimagnets: Ho, Dy, Tb. For most 
materials, the transitions are found continuous but not 
with the same critical exponents. For CsCuCl 3 , the tran- 
sition is found to be weakly of first order, i.e. with small 



discontinuities. The results are summarized in Tables ^ 
and ED 



Compound 


Ref. 


a 


P 


7 


V 


CsMnBr 3 


[86] 
[9TJ 
f92] 
[93] 
[94] 
[951 
[96] 
[96] 
[97] 


0.39(9) 
0.40(5) 
0.44(5) 


0.21(1) 
0.24(2) 
0.21(2) 
0.25(1) 
0.22(2) 


1.01(8) 
1.10(5) 


0.54(3) 
0.57(3) 


CsNiCl 3 


[98] 
[98] 
[99] 
[100] 


0.37(8) 
0.37(6) 
0.342(5) 


0.243 (5) 






CsMnI 3 


[98] 


0.34(6) 








CsCuCl 3 


[101] 
[90] 


0.35(5) 


0.23-0.25(2) 
1 st order 






TABLE I: The critical exponents of the XY STA. 


Compound 


Ref. 


a 


P 


7 


V 


Tb 


[102] 
[103] 
[104] 
[105] 


0.20(3) 


0.23(4) 
0.21(2) 




0.53 


Ho 


[106] 
[107J 
[95] 
[108] 
[108] 
[109] 
[110] 

[mj 

[112] 
[112] 
[113J 

[114] 


l st order 
0.27(2) 
0.10-0.22 


0.30(10) 
0.37(10) 
0.39(3) 
0.39(2) 
0.39(4) 
0.39(4) 
0.41(4) 

0.38(1) 


1.24(15) 
1.14(10) 


0.54(4) 
0.57(4) 


Dy 


[115] 
[116J 
[110J 
[109] 
[113J 
[117] 


0.24(2) 


0.335(10) 

0.38(2) 
0.39(1) 


1.05(7) 


0.57(5) 



TABLE II: The critical exponents of the XY helimagnets. 



We highlight four striking characteristics |118| of these 
data. Their consequences for the physics of frustrated 
magnets will be discussed in more details in the following. 

i) There are two groups of incompatible exponents. In 
the following discussion, we mainly use the exponent f3 
to analyze the results since, as seen in Tables IJ and |Hl it 
is by far the most precisely measured exponent. Clearly, 
there are two groups of materials, each of which being 
characterized by a set of exponents, (3 in particular. 

In the first one — that we call group 1 — made up of: 
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group 1 : CsMnBr 3 , CsNiCl 3 , CsMnI 3 , Tb (24) 
one has: 

(3 ~ 0.237(4). (25) 

Note that, as far as we know, there is no determination 
of the exponent (3 for CsMnI 3 that, being given its com- 
position, has been included in the group 1 of materials. 
Anyway, our conclusions are not affected by this fact. 
In the second — group 2 — made up of 

group 2 : Ho, Dy (26) 

one has: 

/3~ 0.389(7). (27) 

These exponents are clearly incompatible. Actually, we 
find for the average exponents of CsMnBr 3 alone — the 
most and best studied material of group 1 — : 

(3 = 0.228(6), v = 0.555(21), a = 0.416(33), 7 = 1.075(42) . 

(28) 

If we consider all the materials of group 1 (except Tb for 
which the results are not fully under control, but perhaps 
(3) we find: 

(3 = 0.237(4), v = 0.555(21), a = 0.344(5), 7 = 1.075(42) . 

(29) 

For materials of group 2 (Ho and Dy) we find: 

f3 = 0.389(7), v = 0.558(25), 7 = 1.10(5) . (30) 

We do not give a value for a which is poorly determined. 

Let us indicate that the exponents vary much from 
compound to compound in group 1. Although less ac- 
curately determined than (3, a is only marginally com- 
patible between CsNiCl 3 and CsMnBr 3 . Note moreover 
that, even for the same material, the data are not fully 
compatible among themselves: (3 in CsMnBr 3 shows a 
somewhat too large dispersion. 

ii) The anomalous dimension rj is negative for group 1 
which is impossible. If we assume that the transition is of 
second order for group 1, we can use the scaling relations 
to compute 77. In particular, the precise determination 
of (3 allows to use 77 = 2(3 /v — 1 to determine rather 
accurately rj. The exponent v itself can be obtained di- 
rectly from the experiments or deduced using the scaling 
relation: 

u=(2-a)/3. (31) 

The large number of experiments devoted to the determi- 
nation of a allows a precise determination of v. By using 
the scaling relation Eq. I|31|) . we find v = 0.528(11) if 
we consider the experimental results for CsMnBr 3 alone 



and v = 0.552(2) if we consider CsMnBr 3 , CsNiCl 3 and 
CsMnI 3 . By using the relation 77 = 2/3/f — 1 together 
with Eq. H31|) or the relation 77 = 6/3/(2 — a) — 1 and 
by considering the data of CsMnBr 3 alone or the data of 
the materials of group 1 (except Tb for which it is not 
sure that the data are reliable) we can obtain four de- 
terminations of 77. In the four cases, we find 77 negative 
by at least 4.1 standard deviations and the probability 
to find it positive always less than 10 -5 . In fact, the 
most precise determination is obtained by combining all 
the data of group 1, Eq. and by using the relation 

rj = 6/3/(2-a)-l. In this case, we obtain r\ = -0.141(14) 
and thus a (almost) vanishing probability to find it posi- 
tive. Note also that, although (3 and v are less accurately 
known in Tb — for which experiments are anyway deli- 
cate — , rj is also found negative. 

However, we stress that rj cannot be negative in a true 
second order phase transition. This is a general result, 
based on first principles of field theory, that 77 is always 
positive if the theory des cribi ng the transition is a uni- 
tary GLW c/> 4 -like model EH 

as it is the case here (see 

Appendix [XJ . 

Hi) For group 2, the scaling relation 7 + 2(3 — 3v = is 
violated. From Eq. 113011 it is possible to check the scaling 
relations. We find 7 + 2(3 - iu = 0.202(92) and thus a 
violation by 2.2 standard deviations. 

iv) CsCuCl^ undergoes a weak first order phase tran- 
sition. Until recently, CsCuCl 3 was believed to undergo 
a second order phase transition with exponents compati- 
ble with those of group 1, see Table [I] It has been finally 
found to display a weak first order phase transition [90| . 



2. The numerical situation 

Monte Carlo simulations have been performed on five 
different kinds of XY syst ems. The first one is the STA 
itself EM EM EH EH EM EM Ell- The second 
model is the STAR (where R stands for rigidity) which 
consists in a STA for which the local rigidity constraint 

— Eq. (0J — has b een imposed on each plaquette at all 
temperatures |l27f . The third model is the Stiefel V2 2 
model whose Hamiltonian is given by Eq. © EME3- 
This is a hard spin, discretized version of the NLct model 
relevant to frustrated magnets. Note that, for this last 
model, the triangular structure is irrelevant since the in- 
teraction is ferromagnetic; a cubic lattice can be chosen. 
Also a soft spin, discretized version o f the GLW model 
has recently been studied by Itakura |l26j who also re- 
studied the STA model for large sizes. Finally, a heli- 
magnetic system defined on a body-centered-tetragonal 
(BCT) lattice — the "BCT model" — has been investi- 
gated [l29j . 

Here, we emphasize that the local rigidity constraint 

— Eq. (0J — as well as the manipulations that lead to 
the STAR, Stiefel V 2 , 2 , GLW and BCT models only af- 
fect the massive — noncritical — modes. Thus, all the 
STA, STAR, Stiefel V 2 , 2 , GLW and BCT models have the 
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same critical modes, the same symmetries and the same 
order parameter. Therefore, one could expect a common 
critical behavior for all these systems. 

Let us comment the results of the simulations given 
in Table UTTl Note that, due to the its novel char acter, 
we shall comment the recent work of Itakura |126| sepa- 
rately. 



System 


Ref. 


a 





7 


V 


V 


STA 


[120, 123] 

[124] 

[125] 


0.34(6) 

0.46(10) 

0.43(10) 


0.253(10) 
0.24(2) 


1.13(5) 
1.03(4) 


0.54(2) 
0.50(1) 
0.48(2) 


-0.09(8) 
-0.06(4) 


STA 


[126] 


1 st order 


STAR 


[127] 


1 st order 


V 2 ,2 


[127] 


1 st order 


BCT 


[129] 


1 st order 


GLW 


[126] 


1 st order 



TABLE III: Monte Carlo critical exponents of XY systems. 
Note that the exponent n is computed from y/u = 2 — n. 

i) For STA, scaling laws are found with exponents com- 
patible with those of group 1. Let us however notice that, 
similarly to what happens for the materials of group 1 
there exists, in the numerical simulations of STA, a rather 
large dispersion of the results. For instance, the two ex- 
treme values of v differ by 2.1 standard deviations. 

Let us make two other remarks. First, the good agree- 
ment between the numerical results for STA and the ex- 
perimental ones for materials of group 1 has been re- 
peatedly interpreted in the literature as a proof of the 
existence of a second order transition and even as an ev- 
idence of the existence of the chiral fixed point of the 
GLW model [jjj . We emphasize here that the fact that a 
Monte Carlo simulation reproduces experimental results 
only means that the Hamiltonian of the simulated system 
is a good approximation of the microscopic Hamiltonian 
describing the physics of real materials. However, this 
neither explains nor proves anything else — and certainly 
not the existence of a second order phase transition — 
since Monte Carlo simulations suffer from problems anal- 
ogous to those encountered in experiments: a weakly first 
order phase transition is very difficult to identify and to 
distinguish from a second order one. 

Let us now come to our second remark. In a beauti- 
ful experiment, Plakhty et al. [8(| have meas ured the 
so-called chiral critical exponents (3 C and y c |l2dj in 
CsMnBr3 . They have foun d va lues compatible with those 
found numerically in STA |l20j| . Let us emphasize, again, 
that this agreement simply means that the parameters 
characterizing the numerical simulations are not too far 
from those associated to the experiments. By no means it 
implies — or gives a new indication of the existence of — 
a second order transit ion. Let us notice that j3 c has also 
been measured in Ho |ll4| . The value found completely 
disagrees with the result found in STA and in CsMnBr3. 

ii) The anom alous dimension r\ is negative for STA. As 
shown in |l27j , rj is found negative using the two scaling 



relations r) = 2(3 Jv — 1 and 77 = 2 — 7/1/ for the two 
simulations where these calculations can be performed. 

Hi) The simulations performed on STAR, V%fl and 
BCT models give first order transitions. Therefore, the 
modifications in the microscopic details which change 
STA into STAR, V 2 . 2 and BCT affect drastically the scal- 
ing behavior. 

iv) In a remarkable work, Itakura has recently per- 
formed Monte Carlo and Monte Carlo RG approaches 
of the STA and its GLW mode l version that has led to 
a clear first order behavior X12&I . Itakura has performed 
standard Monte Carlo simulations of the STA involving 
sizes up to 126 x 144 x 126 leading to clear first order tran- 
sitions. In particular, for these lattice sizes, the double- 
peak of the probability distribution of the energy at the 
transition is clearly identified. Itakura has also used an 
improved Monte Carlo RG simulation of the STA and 
its GLW model version. One advantage of this approach 
compared with previous RG Monte Carlo studies is that 
it allows to reach the asymptotic critical behavior using 
systems of moderately large lattice sizes. Within this 
approach, Itakura has found several evidences for a first 
order behavior with, notably, a runaway behavior of the 
RG flow and the absence of any nontrivial fixed point. 



3. Summary 

We now summarize the results of our analysis of both 
experiments and numerical simulations for XY frustrated 
magnets. 

1) Scaling laws are found in STA and helimagnetic ma- 
terials on a rather wide range of temperature. Thi s is 
also the case within all — but an important one |l26j 
numerical simulations of the STA. 

2) There are two groups of systems that differ by their 
critical exponents. The first one includes the group 1 
of materials and the numerical STA model. The second 
one corresponds to the group 2 of materials. One also 
observes variations of critical exponents inside a given 
group of exponents. 

3) The anomalous dimension r) is negative for the ma- 
terials of group 1 and for the numerical STA model. This 
is very significant from the experimental results, less from 
the numerical ones. 

4) For group 2, the scaling relations are violated by 2.2 
standard deviations. 

5) CsCuCl3 is found to undergo a weak first order tran- 
sition. 

6) STAR, V2,2 an d BCT models undergo strong first 
order transitions. 

7) Recent Monte Carlo and Monte Carlo RG ap- 
proaches of STA and the soft spin discretized version 
of the GLW model give clear indications of first order 
behaviors. 
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4- Conclusion: five possible scenarios 

Let us now propose five possible scenarios to explain 
the phenomenology of XY frustrated systems. 
Scenario I. 

This scenario is — together with the second one — the 
most often invoked: the critical behavior of frustrated 
magnets, when they display scaling, is controlled by a 
unique fixed point of the RG flow whi ch is associated 
to a new universality class [63, EH, I122L Il23t Il30j . 
Although, from point 1) above, XY frustrated magnets 
appear to display rather generic scaling behaviors, the 
examination of the experimental and numerical data pro- 
vides clear indications against this first scenario. Indeed, 
from point 2), there is a manifest lack of universality in 
the scaling behavior of frustrated magnets. Also several 
points, from 3) to 7), strongly militate in favor of first 
order behaviors. 

Scenario II. 

In the second scenario, the two sets of exponents cor- 
responding to groups 1 and 2 are, in fact, associated with 
two true second order phase transitions from which result 
two distinct universality classes. This scenario is ruled 
out by the fact — see point 3) — that the anomalous di- 
mension 7/ is negative for group 1 and for the numerical 
STA model. Thus, provided i) the quoted error bars in 
the literature are reliable, ii) our hypothesis of a purely 
statistical origin of the errors does not completely bias 
our analysis and in) corrections to scaling do not alter 
drastically all the results, we are led to the conclusion 
that the behavior of the materials of group 1 and of the 
numerical STA model cannot be explained by the exis- 
tence of a fixed point in the GLW model. In the simplest 
hypothesis, these systems must undergo first order phase 
transitions. This last hypothesis seems to be confirmed 
by several other facts. Firstly, CsCuCl3, whose expo- 
nents are close to those of group 1 has been finally found 
to undergo first order phase transitions, see point 5) of 
the summary. Secondly, point 6), numerical models very 
close to STA — STAR, V 2 ,2 and BCT — also undergo 
first order phase transitions. Finally the hypothesis of a 
first order phase transition for STA itself is corroborated 
by the fact, point 7), that recent Monte Carlo and Monte 
Carlo RG simulations of this system predict a first order 
phase transition |126| . 

Scenario III. 

In the third scenario, materials of group 2 undergo a 
second order phase transition — 77 is found positive there 
— while those of group 1 as well as the numerical STA 
model all undergo weakly first order phase transitions. 
Within this scenario, the critical exponents of materials 
of group 1 should be considered as effective or pseudo- 
critical exponents, characterizing the pseudo-scaling ob- 
served, valid for temperatures far enough from the critical 
temperature. There is no direct and definitive argument 
against this scenario. Of course, violation of the scaling 
relations for materials of group 2, point 4), makes doubt- 
ful a second order behavior. However, this violation is 



too small to definitely reject it. Actually, the drawback 
with this third scenario is its lack of naturalness. Indeed, 
it implies a very specific fine-tuning of the microscopical 
coupling constants — i.e. of the initial conditions of the 
RG flow — for materials of group 1. Their representative 
points in the coupling constant space must lie outside the 
basin of attraction of the fixed point governing the criti- 
cal behavior of materials of group 2 but very close to its 
border so that the transitions are weakly of first order. 
Scenario IV. 

In the fourth scenario, all frustrated magnets undergo 
first order phase transitions that almost generically ap- 
pear to be weak or very weak and are characterized 
by pseudo-scaling and pseudo-critical exponents. This 
fourth scenario, compared with the third one, could thus 
seem even more unnatural. This is true, but only within 
the usual explanation of weak first order phase transi- 
tions where the weakness of the first order transition is 
obtained by fine-tuning of parameters. Actually, we shall 
provide arguments in favor of the present scenario and 
shall show that the genericity of pseudo-scaling has, in 
fact, a natural explanation relying neither on the exis- 
tence of a fixed point nor on a fine-tuning of parameters. 

Scenario V. 

Finally, one can imagine several variants of these sce- 
narios. For instance, we have adopted the standard posi- 
tion that consists in associating a unique set of critical ex- 
ponents t o a fixed point. On the contrary, Calabrese et al. 
|l3ll Il32fl have suggested that a unique fixed point could 
lead to a whole spectrum of critical exponents. This sce- 
nario, which would explain the occurence of a spreading 
of critical exponents in the experimental and numerical 
contexts, will be discussed in details in the following. 

We now review the experimental and numerical results 
obtained for the Heisenberg systems. 

C. The Heisenberg systems 

1. The experimental situation 

Contrarily to the XY case, there is no Heisenberg he- 
limagnets (see however |67|). Therefore there remain, a 
priori, only the Heisenberg STA materials. In fact, the 
A/B phase transition of He3 can be described by the 
same GLW Hamiltonian as the Heisenberg STA H,|g|. It 
is thus a candidate. Unfortunately, the narrowness of the 
critical region of this transition does not allow a reliable 
study of the critical behavior of this system and there are 
no available data about it. 

Three classes of Heisenberg STA materials have 
been studied. First, systems like VCI2, VB1-2, 
Cu(HCOO)22CO(ND 2 )22D 2 an EefSaCNfCaHsJakCl 
which are generically quasi-XY except in a particular 
range of temperature where their anisotropics are ir- 
relevant. Second, those which become isotropic thanks 
to a magnetic field that exactly counterbalances the 
anisotropies. This is the case of CsNiCLj and CsMnl3 
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at their multicritical point. Finally, those which be- 
come isotropic because they have been prepared in a fine- 
tuned stcechiometry such that the Ising-like and XY-like 
anisotropics cancel each other to form an isotropic mate- 
rial. This is the case of CsMn(Bro.igIo.8i)3- 

Let us comment the experimental results summarized 
in Table HVI 



Compound 


Ref. 


a 





7 


V 


VC1 2 


[133] 




0.20(2) 


1.05(3) 


0.62(5) 


VBr 2 


[134] 


0.30(5) 








A 


1135] 




0.22(2) 






B 


[87, 136, 137] 
[138] 


0.244(5) 


0.24(1) 


1.16(3) 




CsNiCla 


198, 139] 

[99] 

[100] 


0.25(8) 
0.23(4) 


0.28(3) 






CsMnI 3 


[98] 


0.28(6) 








C 


1140] 
[141] 
[142] 


0.23(7) 


0.29(1) 
0.28(2) 


[0.75(4)] 


[0.42(3)] 



TABLE IV: The critical exponents of the Heisen- 
berg STA. The abbreviations A, B and C stand 
for Cu(HCOO)22CO(ND 2 )22D 2 0, Fe[S 2 CN(C 2 H 5 )2]2Cl and 
CsMn(Bro. 1910.81)3 respectively. The data in brackets are sus- 
pected to be incorrect. They are given for completeness. 

i) As in the XY case, the Heisenberg materials fall into 
two groups. The group 1, made up of: 

group 1 : Cu(HCOO) 2 2CO(ND 2 ) 2 2D 2 0, 
Fe[S 2 CN(C 2 H 5 ) 2 ] 2 Cl, 
VCl 2 ,VBr 2 (32) 

is characterized by: 

(3 = 0.230(8) (33) 
while for group 2, made up of: 

CsMCl 3 , CsMnI 3 , CsMn(Br Q ig , I . 81 ) 3 (34) 

one finds: 

13 = 0.287(8) . (35) 

Note that, strictly speaking, the values of (3 for VBr 2 and 
for CsMnl3 are not known and, thus, our classification is 
somewhat improper. It seems however logical to suppose 
that VBr 2 is close to VC1 2 and CsMnI 3 close to CsNiCl 3 . 
Anyway, it will be clear in the following that our analysis 
is almost insensitive to this point. 

For group 1, the average values of the critical expo- 
nents are given by: 

(3 = 0.230(8), a = 0.272(35), v = 0.62(5), 7 = 1.105(21) . 

(36) 

A very severe difficulty in the study of the materials of 
group 1 is their two-dimensional character and Ising-like 



anisotropics. The temperature range where the systems 
behave effectively as three-dimensional Heisenberg sys- 
tems is narrow. This is the case of VC1 2 where this range 
is less than two decades and where, closer to the criti- 
cal temperature, the system becomes Ising-like. For this 
group of mater ials the exponent (3 is very small and the 
authors of |l33fl have noticed that such small values have 
also been found in materials where dimensional cross-over 
is suspected. Thus, it is not clear whether the whole 
set of results really corresponds to a three-dimensional 
Heisenberg STA. 

For group 2, the experimental situation seems to be 
better under control. The average values of the critical 
exponents are given by: 



(3 = 0.287(9), a = 0.243(3), v = 0.585(9), 7 



1.181(33) 
(37) 



where the scaling relations have been used to compute v 
and 7. Note that the values of v and 7 thus obtained dif- 
fer significantly from those of CsMn(Bro.igIo.8i)3 whose 
critical behavior has been claimed to be perturbed by 
disorder (see however [209] ). 

ii) For group 1, the anomalous dimension r\ is signifi- 
cantly negative. Using the two exponents that have been 
measured at least twice in group 1 — [3 and 7 — we can 
compute the anomalous dimension from the scaling re- 
lation X] = (4/3 - 7)/(2/3 + 7). We find n = -0.118(25) 
which is thus negative by 4.8 standard deviations. 

Hi) For group 2, the anomalous dimension r] is 
marginally negative. Using the critical exponents given 
in Eq. (|37|l . one obtains, for the anomalous dimension: 
i] = —0.018(33). Thus r\ is found negative but not signif- 
icantly, contrarily to what happens in group 1. 

iv) For group 1, the scaling relations 7 + 2/3 — 2 + a = 
= 2/3 + 7 — 3z^ are violated. Indeed, 7 + 2/3— 2 + a = 
-0.135(56) and 2/3 + 7 - iv = -0.29(15). Of course, 
none of these violations is completely significant in itself 
because of the lack of experimental data. However, since 
they are both independently violated it remains only a 
very small probability that the scaling relations are ac- 
tually satisfied. 



2. The numerical situation 

In the Heisenberg case, as in the XY case, five different 
kinds of systems: STA, STAR, Stiefel (V 3j2 in this case), 
BCT and GLW models have been studied. The results 
of the simulations are given in Table Ivl 

Let us c omm ent them. Again, we put aside the work 
of Itakura |l26j . 

i) For the STA, scaling laws are found with an exponent 
(3 close to that of group 2. The average values for the 
exponents of STA are: 



[3 = 0.288(6), 7 = 1.185(3), = 0.587(5) 



(38) 



/3 is thus extremely close to the experimental value of 
group 2 while v and 7 are extremely close to the ex- 
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System 


Ref. 




P 


7 


V 


V 


STA 


[120, 130] 
443] 
444] 
445] 
1140 1 


0.240(80) 
0.242(24) 
0.245(27) 
0.230(30) 


0.300(20) 
0.285(11) 
0.289(15) 
0.280(15) 


1.170(70) 

1.185(3) 

1.176(26) 


0.590(20) 

0.586(8) 

0.585(9) 

0.590(10) 

0.589(7) 


0.020(180) 
-0.033(19) 
-0.011(14) 

0.000(40) 


nrri* 
BlA 


|12t>| 


1 st order 


b 1AK, 


r i a ti 
|i4rj 


0.488(30) 


0.221(9) 


1.074(29) 


0.504(10) 


-0.131(13) 


V'i 9 
v o,Z 


[147] 
[128] 


0.479(24) 
0.460(30) 


0.193(4) 


1.136(23) 
1.100(100) 


0.507(8) 
0.515(10) 


-0.240(10) 
-0.100(50) 


V 3 ,2 


[126] 


1 st order 


BCT 


[148] 


0.287(30) 1 0.247(10) 1 1.217(32) 1 0.571(10) |-0.131(18) 


GLW 


[126] 


1 st order 



TAB LE V: Monte Carlo critical exponents of the Heisenberg systems. r\ is computed by j/v = 2 — n and, apart in | l20l ] and 
128], a is computed by Sv = 2 — a. 



perimental values deduced from the scaling relations, 
Eq. The scaling relation 7 + 2(3 — 2>v = is very 

well verified since 7 + 2/3 - 3v = 1(T 4 ± 6.10 -2 . 

ii) For the STA, rj is negative. Using the values of 
(3/v and 7/z/ obtained directly in the simulations, one 
can compute the average value of r\: —0.0182(89). The 
probability of it to be positive is 0.02 and is thus small 
although not vanishing. 

Hi) For the STAR, VJj^ and BCT models, the values of 
(3 are all incompatible with that of STA (three standard 
deviations at least) and are all incompatible among each 
others. This has been interpreted as an i ndication of 
very weak first order phase transitions |l47j ]. This is to 
be compared with the XY case, where the transitions for 
STAR and the V2.2 model are strongly of first order. 

iv) For the BCT, STAR and models, r\ is always 
found significantly negative, see TablefVlwhere r\ has been 
calculated from j/v. 

v) The Monte Carlo and Monte Carlo RG approaches 
of the STA, V-i,2 and GLW model pe rform ed by Rakura 
has led to clear first order behaviors |l26| . For Heisen- 
berg STA, contrarily to the XY case, even for the largest 
lattice sizes — 84 x 96 x 84 — the double-peak of the 
probability distribution of the energy is not observed. 
However, the model displays a clear double-peak. 
Moreover, for the STA and the model, the RG flow 
clearly does not exhibit any fixed point. Instead, a run- 
away of the RG flow toward the region of instability is 
found which indicates first order transitions. The transi- 
tions are thus — weakly — of first order. The transition 
is also weaker of first order for Heisenberg than for XY 
spins. 



3. Summary 

We now summarize the experimental and numerical 
situations for frustrated magnets with Heisenberg spins. 
Here, the experimental situation is much poorer than in 
the XY case and is still unclear on many aspects. On the 



contrary, the numerical results are numerous and more 
precise than in the XY case. 

1) Scaling laws are found in STA materials on a rather 
wide range of temperatures as well as in all Monte Carlo 
simulations — apart that based on Monte Carlo RG. 

2) There are two groups of materials that do not have 
the same exponents. The exponent j3 of the numerical 
STA model agrees very well with that of group 2. 

3) The anomalous dimension is manifestly negative for 
group 1 and marginally negative for group 2. For the 
numerical STA model, rj is found negative although not 
completely significantly. For STAR, and BCT, r\ is 
found significantly negative. 

4) For group 1, the scaling relations are violated. 

5) STAR, V3,2 and BCT exhibit scaling behaviors with- 
out universality. Also, the results are incompatible with 
that of the numerical STA model. 

6) A Monte Carlo RG approach of the STA, V 3y2 and 
GLW models has led to clear first order behaviors. 



4- Conclusion 

Let us now draw some conclusions about the Heisen- 
berg case. The experimental and numerical data reveal 
the same problems as those encountered in the XY case: 
the different materials split into two groups, the anoma- 
lous dimension is found negative in many materials and 
in most numerical simulations, the scaling relations are 
violated in some materials and there is no universality 
in the exponents found in the simulations. The same 
kind of conclusion as in the XY case follows (see Section 
IIII B 4p : the first scenario, that of an explanation based 
on the existence of a unique fixed point appears unlikely. 
There are also signs of first order behaviors but less sig- 
nificantly than in the XY case. Thus, at this stage, it 
is impossible from the experimental and numerical data 
alone to discriminate between the different scenarios II, 
III, IV and V. It is therefore important to gain insight 
from the theoretical side. 
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Before discussing this, let us mention another interest- 
ing numerical result. 



D. The N = 6 STA 

Let us quote a simulation of the STA with six- 
component spins that has been performed by Loison et 
al. • The results are given in Table lvTl Six-component 
spins were chosen since it was expected that the tran- 
sition was of second order. Loison et al. have clearly 
identified scaling laws at the transition with a positive 
anomalous dimension. Let us emphasize that, even if the 
transition is actually of first orde r, as suggested by the 
recent results of Calabrese et al. [l32j . it should be ex- 
tremely weakly of first order — see the following. Thus, 
scaling laws should hold for all temperatures but those 
very close to T c . In this respect, the exponents for N = 6 
are therefore very trustable so that reproducing them is 
a challenge for the theoretical approaches. 



I System | a \(3 1 7 | v | 

\N = 6 STA |-0.100(33) 1 0.359(14) 1 1.383(36) 1 0.700(11) | 

TABLE VI: Monte Carlo critical exponents for six-component 
spins in the STA system (2^|. 

Note that using the results of Loison et al. and the relation 
?7 = 2 - j/u, one finds n = 0.025(20). 



IV. A BRIEF CHRONOLOGICAL SURVEY OF 
THE THEORETICAL APPROACHES 

Let us briefly review the most important theoretical 
developments concerning this subject. 

The first microscopic derivation and RG study — at 
one- and two-loop order in d = 4 — e — of the effec- 
tive GLW model relevant for the STA — see below - 
was performed for He3 byJones et al. in 1976 @ and 
by Bailin et al. in 1977 [9j. The model was re-derived 
and re-studied in the context of hclimagnets (for gen- 
eral N) by several groups including Bak et al. (1976) 
p9j . Garel and Pfeuty (1976) H3 and Barak and Walker 
(1982) |150| . It was established at that time that, around 
d = 4, the transitions for Heisenberg spin systems had to 
be of first order. More precisely, these authors found that 
there exists a critical value N c (d) of the number N of spin 
components above which the transition is of second order 
and below which it is of first order. They found : 

N c {d = 4 - e) = 4(3 + V6) - 4 ^3 + e + 0(e 2 ) 

~21.8-23.4e + 0(e 2 ) 

(39) 

with e = 4— d. A first large- N expansion was also studied, 
in particular for d — 3, by Bailin et al. in 1977 0. 



A group theoretical derivation of the GLW model rel- 
evant to the XY STA was performed by Yosefin and 
Domany (6^| in 1985. They found the same Hamilto- 
nian as f or he l imagnets . Between 1985 and 1988 Kawa- 
mura |7Ctll22L ll23. 130j have performed this analysis for 
iV-component STA. He has shown that the Hamiltonian 
is the same as for He3 or helimagnets, the RG analysis 
giving obviously the same results. This author has also 
extrapolated the two-loop result for N c (d) of Eq. i|3~9"|) in 
d = 3 and found N c (d — 3) < 2. This led him to con- 
jecture the existence of a second order phase transition 
for frustrated magnets associated with a new universality 
class. However, as well-known, this direct extrapolation 
cannot be reliable since it is notorious that the pertur- 
bative series must be resummed. 

In 1988, Dombre and Read |6j| derived, in the quan- 
tum case, the Nonlinear Sigma (NLct) model relevant to 
frustrated magnets. In 1990, Azaria et al. studied the 
classical version of this NLct model around d = 2. They 
found a fixed point of the RG flow in a two- loop calcu- 
lation for any N > 3 [l5l| . For N = 3, they found the 
phenomenon of enlarged symmetry: at the fixed point 
the symmetry becomes SO (3) x 50(3) ~ 50(4) instead 
of 50(3) x 50(2). Thus, their conclusion was that, if 
the transition is of second order, it is characterized by 
0(4)/0(3) critical exponents — at least for v. Another 
possibility proposed by these authors was that the tran- 
sition could be also mean-field tricritical or of first order. 
However, none of the experimental or numerical results 
are compatible with the 0(4)/0(3) or mean-field tricrit- 
ical exponents. Note finally that these authors supposed 
that, if tricritical, the behavior at the transition should 
be mean- field tricritical in d = 3, something which is 
mandatory only for 0(N)/0(N — 1) models, but not for 
more complex models. 

The first nonperturbative approach to frustrated mag- 
nets was performed by Zumbach in 1993 [64L l65l l66| . 
He wrote down the NPRG equations for the GLW mod- 
els suited to the description of these systems. He 
studied them within the Local Potential Approximation 
(LPA) of the Wilson-Polchinski equation — analogous to 
the Wegner-Houghton approximation [30j — and found 
N c (d = 3) ~ 4.7. Since he found no fixed point for N — 2 
and N = 3 he claimed that the transition is of first order 
in these cases. In the case N = 3, there is a minimum 
in the RG flow, a pseudo-fixed point, that fakes a true 
fixed point (see below for details). The transition was 
thus conjectured to be weakly of first order with pseudo- 
scaling characterized by pseudo-critical exponents. Note 
that, within the LPA, all derivative terms in the Hamil- 
tonian are neglected so that the anomalous dimension is 
vanishing. This has two important consequences. First, 
the pseudo-critical exponents found by Zumbach were 
not very reliable and thus difficult to compare with the 
experimental and numerical results. Second, this ap- 
proach neglects terms — the so-called current-term (see 
below) — that are fundamental within the perturbative 
approach of the NLcr model performed around two di- 
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mensions. Thus, within Zumbach's approach, it was not 
possible to match with these results. Finally note that, 
in the N = 2 case, no minimum in the RG flow was found 
and, thus, no pseudo-critical exponent was obtained, in 
contradiction with the scaling behaviors observed in the 
experimental and numerical contexts. 

Then, three-loop calculations have been performed by 
Antonenko et al. in 1994 and 1995 on the GLW model. 
In d = 3 this has led, afte r Pade-Borel resummation, 
to N c (d = 3) = 3.91 In d = 4 - e, they have 

determined the three- loop contribution — 7.1 e 2 — to 
N c (d = 4 - e) see Eq. (J3SJ). This has led to N c (d = 
3) = 3.39 |l53j . These authors have mentioned that, 
contrarily to the O(N) models, their three- loop results 
were not well converged. 

In 1996, Jolicceur and David studied a generalization 
of the Stiefel model that involves N vectors with N com- 
ponents |l54j . They showed within a mean field approxi- 
mation and a one-loop calculation performed in d = 2 + e 
that a first order line should appear in a nontrivial dimen- 
sion above two. It should isolate the chiral fixed point 
in the metastability region in such a way that this point 
should no longer play any role. Above this dimension, 
the transition should therefore be of first order. 

Then, in 2000-2002, using the technique of the ef- 
fective average action, including derivative terms, the 
present authors performed a nonperturbative study of 
frustrated magn ets for any dimension between two and 
four |6^l63llll8j . They recovered all known perturbative 
results at one-loop in two and four dimensions as well as 
for N — > co. They determined N c (d) for all d and found 
N c (d = 3) = 5.1. Accordingly, for N = 6, they found a 
second order phase transition. Their exponents were in 
very good agreement with those found numerically. For 
-/V = 3 |63j . they recovered Zumbach's results — pres- 
ence of a minimum in the RG flow — and improved his 
approach: they found pseudo-critical exponents in good 
agreement with some experimental realizations of frus- 
trated magnets. However, regarding the spreading of the 
experimental and numerical data, the recourse to a min- 
imum, leading to a unique set of pseudo-critical expo- 
nents, was clearly not th e end of the story. During the 
study of the N — 2 case |118| . the present authors real- 
ized that the property of pseudo-scaling and even more, 
generic pseudo-scaling, does not rely on the concept of 
minimum of the flow. Pseudo-scaling appears as a con- 
sequence of the existence of a whole region in the flow 
diagram in which the flow is slow. This allowed them to 
account for the nonuniversal scaling that occurs in XY 
as well as in Heisenberg frustrated magnets. The present 
article accounts for these last developments. 

In 2001, Pelissetto et al. |l55l | derived the six- loop se- 
ries for the GLW model. They used sophisticated re- 
summation methods in order to find the fixed points and 
to determine the critical exponents of the model. For 
N > 7, they found a fixed point of the same nature as 
the one obtained at large- TV and in the 4 — e expansion. 
Thus, a second order phase transition is expected in this 



case. For 5 < N < 7, they considered that their re- 
summed series were not well converged, the number of 
fixed points depending strongly of the number of loops 
considered. This led them to interpret this result as an 
indication that N c (d = 3) ~ 6. Finally and surprisingly, 
for N < 5 and, in particular, for the physically relevant 
cases N = 2 and N = 3, they found stable fixed points. 
Thus, a second order phase transition was also predicted 
in these cases. However, the critical exponents found 
were far from all experimental and numerical data (see 
the following). Moreover, regarding again the spreading 
of these data, an interpretation in terms of a unique set 
of exponents was cl early insufficient. 

In another work 156], assuming that N c (d = 2) = 2, 
Pelissetto et al. have reformulated the three-loop ver- 
sion of the series of Eq. I|39|l — see below — to make 
it compatible with this last guess. The series seemed to 
have better convergence properties — see however be- 
low — and allowed Pelissetto et al. to compute N c (d). 
They found N c (d = 3) = 5.3, in good agreement with the 
value — NJd = 3) = 5.1 — obtained from the NPRG 
approach |63Lll57l |. 

Recent re-i nvestigat ions of the five and six- loop pertur- 
bative series 

ESSIES 

have led Calabrese et al. (see also 
|158( 1 to conjecture that the fixed point found by Pelis- 
setto et al. — that corresponds to a focus fixed point 
- could explain the existence of the spreading of crit- 
ical exponents encountered in frustrated magnets. In- 
deed, they observed that, due to the specific structure 
of the fixed point, the critical exponents display strong 
variations along the RG trajectories that could explain 
the lack of universality observed experimentally and nu- 
merically. They have also given estimates of the critical 
number of spin components for which there is a change 
of the order of the phase transition. They have found 
that there is a first order phase transition in the whole 
domain 5.7(3) < N < 6.4(4) and a second order phase 
transition for the other values of N and, in particular, 
for N = 2 and N = 3. 

Finally, a very recent computation of the five-loop (3 
function of the GLW model in a 4 — e expansion has lead 
to a nov el estimation of N c (d). Calabrese and Parruc- 
cini [Hit have found the value N c (d = 3) = 6.1(6) which 
is compatible with the value N c (d — 3) = 6.4(4) found 
within th e six - loop computation performed in three di- 
mensions |l3llll32| . 

Since several aspects of the recent perturbative and 
nonperturbative approaches differ, in particular in their 
interpretations of the origin of the nonuniversal scaling 
found in frustrated magnets, we postpone the detailed 
discussion of these last developments of both methods to 
the following sections. 



V. THE PERTURBATIVE SITUATION 

Let us discuss in more details the perturbative ap- 
proaches that have been used to investigate frustrated 
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magnets. 

There are essentially two different methods to ana- 
lyze the critical behavior of the system described by the 
Hamiltonian @. They correspond to two different meth- 
ods to deal with the constraints obeyed by the micro- 
scopic degrees of freedom, Eq. © . They lead to the NLcr 
and GLW models that have been both perturbatively an- 
alyzed around their respective critical dimension as well 
as, for the GLW model, directly in three dimensions. Let 
us review the results of these approaches. 



A. The Nonlinear Sigma (NLcr) model approach 

The idea underlying the construction of this model is 
to consider the system in its low-temperature — sym- 
metry broken — phase and to take into account small 
fluctuations of the fields around the direction of the or- 
der parameter. The corresponding treatment is thus, by 
construction, a low-temperature expansion. Its actual va- 
lidity is in fact less stringent than that: it is sufficient that 
the system is locally ordered and that the temperature is 
small. This explains why this approach is valid even in 
two dimensions for systems obeying the Mermin- Wagner 
theorem. Note that this approach applies — a priori (see 
Section III Bl and the discussion at the end of this sec- 
tion) — only for N > 3. Indeed, in the N — 2 case, the 
low-temperature expansion of the NLcr model leads to a 
trivial result, i.e. the theory is perturbatively free. This 
result is however not reliable since there exist topological 
as well as Ising-like degrees of freedom in the XY frus- 
trated case (see Section [H}. These degrees of freedom, 
that are completely missed within the low-temperature 
perturbative approach, drastically affect the physics at fi- 
nite temperature as in the famous Berezinskii-Kosterlitz- 
Thouless phase transition jSJ E3| • 

Within the NLcr model approach, the partition func- 
tion of the 50(3) x 50(2)-symmetric model follows from 
the Hamiltonian Eq. Q together with the constraints of 
Eq. © Hi: 



z= fvfavfa Y[ 

i<3 



■fa ~ S ij) 



(40) 



The delta-functionals allow the integration of the three 
massive modes among the six degrees of freedom of fa 
and (f>2- Therefore, only the three — Goldstone - 
modes 7? re main, in term s of which the partition func- 
tion writes pMll60L IT6l | : 

Z = J^D7f CX P (-^ J d d X ff y(7T)fl7rWj . (41) 

The Eq. I|41|) . where gij(n) embodies the interaction, is 
the suitable expression for a low-temperature expansion 
of the SO (3) x SO(2)/SO(2) NLcr model. 



The low-temperature expansion of such NLcr mod- 
els has been studi ed i n general but rather abstract 
terms by Friedan |160| . The specific study of the 
SO (3) x 50(2) /SO{2) model and its generalization to TV- 
component spins — the O(N) x 0(2) /{0(N - 2) x Q(2)) 
model — has been performed by Azaria et al. |l5ll Il6lj 
(see also ^56]). The RG analysis requires to consider the 
most general Hamiltonian invariant under O(N) x 0(2) 
and renormalizable around d — 2. This Hamiltonian in- 
volves not only the usual kinetic terms for fa and fa, 
Eq. 1)41) [I. but also a nontrivial derivative term, called the 
"current-term", which reads: 



i d x (fa 



(42) 



This term must be included in the model since it has the 
right symmetry, is power-counting renormalizable around 
d = 2 and is thus generated during the RG flow. The 
correct NLcr model — in the sense of stability u nder RG 
transformations — is given by (for any N > 3) |l5l| : 




(Wi) 2 + (dfaj 



+ (f -f) ($i-dfa-fa.dfa)'' 



(43) 



where we have chosen to reparametrize the coupling con- 
stants in a way convenient for what follows. Now, the 
Hamiltonian of the naive continuum limit Eq. (|40|l is 
just the initial condition of the RG flow corresponding 
to rji = IJ2/2 = 1/T. Note that we have included the 
temperature in the coupling constants. 

For the special case N = 3, it is convenient to rewrite 
the model differently. We define, as in Eq. I|l(jf> . a third 
vector fa by: 



fa = fa A fa 



(44) 



With this expression, it is easy to verify that the current- 
term, Eq. 142|) , is nothing but a linear combination of the 
kinetic terms of fa, fa and fa: 



d d x($ 1 .dfa- 



d d *((dfa)\(dfa) 2 -(dfa)' 



(45) 



One can then gather the three vectors fa , fa and fa into 
a 3 x 3 matrix: 



(fa,fa,fa^j 



(46) 



Since (fa, fa, fa) are three orthonormal vectors, one has 
*$$ = ]! an d $ is therefore a 50(3) matrix. 
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The partition function thus reads: 



d d xTr(V <9*$<9$) 
Z = I V<5> <5('$$ - l)e J (47) 



where V is a diagonal matrix of coupling constants: V = 
diag(pi = p 2 = m/4:,P3 = Vi/ 2 ~ W 4 )- 

It is easy to check on Eq. I|47|l that the model is invari- 
ant under the right transformation: 



(48) 



with V being the subset of 50(3) matrices that com- 
mute with V. When p% ^ p%, i.e. r\\ rj2, V is 
isomorphic to 50(2). When V is proportional to the 
identity, V is isomorphic to the whole 50(3) group. In 
this last case, the high-temperature symmetry group is 
G = 50(3) x 50(3) - S0(4). Note that this identity 
has to be understood at the level of the Lie algebras since 
50(3) x 50(3) and 50(4) are locally isomorphic but 
differ globally and have different topological properties. 
This fact will be important in the following. 

The RG equations for the 0(N) x 0(2)/(0(JV - 2) x 
0(2)) mo del have b een computed at two-loop order in 
d — 2 + e |l5ll Il6lj . We recall here the one-loop result 
that will be useful in the following: 



/3 ni = -(d-2)m + N-2-£- 



p va = -(d-2)ri 3 



N-2 flfr 



(49) 



A fixed point is found for any N > 3. For N = 3, it corre- 
sponds to p* = p\, i.e. r/* — r\\ and, thus, to an enlarged 
symmetry 50(3) x 50(3)/50(3) ~ 50(4)/50(3). This 
fixed point has only one direction of instability — the 
direction of the temperature — and thus corresponds to 
a second order phase transition. Surprisingly, the critical 
behavior is thus predicted to be governed by the usual 
ferromagnetic Wilson-Fisher fixed point with the subtlety 
that it corresponds to ] our- component spins. Note that 
this precisely corresponds to the particular case consid- 
ered in Section llTni Another subtlety is that since, here, 
the order parameter is a matrix instead of a vector - 
it is a 50(4) tensor — the anomalous dimension is dif- 
ferent from the usual anomalous dimension of the four- 
component vector model. Only the exponent v is inde- 
pendent of the nature of the order parameter and is thus 
identical to the usual v alue of v of the Wilson-Fisher 
N = 4 universality class |l5lL Il6l| . 

In fact, it is easy to convince oneself that the fixed 
point found exists to all order of perturbation theory. Ac- 
tually, the crucial fact is that, in d = 2 + e, the perturba- 
tive (5 functions of a NLcr model associated with the sym- 
metry breaking scheme G — > H only depend on the local 
geometrical structure of the manifold G/H wh ich is itsel f 
determined by the Lie algebras of G and H |l5ll Il6l| . 
Since the Lie algebras of 50(3) x 50(3) and of 50(4) 



are identical, the perturbative f3 function for the — re- 
maining — coupling constant of the model with p\ = p% 
is identical at all orders to the perturbative function 
of the usual 50(4)/50(3) NLcr model. The existence of 
a fixed point for the 50(3) x 50(3)/50(3) NLcr model 
at all order of perturbation theory follows from the fact 
that its existence makes no doubt for the 50(4)/50(3) 
NLcr model. 

At the time of the first investigation of the O(N) x 
0(2)/(0(N - 2) x 0(2)) NLcr model, the most natural 
position was to extend this equivalence beyond perturba- 
tion theory and to assume that the 50(3) x 50(3)/50(3) 
fixed point exists everywhere between two and four di- 
mensions, as it is the case for the 50(4)/50(3) fixed 
point. This wa s, in parti cular, the position advocated by 
Azaria et al. |l5lL Il6lj . The outstanding fact is that 
although the 5Q(4) behavior has indeed been seen nu- 
merically in d = 2 0, |?3 , it actually does not exist far 
from two dimensions. This is clear since no such fixed 
point is found in d — 4 — e and since, as already empha- 
sized, the 50(4) behavior is not seen in any numerical or 
experimental data in d = 3. It is thus extremely proba- 
ble that either the fixed point disappears in a nontrivial 
dimension smaller than 3 or it survives in d = 3 while 
being no longer the usual N = 4 fixed point. Note that, 
in the first case, its 50(4) nature can also change before 
it disappears. Anyway, this fixed point must disappear 
below d — 4. The situation is thus more involved than 
in the "usual" 50(4)/50(3) model. There must exist 
nonperturbative reasons explaining the disappearance of 
the fixed point and/or the loss of its 50(4) character. 

Actually, it is clear that the perturbative low- 
temperature expansion performed on the NLcr model 
misses several nonanalytic terms in T — typically terms 
that behave as exp(— 1/T) — that could be responsible 
for the disappearance of the fixed point and/or its change 
of nature. There are, at least, two origins for such terms. 

1) The first one consists in the nontrivial topological 
configurations — see the discussion in Section lll Bl follow- 
ing Eq. I|15(l — that are completely neglected in the low- 
temperature expansion of the NLcr models. Indeed this 
expansion relies, by construction, on the local geomet- 
rical properties of the manifold G/H and is insensitive 
to its global — topological — structure. Thus it ignores 
vortex-like configurations that likely play an important 
role in three dimensions. 

2) The second origin of nonanalytic corrections to the 
low-temperature (3 function is more technical. The low- 
temperature expansion is performed in terms of the Gold- 
stone — or pseudo-Goldstone in d = 2 — modes that 
are represented by fields constrained to have a modulus 
less than one, see Eq. (|41[l . This inequality cann ot be 
taken into account in the perturbative treatment |162| 
and is thus relaxed, leading to neglect terms of order 
exp(— 1/T). All these terms are negligible for the criti- 
cal behavior when the critical temperature is very small, 
which is the case near d — 2. However, they become 
important when T c ~ 1 which is typically the case in 
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d = 3. 

Only a nonpcrturbative treatment can take into ac- 
count these nonanalytic terms and thus allows to follow, 
when the dimension is increased, the fate of the 0(4) 
fixed point. 



B. The Ginzburg-Landau- Wilson (GLW) model 
approach 

The GLW model for the O(N) xO(2)/(0(N-2) x 0(2)) 
model can be deduced from a generalization of Eq. (|40|l to 
iV-component vectors, by replacing the functional delta- 
constraint by the most general potential that favors the 
field configurations obeying the initial constraint. For 
convenience, we choose to parametrize it by: 



! J - Sij) — > e" 



(50) 



with: 

U = 



2 i 72 
+ <f>2 



X+jj, 
16 



($1 



(51) 



where, as usual, r is proportional to the reduced temper- 
ature while A and p are (/> 4 -like coupling constants. 

All field-dependent terms in Eq. 1)51(1 can be rewritten 
in terms of the rectangular matrix $ defined in Eq. JSJ. 
The corresponding Hamiltonian then reads: 



H = J d d x QlV (9^9$) + r - P + A p 2 + J T 



(52) 



with p = Tr(*$$) and r = ±Tr (*<M> - 11 p/2) being 
the only O(N) x 0(2) independent invariants that can 
be built out the fields, see Appendix IB1 Note that min- 
imizing the term in front of p corresponds to imposing 
*$$ cx 1, i.e. to imposing that 4>i and 02 are orthogonal 
and of the same norm in agreement with the characteris- 
tics of the ground state of frustrated magnets, see Fig.HJ). 



As well known, for any N > N c (d = 4-e) = 21.8+0(e) 
there exist four fixed points: the Gaussian — G — the 
vector 0(2N) — V — and two others, called the chiral 
— C+ — and anti-chiral — C_ — fixed points. Among 
these fixed points one, C+, is stable and governs the crit- 
ical properties of the system and the others are unstable 
(see Fig. I2a|l . When, at a given dimension d close to 
four, N is decreased, C+ and C_ move closer together, 
coalesce at N c (d) and then disappear (see Fig.l2b|l. More 
precisely, for N < N c (d), the roots of the (3 functions ac- 
quire an imaginary part. Since no stable fixed point ex- 
ists below N c (d) and since the flow drives the system in a 
region of instability, it is believed that the transition is of 
first order. Note that for N < N' c {d = 4-e) = 2.2 + 0(e), 
C+ and C- reappear but not in the physically relevant 
region to frustrated magnets. 

For completeness we give the exponent v at one-loop: 



1 ( (N- 3)(N + 4)V48 - 24A + N 2 

2 +£ { 8(144- 24A + 4N 2 + N 3 ) + 



(54) 



N(48 + N + N 2 ) 
+ 8(144- 24A + 4AT 2 +N 3 ) 

and recall that the anomalous dimension vanishes at this 
order. Note that the square root becomes complex for 
2.2 < N < 21.8, which is reminiscent of the critical values 
N c {d) and N' c (d) of the number of spin components, see 
above. 



2. The three and five-loop results in d = 4 — e 

In 4 — e dimensions, the cri tical value N c (d) has been 
computed at th ree-l oop order |l53j and, very recently, at 
five-loop order [l59j: 



N c (d-. 



21.80 



-23.43e 
-0.03e 3 



7.09e 2 
4.26e 4 



0(e 5 



(55) 

In fact, as it is often the case within this kind of ex- 
pansion, the series are not well behaved and it is dif- 
ficu lt to obta in reliable results even after resummation 
[22l Il53l Il59| . We h owever indicate the value found at 
three- loop order |l53j : N c (d = 3) = 3.39 and at five-loop 
order [l5g: N c (d = 3) = 5.45. 



1. The RG flow 

The RG equations for the coupling constants entering 
in Hamiltonian l|52|) have been co mpute d in the e = 4 — 
d-expansion up to five-loop order |159j and in a weak- 
coupling expansion in d = 3 up to six-loop order |l55j| . 
We recall here only the one-loop result of the e-expansion 
to discuss qualitatively the flow diagram: 



(3x 



-ep + 



1 



16tt 2 
1 

16tt 2 



(4Xp + 4p 2 + X 2 (N + 4)) 
(6Xp + Np 2 ) . 



(53) 



3. The improved three and five-loop results 

1 1 

It has been conjectured by Pelissetto et al. 156] that 
N c (d = 2) = 2 , a result which is however somewhat con- 
troversial [l59]. It is possible to use this nonperturbative 
information to reformulate the series obtained within the 
4 — e expansion. Imposing the constraint N c (d = 2) = 2 
to th e three-loop series, Pelissetto et al. have obtained 

-0.16e 2 ) + O(e 3 ) . 

(56) 



N c {d = 4-e) = 2+(2-e)(9.90-6.77e- 
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a) N > N c (d) b) N' C <N< N c (d) 



FIG. 2: Flow diagram for a) N above N c (d) and b) TV be- 
low N c (d). The fixed points C+ and C_ that exist above 
N c (d) coalesce at N c (d) and then disappear. G and V are the 
Gaussian and vector 0(2N) fixed points. 



in d = 3 (l5rj| : 

Around d ~ 4 and d ~ 2 the perturbative results of, re- 
spectively, the GLW and NLct models are recovered once 
the limit N — * oo has been performed. This suggests 
that, at least for sufficiently large N, the two models 
belong indeed to the same universality class in all di- 
mensions. However, wi thin this approach, no N c (d) line 
is found (see however 156]). It is thus impossible to 
extrapolate to finite N the results obtained in this ap- 
proach. 



Reformulated in this way, the coefficients of the series 
decrease rapidly. It is thus reasonable to use this expres- 
sion to es tima te N c (d — 3). Pelissetto et al. have thus 
obtained 156]: N c {d = 3) = 5.3(2) where the error bar 
indicates how N c (d — 3) varies from two to three loops. 
However, Calabrese and Parruccini have shown tha t, ex - 
tended to five loops, the same series behaves badly |l59j : 

N c {d = 4 - e) = 2+ (2 - e)(9.90 - 6.77e + 0.16e 2 + 
+ 0.06e 3 + 2.16e 4 )+O(e 5 ) . 

(57) 

Using different kinds of tricks, notably the inverse of the 
series, they have obtained, from the five-loop series, the 
value N c (d = 3) = 6.1(6). 

From this approach one is strongly tempted to con- 
clude that, in the physical cases N = 2 and N = 3, 
the transitions are of first order, even if it is impossible 
to conclude about the strong or weak character of this 
transition. 



4- The three-loop results in d — 3 

A weak-coupling analysis h as be en performed directly 
in d = 3 at three- loop order in |l52j . This leads to N c (d — 
3) = 3.91. However, as already emphasized, this result is 
not well converged. 



5. The large-N results 



6. The six-loop results in d = 3 

In three dimensions, a six-loop computation has been 
performed by P elissetto et al. |155j and re-examined by 
Calabrese et al. |l3ll Il32j , see below. The results arc the 
following: 

1) For N sufficiently large — N > 6.4(4) — there 
exist four fixed points, one stable and three unstable, in 
agreement with the usual picture given above, see Fig-EJa,. 
The transition is thus of second order. 

2) For 5.7(3) < N < 6.4(4), there is no nontrivial fixed 
point and the transition is expected to be of first order. 

3) For N < 5.7(3) and, in particular, for N — 2 and 
N = 3 a stable fixed point is found and a second order 
phase transition is expected. 

According to Pelissetto et al. |l55l Il56j| . the fixed 
points found for N = 2 and N = 3 should be non analyt- 
ically connected with those found in the 1/N and 4 — e 
approaches. Therefore, it should be impossible to obtain 
them by following smoothly those obtained at large- TV or 
close to d = 4. 

The critical exponents obtained by Pelissetto et al. are 
given in Table WTU Note that the error bars are about 
ten times l arge r here than in the ferromagnetic O(N) 
models [2^, Il62j ] computed by the same method. This is 
an indication that the resummed perturbative series are 
converging much slower than in the vectorial case. 

Let us now discuss these results. 



System 


Ref. 


a 


P 


7 


V 


XY 


[155] 


0.29(9) 


0.31(2) 


1.10(4) 


0.57(3) 


Heis. 


[155] 


0.35(9) 


0.30(2) 


1.06(5) 


0.55(3) 



TABLE VII: The six-loop perturbative results in d = 3. 



The large- N expansion was first performed by Bailin 
et al. 9]. It was th en r e-examined by Kawamura [7(j 
and Pelissetto et al. |l56fl . A fixed point is found within 
this expansion in all dimensions between 2 and 4. The 
exponents v and r\ have been computed up to order 1 /N 2 



The XY case. First, one should indicate that the expo- 
nents 7 and v computed from the six-loop approach com- 
pare reasonably well with the data of group 1. However, 
as already mentioned, the value of r\ found by the scaling 
relations must be positive when there exists a fixed point. 
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One finds, with the data of Table IVm T] ~ 0.08 which 
is significantly positive. Let us recall that this is not the 
case for the experiments performed on the materials of 
group 1 and for the numerical simulations performed on 
STA. Note, moreover, that the value of (3 found within 
the six-loop calculation, is very far — around four theo- 
retical error bars — from the average experimental ones 
which are = 0.228(6) for CsMnBr 3 alone, f3 = 0.237(4) 
for the whole group 1 and far from the numerical values 
obtained for STA (3 = 0.24 - 0.25. Thus, contrarily to 
what is asserted in |l55j , it seems extremely improbable 
that the exponents found at six loops could fit with those 
of group 1 and those of the numerical STA model. Actu- 
ally, this is also the case for the materials of group 2 for 
which the average (3 equals 0.389(7). 

The Heisenberg case. First, one notes that the agree- 
ment between the 7 and v exponents obtained from the 
six-loop approach and from the experimental or numeri- 
cal data is not as good as it is in the XY case. Concerning 
77, one finds, with the data of Table IVTT1 77 ~ 0.08. This 
has to be compared with the value of 77 obtained i) for 
the materials of group 1, which is significantly negative 

— 77 = —0.118(25) — ii) for materials of group 2, which is 
marginally negative — r\ = —0.018(33) — and Hi) in the 
simulations of the STA which is also negative although 
not completely significantly: r\ = —0.0182(89). The neg- 
ativity of 77 is an indication of a mismatch between the 
six-loop results and the data for the Heisenberg systems 
even if it cannot be used as a definitive argument against 
a second order phase transition. The exponent 7 ob- 
tained from the numerical simulations of the Heisenberg 
STA model provides a further information. The average 
value of this exponent — 7 = 1.185(3) — is rather far 

— 2 .5 theoretical error bars — from the six-loop results 

From the previous analysis one can conclude that, as 
such, the fixed point obtained within the six-loop ap- 
proach turns out to be not directly relevant to the phe- 
nomenology of XY materials or simulated systems. This 
seems to exclude the scenarios I, II and III that all as- 
sume that, at least, a certain number of compounds or 
systems are well described by a fixed point. 



7. Critical remarks 

As we mentioned at the beginning of our analysis of 
the experimental results, see Section fill Al we have made 
an assumption on the nature of the experimental errors 
which is not realistic: the systematic errors cannot be 
neglected. We now come back on this point and show 
that the conclusions we have drawn from our analysis 
persist without this assumption. 

Let us consider the XY case, where the symptoms of a 
mismatch between the theoretical and experimental re- 
sults are the clearest. We concentrate on the materials 
of group 1 and on the exponent (3 which is the best mea- 
sured, see Table [I] With our assumption, we have found 



(3 = 0.228(6). Let us suppose that, contrarily to our 
assumption, the systematic error is large and dominates 
the total error. Let us take: 

(3 = 0.23(2) (59) 

so that all experimental and numerical results lie in the 
interval of values, see Tables U and IIIII This estimate has 
to be compared with the six- loop result: 

(3 = 0.31(2) (60) 

where, in this case, the authors indicate that they have 
been very conservative in the estimation of the error bar 
jl55j . Although it is difficult to get fully unambiguous 
conclusions out of these numbers, it is clear that the 
agreement between them is not satisfactory. The same 
considerations on group 2 of XY materials lead to sup- 
pose: 

/3 = 0.39(2) (61) 

which, again, is far from being in agreement with the 
six- loop result Eq. (JSUJ). 

It is also possible to test the negativity of the anoma- 
lous dimension r\ with our new assumption. In the same 
spirit, one estimates v — 0.555(30). We find, in this case: 

-0.28 < 77 < -0.048 . (62) 

Thus, 77 is again found negative even in the most extreme 
hypothesis. 

We thus conclude that, although it is — up to now — 
impossible to estimate rigourously the confidence level 
of our analysis of the experimental data since only one 
error bar is given in the literature, it appears to be very 
difficult to reconcile the experimental and numerical data 
with the six-loop results. 



8. The six-loop results in d — 3 re-examined 

In order to cope with the discrepancy between the six- 
loop results obtained by Pelissetto et al. and the exper- 
imental and numerical data, Calabrese et al . have re - 
considered the resummed six- loop series |l3ll Il32l Il58| . 
They claim that they can account for the unusual prop- 
erties of the critical exponents for XY and Heisenberg 
frustrated systems in d = 3 — negative anomalous di- 
mension and weak universality — by the fact that the 
RG trajectories around the stable — focus — fixed point 
found by Pelissetto et al. are spiral-like. By integrating 
the resummed (3 functions for the two coupling constants 
of the GLW model and computing the effective exponents 
77 and v along the RG trajectories, they have found that 
these exponents display large variations in a transient 
regime. These authors argue that the scaling properties 
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of the system are governed, over several decades of tem- 
peratures, by the preasymptotic regime so that the effec- 
tive exponents observed experimentally can differ signif- 
icantly from their asymptotic values, i.e. those defined 
at the fixed point. 

Let us underline here several drawbacks of the scenario 
of Calabrese et al. . 

First, it is based on the existence of stable fixed points 
that are not related to any already known fixed point. In 
particular, the fixed points found for N = 2 and N = 3 
within this computation in d = 3 are, according to Pelis- 
setto et al. and Calabrese et al. non analytically related 
to those found in the large- N as well as in the 4 — e 
expansions. This means that there is no way to check 
their existence using these perturbative methods. This 
is specifically problematic in the context of frustrated 
magnets where the properties of the fixed points appear 
to be very unusual: i) the existence of the stable fixed 
points strongly depends on the order of perturbation — 
they are not present at three-loop order and only show 
up, as far as we know, at five- loop order ii) the location 
of the fixed points, as N and d are varied, seems to have, 
in the (N, d) plane, a very particular structure since, in 
three dimensions, they only exist when N is below an- 
other critical value of TV — which is found to be equal to 
5.7(3). 

Second, it is very difficult, in the computation of Cal- 
abrese et al. , to relate — even in principle — the initial 
conditions of the RG flow to the microscopic characteris- 
tics of real systems. This would require to handle the in- 
finity of coupling constants entering into the microscopic 
Hamiltonian obtained from the Hubbard-Stratonovitch 
transformation. This is impossible, at least within the 
usual perturbation theory. 

Third, it is very difficult to account, in this framework, 
for the first order behavior deduced from sever al numeri - 
cal simulations of XY and Heisenberg systems |l26t Il27j . 
We have also already noticed that XY-systems have a 
stronger tendency to undergo first order transitions than 
Heisenberg systems. However, there is no natural expla- 
nation for this phenomenon in the scenario of Calabrese 
et al. . 

Fourth, in this scenario, it is also very difficult to ex- 
plain why there is no physical system characterized by the 
asymptotic critical exponents, i.e. those corresponding 
to the fixed point found by Pelissetto et al.. This seems to 
require very unnatural experimental circumstances such 
that the initial conditions of the flow corresponding to 
the physical realizations of frustrated magnets are such 
that their long distance properties are never controlled 
by the nontrivial fixed point. 

Finally, there is no possible explanation of the break- 
down of the NLcr model predictions. 



C. Conclusion 

XY and Heisenberg frustrated systems exhibit the 
kind of problems we have described in the introduc- 
tion: the perturbative results obtained within a low- 
temperature expansion around two dimensions, within 
a weak-coupling expansion around four dimensions or 
within a large- TV expansion fail to describe their critical 
physics in three dimensions. Moreover, these different 
perturbative predictions are in contradiction with each 
other. Contrarily to the O(N) nonfrustrated case, there 
is no possible smooth interpolation of these results be- 
tween two and four dimensions and, at fixed dimension, 
between N = oo and N = 2, 3. More surprisingly and, 
again in contradiction with what happens in the O(N) 
nonfrustrated case, high-order calculations performed di- 
rectly in d — 3 also fail to reproduce the phenomenology, 
at least when they are interpreted in the usual way. This 
situation reveals the difficulties of the conventional ap- 
proaches to tackle with the physics of frustrated magnets. 
Only new interpretations or methods can allow to shed 
light on the problems encountered here. We have pre- 
sented the solution proposed by Calabrese et al. on the 
basis of a high-order pertubative calculation and under- 
lined its difficulties. We now present the nonperturbative 
method we have used to explain the unusual behavior of 
frustrated magnets. This is the subject of the next sec- 
tions. We start by a methodological introduction to this 
method and then apply it to the frustrated systems. 

VI. THE EFFECTIVE AVERAGE ACTION 
METHOD 

We now present the NPRG method we use : the effec- 
tive average action method [H IH 0, El IH^. The 
content of this section is neither original nor exhaustive. 
There exist seve ral well-documented reviews on the sub- 
ject |34l I45L IToH . Our aim here is only to provide some of 
the physical ideas underlying this method — notably the 
block spin concept and its formulation in the continuum 
— as well as its technical implementation on the simple 
example of the O(N) model. 

A. Block spin in the continuum 

The effective average action method, as well as many 
other NPRG techniques, i s bas ed on the well known con- 
cept of "block spin" [23, Il65j |: when dealing with any 
strongly correlated system, it is fruitful to integrate out 
the fluctuations step by step and, more precisely, scale 
by scale. In practice, one first gathers the initial — mi- 
croscopic — degrees of freedom into small "blocks" . It is 
then possible, at least formally, to integrate out, in the 
partition function of the system, the internal fluctuations 
of the blocks. This "decimation" is followed by a rescal- 
ing of length-scales, coupling constants and fields. In this 
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way, starting from a "bare" GLW Hamiltonian, one gets 
an effective Hamiltonian for the block degrees of freedom, 
i.e. for the low-energy modes. By iterating this proce- 
dure, one generates a sequence of — scale-dependent — 
effective Hamiltonians, parametrized by a running scale 
k, that all share the same long-distance physics. This 
sequence defines a RG flow. At a fixed point of this 
flow the system displays scale invariance. This allows to 
obtain the critical quantities through an analysis of the 
neighborhood of the fixed point in the flow of effective 
Hamiltonians pj. 

To illustrate how this concept of block spin is imple- 
mented concretely in the continuum, we consider the 
case of an Ising-like system, initially defined on a lat- 
tice which, in the continuum, is described by a scalar 
field C( x )- If the lattice spacing is given by a, the corre- 
sponding continuum field theory is characterized by an 
overall momentum cut-off A of order a -1 . The partition 
function writes: 



a one-dimensional integral: 



C] 



(63) 



where H 1 ^ [Q stands for the interacting part of the GLW 
Hamiltonian and: 



C-C7 X .C 



Wf C(q)C- (q )C(- 



(64) 



corresponds to the cut-off kinetic part. In Eq. (|63|l and 
(|64|l . Ca(c|) is an ultra-violet (UV) cut-off propagator 
that prevents the propagation of unphysical modes with 
momentum higher than A. One writes it: 



F k=A q 2 



(65) 



where Ffc(q 2 ) is a function of the ratio z = q 2 /k 2 that 
rapidly decreases when z — > oo. One also imposes to 
_Ffc(q 2 ) to be unity at the origin: F fe (q 2 = 0) = 1. A 
typical example of function F is: Fk (q 2 ) = e~( q / fc ) but 
other forms can obviously be considered. 

In Fourier space, the idea of block spin is specified by 
separating the low- and high-momentum modes of the 
spin-field £: 



C(q) = C>(q) + C<(q) 



(66) 



The fields C>(q) an d C<(q) being unconstrained, the sep- 
aration between high- and low-momentum modes is ac- 
tually realized through their respective propagator. We 
thus write: 

C A (q) = G fc (q) + (C A (q)-C fc (q)) = C < (q) + C > (q) (67) 

where k is the typical scale that separates the high- and 
low-momenta. In Eq. (|BTI) . C>(q) (resp. C<(q)) prop- 
agates C> (resp. £<), the high- (resp. low-) momentum 
degrees of freedom of the field £. This comes from a prop- 
erty of the Gaussian integral that can be easily seen on 



dz exp I — 



2{a + /3) 



+ f(z)\ cx 



J dxdy exp (— -|g +/(» + V)) 



(68) 



This result is easily obtained by changing, in the right 
hand side of Eq. H68JI . the integration variables x,y, into 
z = x + y and t = x — y and by integrating on t. 
Thus from Eqs. $57$ and JgHJ one gets: 

Z = Jv( K 2?C> exp (-lc<.C-\C< ~ l^-C; 1 -^ 

-ffHc< + c>]) • 

(69) 

The effective Hamiltonian H^l^] for the low- 
momentum degrees of freedom £< is defined through the 
integration over the high-momentum degrees of freedom 
in Eq. Pfr 



-H k nt [C<] 



2?C>exp(--C>.C; 1 .C> 

-#i nt [c< + c>: 



(70) 



Integrating out the internal degrees of freedom of a block 
spin between the scales a and a' > a corresponds, in this 
language, to the integration of the modes £> with mo- 
menta between k = a -1 and k' = a' . The Equation 
(|70|l implements the block spin procedure in the contin- 
uum which is the starting point of any NPRG approach. 



B. The Polchinski equation 

The effective Hamiltonian -ff™'[C<] follows an exact 
equation describing its infinitesimal evolution when the 
running scale k is lowered. To establish this equation we 
rewrite Eq. lf7UI) as: 

= Jvq exp - c<) • c; 1 ■ (C - C<) - #i nt [C 



= exp 



2<5C< 



5(< 



(71) 



This last functional relation can be inferred from the 
one-dimensional identity: 



7 d 2 



dx exp 



(72) 
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By differentiating each side of Eq. (|71|l with respect to k 
we obtain: 



-(d k H, 



rr int 



1 / Sffff* „ ^gj nt 



Finally, the exact evolut ion e quation fo r ijl n t , known as 
the Polchinski equation |166| (see also |l67j |h writes ex- 
plicitly: 



5 2 H int 



SH^ SH'l nt 



(74) 



^C<(qMC<(-q) <5C<(q)^C<(-q) 
or, in real space: 

W*[C<] - \ fd d x d d y duC^x - y). 



SH k nt 5H{ nt 



SC<(*)SUy) <5C<(x)^C<(y) 



(75) 



with dhC^x - y) = / (0r d k C > (q) e*>(*-»>. Note that 
in the preceding equations we have improperly used the 
same notation for the field C< and regulator C> and for 
their Fourier transforms. Note also that, in the following, 
we shall use the notation £ for . A graphical represen- 
tation of the Polchinski equation is given in Fig. [3] 



2 ^^^^ 2 



FIG. 3: A graphical representation of the Polchinski equation. 
The crosses represent the cut-off factor dkC > (q). The black 
circles with n-external legs correspond to the n-th functional 
derivative of H k nt with respect to the field. 

Let us first make some remarks about Eq. I|74|) . A first 
feature we would like to emphasize is that this equation 
involves the quantity H k nt which is the effective Hamil- 
tonian for the degrees of freedom that have not yet been 
integrated out, namely £<. The drawback with H k nt is 
that it is an abstract object that has no direct physical 
meaning since it is a function of a field C< that eventually 
fully disappears in the physical limit k — > 0, i.e. when all 
fluctuations have been integrated out. In particular, one 
should realize that £< is not a precursor of the order pa- 
rameter, i.e. not a local magnetization at scale k. Indeed 
this magnetization should come from a thermodynamical 
average at scale k while £< is just a stochastic variable 
that represents the low-momentum part of the original 
spin field and is thus, roughly speaking, a spatial average 



of this field. Consequently, the effective Hamiltonians by 
themselves do not contain all the information on the in- 
tegration of the high-momentum degrees of freedom . 
For instance, the computation of correlation functions for 
the high-energy field £> would require to first couple the 
system to a source J — a magnetic field — by adding in 
Z a term exp( J.£) and to follow the flow of this term in 
order to obtain the full J-dependence of Z, a rather dif- 
ficult task. Thus, Eq. lf71|) provides at best a flow of the 
running coupling constants that parametrize the effective 
Hamiltonian H k nt at scale k. 

As shown mainly by Wilson, equations like Eq. Ij74(l 
are, in principle, sufficient to compute the critical ex- 
ponents once a fixed point Hamiltonian H k nt has been 
found. Actually, even for the evaluation of the RG flow, 
Eq. (|74|l suffers from an important difficulty: although 
this equation looks simple — its only nonlinearity is a 
term quadratic in H k — it is nevertheless a functional- 
partial-integro-differential equation that has no known 
solution in general. Therefore, in order to render it man- 
ageable, one has to truncate the Hamiltonian H™ 1 . 



1. Derivative expansion 

A natural truncation consists in an expansion of the 
effective Ham ilton ian in powers of the derivatives of the 
field 0, 0, Il68j . For instance, for a one-component 
scalar field theory one has: 



d d x 



u k (0 + lz k (0(dc) 2 + o(d 4 ) 



where Uk(() stands for the potential - - i.e. the 
derivative-independent part — of the effective Hamilto- 
nian and Zk(C) is the quadratic — field- dependent — field 
renormalization. With such a truncation, one neglects 
higher-order derivative terms. This is justified i) when 
one is interested in the long-distance, low-energy physics, 
since these higher-order derivative terms should corre- 
spond to less important operators and ii) when there is 
no qualitative change of nature between the microscopic 
and macroscopic degrees of freedoms — such as the ap- 
pearance of bound st ates at a finite scale k — that could 
induce non localities |169| . A practical guide to evaluate 
the validity of the derivative expansion is the value of the 
anomalous dimension rj. If this quantity is small, one can 
expect that the inclusion of higher-order derivative terms 
provides small corrections to the results. 

At first order in the derivative expansion one sets 
Zk{() = in Eq. I|76() and derives an RG equation for 
C/fc(C) from Eq. (|74|l . This corresponds to the so-called 
Local Potential Approximatio n (LPA) which has been in- 
tensively explored in the past jH, ll70L Il7ll Il72| . In par- 
ticular, this kind of approach has been used by Zumbach 
to analyze th e physi cs of frustrated magnets in three di- 
mensions 13 Ell • The problem with the LPA is that, 
since by definition it neglects the field renormalization, it 
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leads to a trivial — vanishing — anomalous dimension. 
Consequently: i) this prohibits to compare the results 
obtained within this approach to that of a standard per- 
turbative approach when this last one involves a nontriv- 
ial anomalous dimension (this is, for instance, the case 
of the NLcr model around d — 2 already at one- loop), ii) 
this prevents an accurate evaluation of critical exponents 
for systems for which the anomalous dimension is not ex- 
pected to be small. In the context of frustrated magnets, 
these drawbacks are serious since we are precisely inter- 
ested in relating the different pcrturbative approaches 
and, to some extent, by a satisfactory determination of 
the critical exponents. We thus need to compute the field 
renormalization. 

This kind of computation however encounters several 
difficulties. First, whereas the RG equation for U k (() in 
the LPA of the Polchinski equation is universal — cut- 
off independent — , the RG equation derived for U k (() 
and Zk(C) a t second order in the derivative expansion 
explicitly depends on the regulator C > (q) chosen to sep- 
arate the high- an d low-energy degrees of freedom in 
Eq. (EZJ) |l68t Il73f . Another related problem is that 
of reparametrization invariance. The partition function 
(|firtf) and, thus, the physical quantites like critical expo- 
nents, are invariant under a general change of field of 
the kind £ — > £ + G(C) where G is an arbitrary function 
starting at order £ 2 . As a consequence of this invari- 
ance, the normalization of the field Z k (£ — 0) in the 
Hamiltonian is a priori an arbitrary parameter. Unfor- 
tunately, the reparametrization invariance is broken as 
soon as one performs a truncation of the Hamiltonian. 
As a result the critical exponents and, in particular r], 
depend on the normalization Z k {Q — 0). It follows from 
these considerations that, in any practical computation, 
one encounters the problem that physical quantities de- 
pend on nonuniversal parameters such as cut-off func- 
tions and normalizations. Different techniques, such as 
the Principle of Minimum Sensitivity (PMS), have been 
used to decrease the de pendence o f the critical quantities 
on the cut-off function |l73l 1174) - Also, some criterions 
have been proposed to find the best normalization, i.e. to 
find a value Zk{C = 0) such t hat t he derivative expansion 
converges the most rapidly |174| . These considerations, 
having for aim to exploit the Polchinski equation at the 
next to leading order in derivative expansion, have led 
to the determination of rather satisfactory critical expo- 
nents. 

At the same time, there has been a great activity de- 
voted to the search of other formulations of the RG ideas 
that could avoid some of the troubles encountered in the 
use of the Polchinski equation. The effective average ac- 
tion method is the result of this search. 



C. The effective average action method 



mental object, not the abstract effective Hamiltonian 
£/l nt [C]> functional of the stochastic low-energy field 
but, rather, the Gibbs free energy T — called effective 
action in field theory — functional of the order parame- 
ter (j) = (£}. To implement this idea in the RG context, it 
is necessary to build a running Gibbs free energy Fk for 
the high-energy modes that have already been integrated 
out at this scale k. The argument of Tk is, therefore, the 
order parameter at this scale that eventually becomes, 
when i->0, the true order parameter. 

These requirements imply several constraints on the 
definition of T^. First, at the scale of the lattice spacing, 
k = A = a , T k should correspond to the microscopi- 
cal Hamiltonian H since no fluctuations have been taken 
into account. Second, when the running scale k is low- 
ered to 0, Tfc, which then includes all fluctuations, must 
identify with the standard effective action T from which 
all thermodynamical quantities like magnetization, cor- 
relation length, etc, are computed. To summarize, T k 
must respect the constraints: 

r fc=A = h 

(77) 

r fc=0 = r 

and has to interpolate smoothly between these two limits. 

1. Construction 

Let us again consider, for simplicity, the case of a sys- 
tem described by a scalar field C( x )- The construction of 
the effective average action proceeds in two steps. First, 
one should decouple the low-energy modes — with mo- 
menta q 2 > k 2 — in the partition function in order to get 
a theory involving only the high-energy ones that will be 
summed over. Second, in this modified theory, one builds 
the Gibbs free energy, as usual, by a Legendre transform. 
This gives T^. Let us now study how this is implemented 
in practice. 

The first step is conveniently implemented by changing 
the partition function Z into Z k for which a fc-dependent 
term, quadratic in the fields and thus analogous to a 
mass-term is added to the microscopic Hamiltonian 0, 
Eif . With this "mass-term", the partition function in 
presence of a source J writes: 

Z k [J] = J VC, cxp ( - H[Q - AH k [(] + J.c) (78) 
with J.C = J d d q /(q)C(-q) and 

AHklC] = \j 10^~ ^•( t l' t l')C(q)C(q / ) (79) 



The basic — and physically very appealing — idea with 7\L fe (q,q') = (27r) d <5(q + q')R k (q 2 ). In Eq. (jHUJl, 
of this new formulation is to consider as the funda- R k (q 2 ) is the cut-off function that controls the separation 
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between the low- and high-energy modes. To decouple 
the low-energy modes, it must act as a large-mass term 
for small q whereas it must vanish for large q to keep 
unchanged the high-energy sector of the theory. Thus: 



and 



Rk(q 2 ) ~ k 2 for q 2 < k 2 



i?fc(q 2 ) -> when q 2 > k 2 



(81) 



(82) 



The first constraint means that, for momenta lower than 
k, R k {q 2 ) essentially acts as a mass — i.e. an IR cut- 
off — which prevents the propagation of the low-energy 
modes. The second ensures that the high-energy modes 
are fully taken into account in Z k [J] and thus in the effec- 
tive average action. Moreover, since we want to recover 
the original theory when k — > 0, i.e. when all fluctua- 
tions have been integrated out, R k (q 2 ) must vanish in 
this limit. Thus we require: 



i? fc (q 2 



identically when k — > 



(83) 



which ensures that Zk=o[J] = Z[J}. On the other hand, 
when k — > A, i.e. when no fluctuation has been in- 
tegrated out, Tk should coincide with the microscopic 
Hamiltonian. This is achieved by requiring (see below 
for the proof): 



i? fc (q 2 



identically when k — > A 



(84) 



Note that, since we shall not be interested in the pre- 
cise relation between the microscopic characteristics — 
defined at scale A — of a given system and its critical or 
pseudo-critical properties, we set A = oo in the following. 
A widely used cut-off function is provided by : 



i?fe(q 2 ) 



Zq 2 



e q 2 /fc 2 



1 



(85) 



where Z is the field renormalization. Including it in R k 
allows to suppress the explicit Z dependence in the final 
RG equations — see below. The cut-off function i?^(q 2 ) 
corresponding to Eq. (|85|l is plotted on Fig. 0] Another 




FIG. 4: A typical realization of the separation of high- 
and low-momentum modes provided by the cut-off function 
i?fe(q 2 ). At low momentum, i?fe(q 2 ) acts as an effective mass 
of order Zk 2 while the high-momentum behavior is not mod- 
ified. 

useful cut-off fu nctio n, called theta cut-off, has been pro- 
posed by Litim |l75j . It writes: 



R k (q 2 ) = Z (k 2 - q 2 ) Q(k 2 -q 2 



(86) 



where is the usual step function. 

The second step consists in defining the effective aver- 
age action. The free energy at scale k is given — up to 
a factor —ksT — by: 



W k [J] = \nZ k [J] 



(87) 



From Eq. I|87[l. one defines the order parameter 4> k (q) at 
scale k as the average value of the microscopic field £(q) 
in the modified theory: 



Mq) = (C(q)> = 



SW k [J] 
SJ(-q) 



Thanks to the properties of i?fc(q 2 ), the contribution to 
the average value in Eq. (|88() coming from modes with 
momenta q 2 <C k 2 is strongly suppressed. Also 4>k(q) 
identifies with the true order parameter in the limit k — > 
0. Note that, for simplicity, we omit, in the following, 
the index k to <f> k . 

The effective average action is defined by |43j : 



T k [ct>} = -W k [J} + J4-^H k [<t>] 



(89) 



where J = J [(/>], sec Eq. Thus T k [cj)] essentially cor- 

responds to a Legendre transform of W k [J] for the macro- 
scopic field <f) — up to the mass-like term AH k . The rela- 
tion i|89|) implies several unconventional relations. First, 
taking its derivative with respect to 4>(q) provides the 
relation between the source and T k : 



J(-q) = 



Taking the derivative of this relation with respect to c 
implies a second important relation |21l| : 



where if > (q,q') = 5 2 r fc /<ty(q)<ty(q'). 

Let us now show that the definition of Tk, Eq. (|S3|l . en- 
sures that it satisfies the requirements given in Eq. l(77|l , 
i.e. that it interpolates between the microscopic Hamil- 
tonian H for k = oo and the (true) effective action T for 
k — > 0. This last property follows directly from Eq. H89|) 
and the fact that for k = the IR cut-off i?^(q 2 ) identi- 
cally vanishes. The fact that T k identifies with H when 
k — > oo can be shown in the following way. One deduces 
from Eqs. (J7HJ), l|57)l. l|85) and J^Ol the functional iden- 
tity: 



K k {q,q') 5J(-q) 



(2tt) m 
= (2*)- 



Kq') 

5 2 W k 



(90) 

M) 

(91) 
(92) 



e -T k [(f)} = Jv( exp(-if[C]- 
5T k [4>] 



+ ^.(C-0)-Alfc[C-*]) 



(93) 
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In the limit k — > oo, R k (q 2 ) goes to infinity. In this limit, 
the mass-term exp(— AH k [( — 4>]) acts as a hard con- 
straint on the functional integration — exp(— AH k [( — 
<f)}) ~ 6[( - <j>] — so that r fc=oo [0] = H[<j)]. With these 
properties, Y k [<j>\ has the meaning of a coarse-grained 
Gibbs free energy at scale fc _1 : lowering k corresponds 
to including more and more fluctuations. 



2. The equation 

Let us now derive the exact RG equation for IV We 
start from the expression: 



oW h [J] 



J PC exp(-ff[C]-Aif fe [C] + J.( 



(94) 



which results from Eqs. 1781) and 1)8 7|) . One first writes 
the variation of exp(Wfe[J]) with respect to the scale k 

d k e W k[J] = 

= - J V((d k AH k [C]) exp ( - H[{] - AH k [(] + J.c) 



(d k AH k 
d d q 



_6_ 

SJ 



from which follows 
1 
2 



d k W k [J] 



{2ix) d \5J(q) 
d d q 



,W k [J] 



d k R k (q 2 



SJ(-q) 



d k R k (q 2 ) . 



(2TT) d 

6W k [J]6W k [J] S 2 W k [J] 
6J(q) SJ(-q) + SJ(q)5J(-q) 



e W k [J] 
(95) 



(96) 



which looks like the Polchinski equation (|74|) . 

Let us now differentiate the expression l|89ll with re- 
spect to k, at fixed (f>: 



d k T k [cj>} 



(S.J 



+ (d k J).<t>-d k AH k [<t>] (97) 

= -d k W k [J] - d k AH k [4>] 

Together with Eq. i|96[) this gives: 

A Q o ^ S 2 W k [J] 



(2) 



'5J(q)8J(-q) ' (98) 
Using Eq. (|92|) one obtains an equation involving only 
and its second functional derivative T 



.Tr (27r) 2d r^ 



U k (q,-q) (99) 



with t = lnfc and i? fe = <9 t i?fc- In Eq. (JHU, Tr must be 
understood as a trace on internal indices — vectorial or 
tensorial — if C spans a nontrivial representation of a 
group. 

Let us finally give a form of Eq. Q99JI more convenient 
for practical use |212j : 



d t r k [4>] = 



2(2tt 



^a t Tr(ln(( 27 r) 2 ^[ 



(100) 

where the "time derivative" dt only acts on R k , i.e. 
dt=R k d/dR k and where the trace Tr now also means a 
momentum-integral / d d qd d q' (2n)- d 8(q + q'). Eq. l|9"9"|) 
(or Eq. (11001) ) controls the evolution of T k with the run- 
ning scale k. According to the preceding discussion, it 
describes, when k is lowered, how the running effective 
average action is modified when more and more (low- 
energy) fluctuations are integrated out. 



3. Properties 

We now give some important properties of Eq. I|99|l . 
The reader interested in more details can consult Ref. 
@. 

1) Eq. 1|99|) is exact. It thus contains all perturbative 
|!76lll77| and nonperturbative features of the underlying 
theory: weak- or strong-coupling behaviors, tunneling 
between different minima |l78j| , bound states |l69l Il79| , 
topological excitations [55|, etc. 

2) While it has been derived here in the case of a one- 
component scalar field theory, Eq. I|99|l obviously holds 
for any number of components and, more generally, for 
any kind of order par ameter. The extension to fermions 
is also trivial (see [45j for instance). 

3) With a cut-off function R k (q 2 ) which meets the 
condition l|81() or, more generally, with a finite limit 
when q 2 — > 0, the integral in Eq. (|99|l is infrared (IR) 
finite for any k > 0. This IR finiteness is ensured 
by the presence of the mass-term R k which makes the 
quantity T y k ' [<p] + R k positive for k > even at the 
critical temperature. This allows to explore the low- 
temperature phase even in presence of massless — Gold- 
stone — modes. From the UV side, the finiteness of the 
integral in Eq. I|99|) is ensured by a requirement of fast 
decaying behavior of R k (q 2 ). 

4) One can give a graphical representation of Eq. ij9"9"|) . 
see Fig. |SJ It displays a one-loop structure. Obviously 
this one-loop structure must not be mistaken for that of 
a weak-coupling expansion. Actually, the loop involves 
here the full — i. e. field-dependent — inverse propaga- 



tor rjj. 2 ' [<j)] so that the graphical representation of Fig. |S] 
implicitly contains all powers of the coupling constants 
entering in the model. Note also that this one-loop struc- 
ture automatically ensures that all integrals over inter- 
nal momenta involved in this formalism have a one-loop 
structure and are thus one-dimensional. Thus they can 
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be easily evaluated numerically and, when some partic- 
ular cut-off are used, analytically. This radically differs 
from a weak-coupling expansion which leads to multiple 
loop diagrams and thus, multiple integrals. Another im- 
portant feature of Eq. (|99|l is that very simple ansatze on 
Tk allow to recover in a unique framework the one-loop 
perturbative results obtained by standard perturbative 
calculations around two and four dimensions as well as 
in a large- N expansion. 




FIG. 5: A graphical representation of Eq. I|99^ . The cross 
represents the function Rk and the line the propagator 

(rf (V + R*)- 1 - 

Let us make a final remark. The one-loop structure 
of Eq. H99|) contrasts with the Polchinski equation l|74|) 
which, in addition to the loop term, involves a tree part 
(see Fig. |3J). An interesting consequence of the struc- 
ture of this "Legendre version" of the NPRG equation 
is that reparametrization invariance is preserved by the 
derivative expansion when a power-law cut-off is used. 
This means that with such a cut-off function, the anoma- 
lous dimension is no longer ambiguously defined |l80j| . 
This would apparently select the power-law cut-off as 
the best one. The situation is, in fact, more involved 
since the power-law cut-off is afflicted with bad conver- 
gence properties when used within the derivative expan- 
sion, ft has appeared that for instance, the exponential 
cut-off — Eq. dHSJ — or the theta cut-off — Eq. 
— that do not respect the reparametrization invariance 
of the RG equation, lead to better results when opti- 
mization critera are used. We do not enter into more 
details in th ese proble ms of reparametrization invari- 
ance Ejj. |4?1 1174. Il80l and optimiz ation of the results 
pi l5aTl75 ni8lL 11821 Il8.il ITsl Il8fi | and refer to the lit- 
erature. The main reason for this is that, as we shall see 
in the following, we shall only deal with pseudo-critical 
exponents that, being given their lack of universality, i.e. 
their strong dependence with respect to the microscopic 
parameters, makes superfluous an optimization of the 
computations. 



4- Truncations 

As it is the case for the Polchinski equation, Eq. I|99|) 
is too complicated to be solved exactly. Its nonlinearities 

are even worse than in the Polchinski case since it involves 

(2) 

all powers of T k . As a consequence, the functional 
has to be truncated. Different kinds of expansions have 
been considered 0: 

1) Field expansion where Tk is expanded in powers of 
the order parameter cf>. For a scalar field theory, one has: 



r ^] = E^ / Ii ddx < <Xxi)...<Mxn) ri n) (xi,..,x„) 

n=0 ' •* i=0 

(101) 

where rj^(xi, ...,x„) denotes the n-th functional deriva- 
tive of Tk ■ 

2) Derivative expansion where Tk is expanded in pow- 
ers of the derivatives of the order parameter. For a scalar 
field theory, one has: 

T k [<t>] = J d^fukW + ^Zki&idcf)) 2 +0(d i )\ (102) 

Uk{(f>) being the potential — i.e. derivative- independent 
— part of Tk while Zk(4>) corresponds to the kinetic term. 

3) Combined derivative and field expansions where one 
further expands in Eq. (|1L)2|I the functions Uk{4>) and 
Zk(4>) in powers of <f> around a given field configuration 
(fro- Technically, this kind of approximation allows to 
transform the functional equation l|99|) into a set of ordi- 
nary coupled differential equations for the coefficients of 
the expansion. In practice, it is interesting to consider an 
expansion around (one of) the field configuration (j) that 
minimizes the potential Uk- For the simplest — Ising — 
model, this expansion writes: 

T k [4>] = Jd d x( 

\u' k '{po){p - Po) 2 + ^U'k'(p Q )(p - p Q f + • • • (103) 
+ iz fe ( / 5 )(6V>) 2 + \z' k {po){p ~ Po) ipcpf + ■■■) 

where p = ^4> 2 an d Po — ^<t>o 2 , 4>o being the magnetiza- 
tion at scale k. The rational behind this choice is that the 
minimum of U k is physically the location that we want 
to describe the best since thcrmodynamical quantities 
at vanishing external field are determined from the mini- 
mum of Tk at k = 0. The relevance of such a parametriza- 
tion is confirmed by many works showing that the con- 
vergence of the critical quantities, when more and more 
powers of the field <f> are added in the truncation, is im- 
proved when compared with the same calculation per- 
formed with an expansi on of Uk (4>) and Z k {^>) around 
the 4> = configuration p M ll8 ^ . 

The choice of a good truncation is a complex problem. 
One has to choose a manageable truncation that however 
encodes the relevant physics. In practice, it appears that, 
surprisingly, even at low-orders in powers of derivatives 
and fields, Eq. I|99|l provides correct qualitative features 
of the RG flow. However, the precise determination of 
the critical quantities requires to push the expansion to 
rather large orders in the field and involves a heavy alge- 
bra mug. 

To illustrate how the technique works we now consider 
the simplest case, i.e. the vectorial O(N) model. The 
O(N) x 0(2) model relevant to frustrated magnets is 
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technically more involved but the procedure to derive 
the RG equations follows the same steps. Details of the 
technicalities in this latter case will be given in Section 

rvrri 



D. The O(N) model 

We present here the e ffecti ve average action approach 
of the 0(N) model Hl[l88|. We essentially follow the 
presentation given, for instance, in 44] with some dif- 
ferences, ensuring a self-contained presentation. We use 
a truncation similar to the one we use to deal with frus- 
trated magnets, i.e. where Tk is expanded both in deriva- 
tives and fields. Let us first consider the derivative expan- 
sion of the effective average action for the O(N) model 
at order d 2 : 



= J d d x (u k ( P ) + x -z k ( P ) (d$f+ 



+ -Y k (p) {dpf + Oid') 



(104) 



where is a iV-component vector field and p = (f> 2 /2. 
In Eq. (|104fl . U k (p) is the potential — i.e. derivative- 
independent — part of Y k while Z k (p) and Y k (p) cor- 
respond to kinetic terms. These two last terms embody 
the renormalization for the Goldstone and massive fields, 
respectively. Note that the term proportional to (dp) 2 is 
always absent from the GLW action used for a pertur- 
bative analysis in coupling constant. The reason for this 
is that, in this context, it is power-counting irrelevant. 
On the contrary, in the context of the effective average 
action method, there is no perturbative expansion and, 
thus, no power-counting argument works. One, however, 
expects that the terms of lowest degrees in the field (for 
d > 2) and in the derivative are the most important for 
the long distance physics. 

The case Z k (p) = 1 and Y k (p) = in Eq. (|104|l corre- 
sponds to the LPA. A nontrivial anomalous dimension 
is obtained by going beyond this simplest truncation. 
As said above, we use here a truncation that mixes the 
derivative and field expansions. We thus consider: 



\u' k \p Q ){p - p Q f + ^U k (p )(p - a,) 3 + 

-\z k {p Q ) {d$f + \z' k (p ){ P - ~ Po )(d$) 2 + ■ 

~Y k (p ) (dp) 2 + \Y k \~ Po )(p - po) (dp) 2 + ■ 



(105) 



where po = \<\>\\ parametrizes the ^-dependent field con- 
figuration that minimizes U k . Since our aim here is only 
pedagogical and not devoted to the calculation of pre- 
cise critical quantities, we consider the following ansatz 



which is limited to the smallest expression providing an 
nonvanishing anomalous dimension: 



r fc f 



d d x l -Z(d$) 2 



1 



u 2 {p- po) 



(106) 



where Z=Z k (po) and U2=U k (p~a). This approximation 
looks very much like the GLW Hamiltonian used to 
study perturbatively the O(N) model, up to a trivial 
reparametrization. There is, however, a major difference. 
Here the ansatz (|106fl is not studied perturbatively in the 
(j) coupling constant U2- It is to be inserted in the exact 
RG equation (|9^|) . 

Let us now establish the RG equations for the coupling 
constants entering in Eq. i|10tj|) ■ The calculation proceeds 
in four steps: 

i) We first define the running coupling constants po, 
it2 and Z from functional derivatives of the ansatz of 
T k , Eq. (|106f) . This is analogous to imposing renormal- 
ization prescriptions for the renormalized coupling con- 
stants in usual perturbative calculations. As in this case, 
the coupling constants are defined as (combinations of) 
functional derivatives of — the "vertex functions" - 
taken in a specific field-configuration of the model. How- 
ever, contrarily to the perturbative approach which is 
generally performed in the high-temperature phase and 
thus, around a zero field configuration, we perform this 
expansion around a nontrivial running field configuration 

00- 

ii) We then apply the operator <9 t on these definitions. 
This is implemented by the use of the evolution equation 
(|99|l or l|100|) . The flow equations for the coupling con- 
stants are then expressed as traces of products of vertex 
functions that are evaluated from the ansatz Eq. (|106fl . 

Hi) The flow equations involve integrals over the inter- 
nal momentum. It is convenient to express these integrals 
in terms of dimensionless functions, known as threshold 
functions . The properties of these functions are such 
that they make explicit the phenomenon of decoupling of 
massive modes, see below. 

iv) Also, as usual, one introduces dimensionless renor- 
malized quantities to study the scale invariant solutions 
of the RG equations. 



1. Definition of the coupling constants 

Let us first choose one of the uniform field configura- 
tions that minimize the effective average action T k : 



<j) Mln (x) 



MA 
o 



V / 



(107) 



or equivalently: 



iMin 



(q) = (2^)% SnS(q) 



(108) 
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where 4>q — (2po) 1 ^ 2 is a /c-dependent quantity. Due to 
the O(N) symmetry of r fc , which is preserved at any t by 
the RG flow Eq. , the choice of a particular direction 
for cf> Mln is irrelevant and thus does not affect the RG 
equations. 

Let us now define the coupling constants. To do this 
we introduce the notation: 



5 n T k [ct>\ 



k{ ai ,n},...,{a n , Vm } 5</> ai (sn)...5<t> an (p„) 



(109) 



As said above, po specifies the position of the — run- 
ning — minimum of F^. It is implicitly defined by: 



-(i) 

k {a,0} 







(110) 



Mil 



where the notation "Min" refers to the configuration 
given in Eq. H107JI . Because of our particular choice 
of ^ the previous equality is trivially satisfied for 
a = 2, . . . , N and we shall consider only the case a = 1 
in the following. 

The coupling constant ui is defined along the same line 
as: 



"2 



2p 5(0) 



.(2) 
k {1,0}, {1,0} 



(111) 



Mil 



Finally, the k— dependent field renormalization Z is ob- 
tained by considering the term quadratic in momentum 
of a momentum dependent configuration: 



5(0) 



dp 2 



,(2) 

k {2,p},{2,-p} 



Min 



(112) 



In this last equation, the index 2 specifies a direc- 
tion orthogonal to that defined by the minimum (see 
Eq. (|107[) 1. Note that one could have considered any 
of the N — 1 directions orthogonal to that defined by the 
minimum. Note finally that the 5(0) term appearing in 
Eqs. (|lll|l and (|112|l is the volume of the system and is 
present here since T k is an extensive quantity while the 
coupling constants are not. 



2. The t-denvation 

The flow equations for the coupling constants po, u 2 
and Z are obtained by derivating, with respect to t, the 
previous definitions Eqs. piOfl . (|lll|l and (|112ll . 

Let us start by po- One has to take care of the t- 
dependence of both rjj, 1 ^ j and its argument, the con- 
figuration (f) Mln — Eq. HlU7f) — which has a nontrivial 
t-dependence through that of 0o : 



d, F 



.(i) 

k {a,0} 



Min 



- F> r (1) 

— U * L k {a,0} 



Mil 



,(2) 

k {q,0},{1,0} 



Mil 



(2n) d dtpo (113) 



The RG flow for po follows from this equation, taken 
for a = 1: 



dtPo 



y^> 1 k {i,o}, {1,0} 

In the same way, one obtains: 



°t L k {1,0} 



(114) 



Min 



_ (2n) d (2) 
dtU2 ~ 2pV(0) dtr Mi,0},{i,0} 



Mil 



(2n) d dtpo 
2p 6(0) 

(^) d r (3) 



-r 



Po 



(2) 

ft {1,0}, {1,0} 



(115) 



Min 



Jp^-k {1,0}, {1,0}, {1,0} 



Min 



and: 
d t Z 



<J(0) p^odp 2 



n F (2) 

°t L k {2,p},{2,-p} 



Mil 



(27T) d d t p ( 3 ) 

L k {2,p},{2,-p},{l,0} 



I2fv 



Mi; 



(116) 



The RG flows for the coupling constants po, "2 and 
Z involve successive functional derivatives of dtTk with 
respect to different </>;(qj). These quantities are easily 
obtained from Eq. (|100|) . Let us take its derivative with 
respect to ^(qi). Using: 



ii (qi) 



In 



((27T) 



2d r (2) 
1 ft 



n k ) = 



\2<Z r (3) 



v / ft {u,qi} 
where we have introduced the notation: 



V r , 



v r = ((2^) 2d rf +K k 



(117) 



(118) 



one obtains: 



(i) 



f - i r i 

1 ft {n,qi} 



(2-) rf 5 



^ Tr { r i 3 k,qi}-^} ( 119 ) 



for the one-point vertex function. Note that, in the right 
hand side of the preceding expression, we have only spec- 
ified the external indices {ii,qi} and omitted the inte- 
grated (or summed over) variables. The dot is here to 
remind that these integrations and summations have to 
be performed. The equation (|119H admits a graphical 
representation: 




(120) 



In this representation, the external leg implicitly car- 
ries an index of internal symmetry ii and a momentum 
qi. Now taking the derivative of Eq. H119|) with respect 
to 4>i 2 (t\2) and using: 



-V r .T 



(3) 



fc {*,q} • ' r, 



Vr, 



(121) 
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one obtains: 



3. The renormalization group flow 



(2) 



* k {ii,qi},{«2,q2} 



(2tt 



■(4) 



1 k {ii,qi} ' ' r - 1 k {i 2 ,qi} ■ 1 r j 

which can be graphically represented by: 



(122) 



o F (2) 

U t L k {ii,qi},{j2,q 2 } 



(123) 



We now explicitly write the flow equations 
for the coupling constants. This requires to 
know the vertex funct: 



: ons taken at the minimum 
appearing in Eqs. (|114fl . 

' L J Mm 

( 11151 and (|116fl . To evaluate them, one uses the 
truncation Eq. l|lUtj|) . One obtains: 



k {Ql,Pl},{ a 2,P2}, -,K,Pn} 



r 



(i) 

k {ii.qi} 



Mi. 



r 



(2) 

k {n,qi},{i2,q 2 } 







= [Zqi di li2 + 2p Q u 2 d lll d l2l ) — — ^ — 

Min (27T) a 



.(3) 

k {ii,qi},{i2,q 2 }d i 3,q3} 



Mi: 



2po - "2(^ii 2 ^3l + ^2«3 1 + ^ii^l) 



<5(qi + q2 + q3) 



(124) 



r 



(4) 

{*i.qi},{»2.q2},{«3,q3},{u,q4} 



Min 



^2(^2122^3*4 ~^ ^2123^2244 ^2124^2223) 



(2tt) m 
^(qi + q2 + q3 + q4) 



(2tt 



is of 



Mil 



(2) 

In this last set of equations, T) ', ■ „ x <■ „ x 

^ ' k {2l,qi},{22.q2} 

particular interest since its inverse provides the propaga- 
tor V r at scale k and, thus, the spectrum of excitations 
of the theory, at this scale. We easily get from Eq. (|124fl : 



^!"{ii:qi},{2'2,q2} 
1 



Min 



(2n) d S(q 1 +q 2 ) 



Zqf + R k (qf) 
1 



Zqi 2 + i? fc (qi 2 ) + 2/5ou 2 







if i\ = i 2 = 1 



if ii = 12 7^ 1 



if it ^ i 2 



(125) 



where i?fc(qi 2 ) is the contribution of the regulating term 

ijHDl- 

It is clear on the expression (|125JI that the N x N 
matrix 7 > r{ii,qi},{j 2) q2} ^ s diagonal. This holds inde- 

Min 

pendently of the kind of truncation used. The spectrum 
of excitations around the minimum, at scale fc, is thus di- 
rectly red on Eq. (|125|) . We find — up to the Rk term — 
one massive mode of squared mass 2poU2 in the longitudi- 
nal direction and TV — 1 massless modes in the directions 
orthogonal to the magnetization <f> . The deformations 
of the vector cj> associated to these modes are represented 
in Fig. E| 



a) massive excitation 



b) N - 1 "Goldstone" 
excitations 



FIG. 6: Schematic description of the deformations of the 
vector <j> associated with the proper excitations of the O(N) 
model: dotted arrows display the configuration chosen at the 
minimum of I\ and plain arrows display the relevant defor- 
mations: a) massive singlet, b) "Goldstone" N — 1-uplet. 



It is important to keep in mind that this spectrum 
corresponds to effective masses — at scale k — for which 
only high-momentum fluctuations — higher than k - 
have been considered. It is only in the limit k — > that 
one retrieves the physical spectrum. In particular, we 
stress that a qualitative change in the spectrum can oc- 
cur when k is varied. For instance, the following situation 
can happen: for large k, the minimum of Uk is nonvan- 
ishing so that, at this scale, the system behaves as if it 
was in its broken phase. However, when k is decreased, 
the minimum moves toward zero and eventually vanishes 
for some k > 0. Thus, while the system, for k = A, looks 
as if it was in its broken phase, it is actually, i.e. for 
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k = 0, in the high-temperature phase. This is what hap- 
pens when the temperature lies between the true critical 
and the mean-field critical temperature. Another sub- 
tlety is that, in order to analyze the critical behavior, we 
have to consider the dimcnsionlcss renormalized quanti- 
ties. Again, one has to be careful about the conclusions 
deduced from the behavior of these quantities. For in- 
stance, the dimensionless renormalized position of the 
minimum at k = can be nonvanishing even at the crit- 
ical temperature whereas the "true" magnetization is, of 
course, vanishing at T c . This is possible since the dimcn- 
sionful quantities are the products of their dimensionless 
counterparts and of positive powers of k. 

Using Eqs. (EH, ljTT5|) . (jnni, Hill), (EH and fiSty . 
one obtains the flow for pq: 



dtpo 



for u 2 : 



N-l 



-8 f ^ 

2 '7 (2ir) d \Zq 2 + R k (q 2 ^ 



Zq 2 + R k (q 2 ) + 2p u 2 J ' 



(126) 



a - ui~ [ d d q 

o t u 2 = — —o t 



N - 1 



2 J (2tt)<* \z(q 2 +R k (q 2 )) 2 

9 

+ (Zq 2 + R k (q 2 )+2p u 2 ) 2 

and for the field renormalization Z: 



(127) 



d t Z = -2p u 2 lim —— fi 



d d q 



1 



P^odp 2 ^ J (2ir) d Zq 2 + R k (q 2 ) 
1 \ 



Z(p + q) 2 + R k ((p + q) 2 ) + 2p u 2 ) ' 



(128) 



The search for fixed point requires to introduce dimen- 
sionless renormalized coupling constants. We define: 



Po = Zk z - d p 
u 2 = Z- 2 k d - 4 u 2 



(129) 



These changes of variables are deviced so that k and Z 
disappear from the flow equations of the renormalized 
dimensionless quantities. 

The corresponding flow equations thus write: 

dtPo = -(d~2 + 7])p + 

+ 2v d (N- l)lf(0) + 6v d l d (2u 2Po ) 
d t u 2 = (d - 4 + 2ij)u 2 + 

+ 2v d (N - \)u\ l d (0) + I8v d u\ li{2u 2 p ) 

(130) 

that depend on Z only through rj, the running "anoma- 
lous dimension" : 



V 



-fit log Z 



(131) 



In our truncation, it is given by: 



V = 



16v d 



ulpo mi 2 (2u 2 p ) ■ 



(132) 



The usual anomalous dimension is given by the fixed 
point value of Eq. (|T32j) . In Eqs. (|T3U|) and fEEty . we 

have introduced the — dimensionless — threshold func- 
tions it and TO, 



'2,2- 



ryn i, — d +2n 

itM = ~ ft 



4u d 



d d q 

J2nf 



1 



2 H 



(Zq 2 + R k (q 2 ) 

dZ 2 k 6 - d 



8v d 
1 



+ Zk 2 w) n 
lim 

P 2 ^o dp 2 



d d q 

(2n) d 



Zq 2 + R k (q 2 ) 



1 



' Z(p + q) 2 + R k ((p + q) 2 ) + Zk 2 w J 

(133) 

with v^ 1 = 2 d+1 TT d / 2 T(d/2). Some properties of these 
threshold functions are provided in the Appendix [0 We 
concentrate here on the main physical aspects of the 
threshold functions: 

1) Note first that the arguments of the threshold func- 
tions entering in Eqs. I|130l) and (|132fl are either or 
2u 2 po that are — up to R k (q 2 ) — the dimensionless 
renormalized square masses associated with the excita- 
tions around the minimum. 

2) The threshold functions l d (w) and m d 2 (w) decrease 
as power-laws when their arguments increase: 



l d {w)(xw 
m 22 (w) tx w 



-n-i 

2 



(134) 



for w ^> 1. The RG equation l|99|l makes thus explicit 
the phenomenon of decoupling of massive modes: if the 
renormalized square mass M 2 — here 2u 2 po — of a mas- 
sive mode increases when the scale k is lowered, the con- 
tribution of this mode to the flow becomes negligible be- 
low a scale k c defined by M kc ~ 1. 

3) The threshold functions are nonpolynomial func- 
tions of their arguments. Thus the flow equations l|130|) 
and 1)13211 are nonpcrturbative with respect to the cou- 
pling constant u 2 as well as to 1/po which, as we shall see, 
is proportional to the coupling constant — the tempera- 
ture T — that parametrizes the perturbative expansion 
of NLcr model. 

As we now show, the effective average action approach 
allows to recover the perturbative results obtained at low- 
temperature around d = 2, at weak-coupling around d = 
4 and in a X/N expansion. 
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4- The weak- coupling expansion around d — 4 

Just below four dimensions, the nontrivial fixed point 
governing the phase transition of the O(N) model is very 
close to the Gaussian fixed point. This justifies to expand 
the RG flow equations and (|132|) both in the cou- 

pling constant u 2 and in e = 4 — d. At lowest order, 
the function r\ is vanishing. Since po remains finite, the 
quantity 2u 2 po is of order e and one can perform a small 
mass expansion. The flows of the coupling u 2 and of po 
are obtained using ^(w) ~ ^(0) — nuil^ +1 (0) for id <C 1 
and 4(0) = 1. This leads to: 



dtPo 
d t u 2 



N- 



(N- 



16tt 2 



lf(0) 



8tt 2 



"2P0 



(135) 



-eu 2 



16tt 2 



At leading order, the roots of these equations are the 
gaussian fixed point - u\ — and p$ = (N + 
2)1 * (0)/327r 2 — and a nontrivial fixed point obtained for 
u* = 167r 2 e/(iV + 8) and p* = (N + 2)^(0)/32tt 2 . One 
easily deduces the critical exponent v from Eqs. 
linearized around the nontrivial fixed point: 



1 e 7V + 2 

2 + 4 7V + 8 



(136) 



It coincides with the one-loop expression obtained within 
a perturbative weak-coupling approach of the corre- 
sponding GLW model in d = 4 — e. 



5. The low-temperature expansion around d — 2 

Let us now show that Eq. (|99l) also allows to recover 
the one-loop results obtained around two dimensions in 
a low-temperature expansion of the NLcr model. We first 
make contact between the parameters — essentially the 
temperature — of the O(N) NLcr model and those of the 
effective average action (|106|) . 

The partition function of the O(N) NLcr model is given 
by: 



Z = Jv$s($ 2 -\) exp (~J d d * (defy 



(137) 



Let us replace the delta-constraint by a soft constraint. 
Using the field redefinition <p — > 4>VT one obtains: 



T>(j) exp 



\jd d *((d$? 



\(cf> 2 T-l)' 



where the delta-constraint is recovered when A — > oo. 
Comparing this expression with the ansatz 1106(1 and us- 
ing the relation (|129f) one obtains the relation: 



T 



1 

Wo 



(139) 



As a consequence, the low-temperature one-loop pertur- 
bative results can be recovered from a 1/po expansion. 
In fact, since the dimensionless renormalized mass of the 
massive modes is given by 2u2Pa, one has to perform a 
large-mass expansion. Physically, this corresponds to the 
known fact that, around two dimensions, the longitudinal 
modes of the O(N) linear model are frozen and only the 
— transverse — "Goldstone" fluctuations are activated. 
This phenomenon corresponds to the decoupling of mas- 
sive modes. Technically, this is realized through the be- 
havior of the threshold functions. As already stated, the 
threshold functions decrease as power-laws for large argu- 
ments, Eq. (|134fl . As a consequence, in the flow equations 
(|130fl . the contribution of the massive mode — which 
is proportional to if^poX) — is subdominant compared 
with the contribution of the Goldstone modes — pro- 
portional to ^(0). Now, using the large-mass expansion 



,d=2 



+ 0(uj~ 3 ), one gets from Eq. lfT52)l : 



47TP0 



(140) 



Using this expression of the anomalous dimension and 
the value lf(0) — 1, one gets: 



dt Po ' 
d t u 2 



N 



-epa 



4ir 

-2u 2 + ^—!-u 2 2 l 2 2 (0) 



(141) 



The flow equation for po coincides exactly with the re- 
sult obtained in the one-loop low-temperature expansion 
of the NLcr model for the temperature — which is given 
by Eq. (|139fl . The fixed point coordinates are given by: 
p* = (N- 2)/(47re) and = 8tt/((AT - l)if(0)). This 
leads to the critical exponents: 



n 



N 
1 



(142) 



v = — 
e 



which identify with those given by the low-temperature 
pertu rbative expansion of the NLcr model at one-loop 

Note that the perturbative [3 function for it 2 - 
Eq. (|135f) — and for p — Eq. l|141f) — are universal, 
i.e. independent of the cut-off function JZfe(q). Indeed, 
these j3 functions only depend on the values of the thresh- 
old functions ^ n (w) at u> — and on mf^ 2 ^) at large 
u> that, as shown in Appendix IC 31 do not depend on the 
cut-off function Rj.(<\). The matching with the pertur- 
bative results obtained around d = 2 and d = 4 is a very 
important feature of the effective average action method. 
First, it allows to interpolate smoothly between two and 
four dimensions in a unified framework. Second, it sug- 
gests that it is possible to reliably explore the behavior 
of the system in any dimension d and, in particular in 
d = 3, see below. 
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6. The large-N analysis 

The flow equations (|130[) and (|132[) can easily be ex- 
panded in the large- N limit. The leading contributions 
come from the Goldstonc modes, which appear with a 
multiplicative factor N. The (3 functions then read: 

d tPo = (2 - d) Po + 2Nv d i?(0) 

(143) 

d t u 2 = (d-4)u 2 + 2Nv d ul l*(0) 

where we have anticipated that the anomalous dimension 
vanishes at leading order, see below. The fixed point solu- 
tions are easily found to be: Pq = 2N v^lf{Q)/{d—2) and 
u\ = [d — 4)/(2vdNl%(0)). From these results, we check 
that the anomalous dimension behaves like 1/N and thus 
gives subdominant corrections to the (3 functions. We 
can, finally, compute the critical exponents by diagonal- 
izing the stability matrix. We then find v = l/(d — 2), in 
agreement with the standard leading order result of the 
1/N expansion. 

We now check that the effective average action method 
provides reliable results in three dimensions. 



7. The critical exponents in three dimensions 
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TABLE VIII: The critical exponents in three dimensions for 
the O(N) model, a) corresponds to the truncation where 
only the flow of {Z,po,U2} is considered. In b) one adds 113. 
In c), one adds {113, U4, Y{po)} and in d) {113, 114, Yo, Z'(po)}. 
e) takes into account the full dependence of Uk, Zk and Yk 
in the field |Hsj| . In /), the order d 4 terms of the derivative 
expansion have been included jHflj . g) correspon ds to the five- 
loop resummed perturbative results in in 4 — e Il89||. h ) are 
seven-loop perturbative results in three dimensions 189]. 



One of the main interest of the effective average ac- 
tion method is its ability to tackle with the physics in a 
nonperturbative regime, precisely when there is no small 
parameter, as it is the case in three dimensions. Wc 
provide, in Table IVIIII the critical exponents v and 77 ob- 
tained with this method, as functions of the order of the 
derivative and field expansions of Uk(p), Zk(p) and Yk(p) 
(see Eq. QlUSfl l. We have also included the results of high- 
order perturbative expansion for comparison. The expo- 
nent v is rather poorly determined with our simple trun- 
cation (|ll)6fl . However, the precision improves rapidly 
when more terms are added to the ansatz for Y k - For 
the best truncation, v is determined at less than one per- 
cent compared with the world best estimates. Although 
we shall not be concerned in the following in truncations 
using the full potential Uu(p) and the full kinetic terms 
Zk(p) and Yk(p) entering in Eq. i|lU4f> . we have indicated, 
in Table IVlTTl the critical exponents computed with such 
ansatze. One notes that v is in ver y good agreement 
with seven-loop resummed series |l89j | while the anoma- 
lous dimension is less precisely determined until the order 
<9 4 terms of the derivative expansion have been included 
in the ansatz, see [50J. 



8. The XY and multicritical Ising models in two 
dimensions 

Let us close this section devoted to the analysis of the 
O(N) model by a discussion of the results obtained in d = 
2 for the XY and Ising models. These are, in fact, two of 



its most spectacular successes because they correspond 
to truly nonperturbative systems. 

As well known, the physics of the XY model at fi- 
nite temperature is governed by topologically nontrivial 
configurations — vortices — which are not taken into 
account in a low-temperature treatment. According to 
Eq. (|141fl . the flow for po — or T — vanishes identically 
in d = 2 and A^ = 2 so that the theory is free. However, 
as well known, the model actually exhibits a phase tran- 
sition at a finite temperature Tbkt — the Berezinskii- 
Kosterlitz-Thouless phase transition — induced by the 
deconfinement of the vortices, see Remarkably, 
the simplest RG equations (|130|l and (|132|l already al- 
low to obtain the correct qualitative behavior of the XY 
model at finite temperature: a very small /3-function of 
T is found between T = and a finite Tbkt [54j. Re- 
cently, treating the full field-dependence of Uk, Z^ and 
Yk , von Gersdorff and Wetterich [55| have recovered the 
correct behavior for the correlation length of the XY 
model around Tbkt: 

e ^»(w^tf) ■ (144 > 

The exact results are r = 1 12 an d r\ = 1/4 for the anoma- 
lous dimension at Tbkt j53lll90j ] while von Gersdorff and 
Wetterich have found r = 0.502 and 77 = 0.287. This 
shows that the physics of topological excitations is cap- 
tured by the lowest orders of the derivative expansion, 
without including explicitly these degrees of freedom a la 
Villain [191) . 
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As for the Ising model, it is known that, in two dimen- 
sions, it can undergo infinitely many nontrivial kinds of 
phase transitions assoc iated with infinitely many mul- 
ticritical fixed points 192]. It is shown in Ref.0] 
that they all correspond to strong coupling fixed points. 
They are therefore very difficult to study by perturba- 
tive means. By a systematic search of fixed points in the 
two-dimensional scalar field theory, using an order d 2 
truncation of the derivative expansion, Morris |6l| has 
been able to find explicitly the first ten fixed points of 
this series. He has also shown that no other fixed point 
exists but the multicritical fixed points. 

9. A difficulty related to the field expansion 

Let us finally mention a difficulty linked to the field 
expansion of the potential Uk showing up, for instance, 
when the stable fixed point of the O(N) model is followed 
from d — 4 down to d = 2. When the simplest 4 trun- 
cation, Eq. H1U6|) . is used no problem occurs: one can 
smoothly follow the stable — critical — fixed point from 
d = 4, where it identifies with that found in a weak cou- 
pling expansion of the GLW model, down to d = 2, where 
it coincides with that obtained within a low-temperature 
expansion of the NLcr model, Eq. i|141|) . However, once 
the <fi 6 term is added, a new nontrivial fixed point emerges 
from the Gaussian fixed point in d = 3. This is a tricrit- 
ical fixed point, i.e. a fixed point with two directions of 
instabilities. As d is lowered, the critical and tricritical 
fixed points move closer together and eventually coalesce 
in a dimension 2 < d < 3. Actually, they both become 
complex. Note that when d is further lowered, the two 
fixed points become again real. In d = 2 + e, the stable 
fixed point can be identified with that found within the 
low-temperature expansion of the NLu model with the 
</> 4 truncation. Thus, there exists a small region between 
d = 2 and d = 3 where one fails to correctly describe the 
fate of the stable fixed point of the model using the (f> 6 
truncation. However, this is just an artefact of the field 
expansion, not of the method. To show this, let us de- 
scribe what happens when the order p of the truncation 
is increased. First, when including a new monomial <j) p 
in the effective potential, a new — multicritical — fixed 
point emerges from the Gaussian fixed point in the di- 
mension 2p/(p—2). Again the stable fixed point coalesces 
with one of these multicritical points and reappears close 
to d = 2. Second, as p increases, the coalescence of the 
stable fixed point occurs at smaller and smaller dimen- 
sions. Thus, one recovers a better and better description 
when increasing the order of the truncation. Also, it 
has been checked that when the full field-dependence of 
the potential is kept, the problem fully disappears and 
the stable fixed point can be followed smoothly between 
d — 4 and d — 2 |l93fl . Finally, it is important to indicate 
that, in the whole range of dimensions where the stable 
fixed point exists within a field expansion, the critical ex- 
ponents found within this approach at sufficiently large 
order p {p > 10) and those found within a full poten- 



tial computation are very close. The artefact of the field 
expansion described here can be bypassed using either a 
full potential computation or using a field expansion at 
sufficiently high order. Actually, it is not surprising that 
difficulties occur with the field expansion at low dimen- 
sions since the engineering dimension of the field vanishes 
as d — > 2. This strongly suggests that no power of the 
field can be safely discarded in Uk when d — > 2. This is 
confirmed by the fact that the effective potential, which 
is exactly known at N = oo for d — 3 and d = 2, is re- 
spectively a polynomial of order six and an infinite series. 



E. Conclusion 

We have described, in this section, the main features of 
the effective average action method. We now summarize 
them: 

1) The effective average action method allows trivially 
to recover the perturbative results around the upper — 
d = 4 — and lower — d = 2 — critical dimensions and 
thus to make contact with these approaches. 

2) The results obtained via this method are nonper- 
turbative in the different parameters: coupling constant 
and temperature. In this sense, it provides an alterna- 
tive approach to the usual perturbative methods. This 
is of great interest, especially when one suspects that the 
perturbative series could be not reliable as it is the case 
for frustrated magnets. 

3) Even with a very simple truncation of the effective 
average action, it is possible to capture some genuine non- 
perturbative features — like nontrivial topological con- 
figurations — that are unreachable from a conventional 
low-tempcrature expansion. This aspect is particularly 
important in the context of frustrated magnets since one 
knows that the low-temperature expansion performed in 
d = 2 + e does not provide the correct physics in d = 3, 
a possible explanation being the presence of vortex-like 
configurations in these systems. 



VII. THE O(N) x 0(2) MODEL 

We now come back to the study of frustrated magnets. 
We derive the flow equations relevant to the study of 
frustrated magnets. The derivation follows the same lines 
as in the O(N) case (see Section IYl Dl above') . 



A. Truncation procedure 

As emphasized previously, since the NPRG equation 
(I1UUI) cannot be solved exactly, a truncation for Tk is 
needed. We consider here a truncation involving only 
terms having at most two derivatives. At this order, the 
most general form of the O(N) x 0(2) effective average 
action writes: 
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= J d d x \U k (p,T) + \z k {p,r) ((d$i) 2 + (d$ 2 ) 2 ^j + \y1 1] { P ,t) (V x • 50 2 - 2 ■ d$ x y 
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We recall that 0i and 02 are the two TV-component vec- 
tors that constitute the order parameter, Eq. JSJ) while 



P 



= Tr('$.$) and r = |Tr (*$.$- 1 p/2) 



with 



$ = (0i,02) — are the two independent O(N) x 0(2) 
invariants (see Appendix [5] for a more detailed discus- 
sion). The truncation (|145fl is the analogue of Eq. (|104f) . 
in the case of matrix fields. Here Uk{p,r) is the poten- 
tial part of the effective average action while Zk(p, r) and 
Y^\p,r), i = 1,2,3, are kinetic functions. 

At this level of approximation, the RG analysis re- 

(i) 

quires to specify the five functions [/&, Zk and Y fc , 

1 = 1,2,3. This is to say an infinite number of cou- 
pling constants. As in the case of the vectorial 0(N) 
model, we further simplify the ansatz by expanding these 
functions in powers of the fields. Again, we choose to 
expand around a nonvanishing field configuration which 
minimizes Tk- This constraint is satisfied when 0i and 

02 are orthogonal, with the same norm. We choose: 
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the physical results being independent of this particular 
choice. The quantity *J~k entering in Eq. (|146[1 is analo- 
gous to the quantity 0o in the O(N) case, see Eq. <|107[) 
and we refer to it in the following as the magnetization. 

While studying the critical properties of the system, 
we have considered various truncations differing by the 
number of monomials in p and r included in the field ex- 
pansion. Our largest truncations consist either in keeping 
all terms in Uk up to the eighth power of the fields and 

(i) 

all terms in Zk and in the K; s including four powers of 
the fields or all terms in Uk up to the tenth power of the 
fields and the first term of the expansions of Zk and of 
Yu . With these truncations, we have verified that our 
results are stable with respect to addition of higher pow- 
ers of the fields. However, in order to keep our presenta- 
tion concise, we have chosen here to consider a reduced 
truncation that already enables to recover the different 
perturbative results — in 4 — e, 2 + e and 1/N — in their 
respective domains of validity. Within this truncation, 
we expand Uk up to terms containing four powers of the 

(i) 



fields and keep only the leading terms of Zk and Y h 
We also completely discard the two other functions Y£ 
and Y,. . This choice is motivated by the fact that, as we 



(2) 



will see in the next section, only the function Y^ con- 
tributes directly to the physics of the Goldstone modes 
and is thus important around two dimensions. Since one 
of our aims is to recover the results obtained around two 
dimensions, we keep this term in our ansatz. We are then 
led to the simple truncation: 
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Let us now discuss the different terms appearing in this 
expression. The coupling constants A and p have been 
already introduced in the GLW approach (see Eq. J5TJ). 
The coupling constant k describes the position of the 
minimum of the potential and appears in the truncation 
because we expand around the nonvanishing field con- 
figuration <J> Mm . As in the O(N) case, Z corresponds to 
the field renormalization. Finally, the unusual kinetic 
term with coupling u) corresponds to the current-term of 
Eq. H42fl introduced in the discussion of the NL<r model 
approach. This term is irrelevant by power counting 
around four dimensions since it is quartic in the fields 
and quadratic in derivatives. However its presence is 
necessary around two dimensions to recover the results 
of the low-temperature approach of the NLcr model since 
it contributes to the field renormalization of the Gold- 
stone modes. Being not constrained by the usual power 
counting one includes this term in the ansatz. 

The above effective action has all the ingredients to de- 
scribe accurately the physics at low-temperature around 
two dimensions as well as at weak-coupling regime 
around four dimensions. We can therefore anticipate that 
this simple truncation is actually rich enough to recover 
the perturbative results around d = 2 and d — 4. Of 
course, since our main goal is to go beyond the usual 
perturbative expansion, we have studied larger trunca- 
tions and have controlled the convergence of our results. 



1. The spectrum 

We now discuss the spectrum of excitations around the 
minimum l|146|) . The spectrum is given by the eigenval- 
ues and eigenvectors of the matrix S 2 Tk/S4> J i S4> l k — where 
i, k £ {1,2} and j, I S {1, . . . , N} — evaluated in the con- 
figuration (|146fl . We find that the 2N degrees of freedom 
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of the order parameter $ divide in four types that are 
described in Fig. [7J 

1) a family of 2N — 4 massless — Goldstone — modes 
which correspond to rotating rigidly the dihedral (0i, 02) 
by keeping either 0i or 02 unchanged, see Fig. I7al 

The four remaining modes correspond to the situations 
where the two vectors (0i, 02) remain in the same plane: 

2) a massless — Goldstone — singlet mode correspond- 
ing to rotating the dihedral within its plane, see Fig. I7bl 
Together with the 27V— 4 other ones, this gives the 27V— 3 
Goldstone mode of the model. 

3) a massive singlet of square mass Xk corresponding 
to a dilation of the two vectors, see Fig. [7c] 

4) a massive doublet of square mass jlk corresponding 
to fluctuations of each vector of the dihedral, with the 
constraint that the sum of the lengths of the vectors 1 0i | + 
1 02 1 remains unchanged, see Fig. I7dl 



model. They are given by: 
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a) Massless 
(2N-4)-uplet 




02 



>:>, 



b) Massless singlet 



c) Massive singlet 




d) Massive doublet 



FIG. 7: Schematic description of the deformations of (0i, <f>2) 
associated with the four types of proper excitations of the 
model. The dotted arrows display the ground state configu- 
ration and the plain arrows display the relevant deformations: 
a) massless (27V — 4)-uplet, b) massless singlet, c) massive sin- 
glet, d) massive doublet. 

In practice, it is very useful to work in the basis of the 
proper excitations of the model since then, the propa- 
gator being diagonal, the calculations are greatly simpli- 
fied. We therefore introduce 27V directions in the internal 
space, corresponding to the 27V proper excitations of the 
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In this basis, the two-point vertex function, i.e. the 
inverse propagator — up to the Rk term — has the form: 
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In the matrix (|149[1 . the first three lines correspond to 
the massive modes and the last 2N — 3 to the massless 
modes. Note that a nonstandard kinetic term appears 
on the fourth line through an additional field renormal- 
ization ujk for the Goldstone singlet. Let us add that if 
one keeps, in the truncation (|147J) . contributions from the 
functions and (see Eq. (|145|l h the field renor- 
malizations in the first three lines get extra contributions 
similar to what is obtained in the fourth line. Note that 
y( 2 ) and Y^ 3 - 1 affect only the field renormalization of mas- 
sive modes. It is therefore not necessary to take them into 
account in order to retrieve the leading order behavior in 
a low-temperature expansion around d = 2 which is en- 
tirely governed by Goldstone modes. This is why we do 
not keep them in our simplest truncation i|147|) . 
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B. The flow equations 

We now display the flow equations for the coupling 
constants entering in the truncation Eq. (|147|) . We recall 
the four major steps of this procedure (see Section lYl D) : 

i) The running coupling constants are defined as func- 
tional derivatives of the ansatz of T^, Eq. (|147|l . 

ii) The operator dt is then applied on these definitions. 
By making use of the NPRG equation 1(100(1 , flow equa- 
tions for the coupling constants are obtained as traces 
of vertex functions. These expressions are evaluated by 
using the truncated form of I\- Eq. 1(147(1 . 

Hi ) The flow equations are expressed in terms of thresh- 
old functions. 

iv) Dimensionless renormalized quantities are intro- 
duced. 



1. Definition of the coupling constants 

As in the vectorial model, the coupling constants are 
defined as values of the vertex functions in the specific 
configuration $ Mm around which is made the field ex- 
pansion, Eq. 1(146(1 . 

Let us start with the definition of k. This coupling con- 
stant parametrizes the ground state configuration $ Mln . 
One has, as in the O(N) case, an implicit definition of R: 







(150) 



with 5 a ,p given by Eq. 1(148(1 . In the following, as in the 
O(N) case, we shall consider only the case a = 1. 

The other coupling constants are defined using the two- 
point vertex function in different directions: 
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spectrum discussed previously (see Eq. 1(149(1 1. 

We finally define the coupling constants associated 
with the momentum-dependent part of our truncation 
1(147(1 . i.e. the field renormalization factor Z and the 
current-term coupling constant u>: 
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2. The t-derivation and the flow equations 

We now apply the operator d t to the definitions ((1501 - 
11521) . In order to derive the flow equations, we have to 
compute the functional derivatives of dtTk with respect 
to the fields. This is similar to what has been done pre- 
viously in the context of the O(N) model (see Section 
IVID 21 above), except that the tensorial structure in the 
internal space is more involved so that the computation 
of the traces is more cumbersome. We do not give the 
details here. We now introduce the dimensionless renor- 
malized quantities defined as: 



A = Z- 2 k d - A \ 

n = z- 2 k d -^fi 

uj = Z-' 2 k 
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as well as the threshold functions which are defined and 
discussed i n Ap pendix IU1 We then get the following flow 
equations |l!8j : 



These two definitions come directly from the study of the 
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= (d - 4 + 2r?)/x- Ivdii 



dui 
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VIII. TESTS OF THE METHOD AND FIRST 
RESULTS 

This section is devoted to all possible tests of our 
method in the O(N) x 0(2) case. We show, in particular, 
how the various perturbative results are recovered as it 
was already the case in the O(N) model. We also give 
our determination of N c (d) which is compared with the 
three-loop improved perturbative computation. Finally, 
we give our determination of the exponents in the N =6 
case and we compare them with those of the Monte Carlo 
simulation. 



One can also expand the (3 function for k, Eq. (|154afl : 

^_ (2 _ e)s+ M( w+1 + f )_ 

* \ / 9 \ (156) 

3kA k/j, i 2/i\ 
~ T^r 2 ~ 8^ \ + T J 

from which we can deduce the expression of v at order e, 
which coincide with the one- loop result of Eq. it5~4"|) . 

B. The low-temperature expansion around d — 2 



A. The weak-coupling expansion around d = 4 

Around d = 4, wc expect a nontrivial fixed point close 
to the gaussian. One can expand the flow equations at 
leading order in the quartic coupling constants and in 
e, as we did in the O(N) case (see Section IV I D 4|) . As 
expected from power counting, we find that the fixed 
point value of the coupling constant uj associated with 
the current-term is vanishing at leading order. This is 
also the case of rj. As in the O(N) case, the square masses 
Ak and \xk are of order e so that the threshold functions 
can be expanded in powers of their arguments. Once 
this expansion is performed one recovers the standard 
one-loop /3-functions for the coupling constants A and /i 
given in Eq. I|53|) that we recall here: 
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As explained in the context of the 0(N) model (see 
Section TVI D 5|) . in order to recover the NLer model re- 
sults, we need to expand the flow equations at large 
masses. Using the behavior of the threshold functions 
for large arguments (see Appendix EJ), we get |63| : 
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By making the change of variables: 
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we recover the /3-functions found in the framework of the 
NLer model at one-loop order (see Eq. (|49[l ). 



C. The large-TV expansion in d = 3 



As in the O(N) case (see Section lVlD 6|) . our equations 
allow to recover the critical exponents at leading order 
in 1/TV. We have computed rj and v for a large range of 
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FIG. 8: The exponents r\ and v as functions of TV in d — 3. 
The crosses represent our results and the full line the values 
obtained from the 1/TV expansion, Eqs. 1581 . The circles and 
error bars are the Monte Carlo results obtained for TV = 6 

na 

values of TV and have compared our results with those 
caculated perturbatively at order l/N 2 , Eq. J5SJ|. We 
find an excellent agreement — better than 1% — for 
v, for all TV > 10, see Fig. a domain of values of TV 
where one expects the 1/TV 2 results to be very close to 
the exact values. We also quote in Table Hxl our results 
and those obtained by the six-loop calculation for TV = 16 
and TV = 32. 



D. The determination of N c (d) 

Let us now interpolate between the results we have 
obtained around d = 2 and d = 4 and discuss, in par- 
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six-loop 


0.858(4) [155], 0.863(4)[132] 


0.0246(2)[132] 


32 


1/N[156] 


0.946 


0.0125 




NPRG 


0.950 


0.0134 




six-loop 


0.936(2)[155], 0.936(1)[132] 


0.01357(1) [132] 



TABLE IX: Exponents v and r\ computed from the 1/TV ex- 
pansion I156J. by our method (NPRG) and from the six-loop 
calculation |l32l 1 155|1 . 



ticular, the curve N c (d) that separate the regions of first 
and second order. 

We have computed N c (d) with our best truncation and 
with the cut-off function (|5rJj) . In Fig.|§l we give our re- 




_j i i i i i 

2.8 3 3.2 3.4 3.6 3.8 4 



d 

FIG. 9: The full line represents the curve N c (d) obtained by 
the three-loop results improved by the constraint TV c (d = 2) = 
2, Eq. 1561 . The crosses represent our calculation. 

suit (crosses) from d = 4 down to d = 2.7. We also indi- 
cate the improved three-loop results given by Eq. (|56*|> for 
comparison. The two curves agree very well. Note that it 
is probably a coincidence that the curves cross very close 
to d = 3. In this dimension, the NPRG method leads 
to N c (d — 3) ~ 5.1 and the improved three loop result: 
N c (d = 3) ~ 5.3(2). Let us emphasize that, within the 
NPRG method, the quan tity N c (d) is very sensitive to the 
order of the truncation 60], much more than the criti- 
cal exponents. This means that one probably should not 
consider the previous results as very reliable. In this re- 
spect, we recall the results obta ined by means of the six- 
loop calculation in d = 3 et al. |132| : N c (d = 3) = 6.4(4 ) 
and by means of the 4 — e expansion at five-loops |159| : 
N c (d = 3) ~ 6.1(6). 

Let us finally mention that, for the reason already ex- 
plained in Section iyi D 91 the field expansion we have per- 
formed at order <f> 10 forbids us to follow the chiral fixed 
point C+ in dimensions typically between d ^ 2.5 and 
d = 2.1 and thus to determine reliably the curve N c (d) 
in these dimensions. As in the O(N) case, this artefact 
could be overcome by keeping the full field dependence 
of the effective potential Uk (p, r) . 
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E. The critical exponents for N = 6 

As already said, for N = 6, the transition is either 
of second order or extremely weakly of first order. In 
both cases scaling should exist on a large domain of tem- 
perature. The critical exponents obtained with our best 
truncation are given in Table Ixl Note that v and r\ are 
computed directly while 7, (3 and a are computed using 
the scaling relations. Our results agree very well with 



Method 


a 


P 


7 


V 


V 


NPRG 


-0.121 


0.372 


1.377 


0.707 


0.053 


MCf22] 


-0.100(33) 


0.359(14) 


1.383(36) 


0.700(11) 


0.025(20) 



TABLE X: The exponents for N = 6 obtained from the NPRG 
— first line — and from the Monte Carlo (MC) simulation — 
second line. 

the numerical ones |22j . In particular, the error on v 1 
which is as usual the best determined exponent, is only 
1%. This constitutes a success of the NPRG approach 
from the methodological point of view. 



F. Conclusion 

Our method has successfully passed all possible tests. 
This gives us a great confidence in the reliability and the 
convergence of our results. We are now in a position to 
explore the physics of frustrated magnets in three dimen- 
sions. 



IX. THE PHYSICS IN d = 3 

We now tackle with the physics in three dimensions. 
Before embarking in this discussion, two points need to 
be clarified. The first concerns the existence of a fixed 
point for TV < N c (d = 3). The second one concerns the 
situation just below N c (d = 3). 



A. The search of fixed points for N < N c (d) 

Let us first discuss the search, within the NPRG 
method, of fixed points in d = 3 and for N < N c (d — 
3) ~ 5.1. We recall that, for this critical value of N, 
the two fixed points C + and C_ — see Fig. |21 — coa- 
lesce. This means that these fixed points — that can be 
followed smoothly in the (d, N) plane from the gaussian 
in d = 4 — cease to be real below this value. However, 
this does not imply the absence of other real fixed points. 
One has to test the existence of fixed points non trivially 
con nect ed with C + and C_ , as advocated by Pelissetto et 
al. |l55j| . We have thus searched such fixed points both 
by directly looking for zeroes of the /3-functions and by 
integrating numerically the RG flow — see below. After 
an intensive search, we have found no such fixed point. 



This result and its relation with that of Pelissetto et al. 
will be discussed in the following. 



B. The physics in d — 3 just below N c (d): scaling 
with a pseudo-fixed point and minimum of the flow 

In a fixed dimension d, the disappearance of the non- 
trivial fixed points C+ and C_, when N crosses N c (d), 
could seem to be an abrupt process: the two fixed points 
collapse and disappear. Actually, when extended to the 
space of complex coupling constants, this process is con- 
tinuous since the only change is that, when going from 
N > N c (d) to N < N c (d), the fixed points acquire a small 
complex part. This continuous character manifests itself 
as smooth changes of the RG flow that can be explained 
thanks to continuity arguments. 




N > N c {d) 




N just below N c (d) 



FIG. 10: Schematic representation of the flow a) for N above 
N c (d) and b) for N just below N c (d). For the sake of clarity, 
we have represented £' outside £ while it can be included in 
it. 

To understand the evolution of the RG flow as N is 
decreased, we need to consider the space of all coupling 
constants, i.e. the space such that to each point corre- 
sponds a microscopic Hamiltonian of a system. In this 
space, we focus on the subspace £ containing the repre- 
sentative points, at T = T c , of STA, STAR, Vjv,2, BCT 
and of all real materials studied experimentally and, more 
generally, of all systems of physical interest. Let us now 
describe qualitatively the change of the RG trajectories 
as N crosses N c (d). 

i) When N is larger than N c (d), there exists a true 
stable fixed point of the RG flow so that all trajectories 
emerging from £ eventually end on this fixed point, see 
Fig. HOal All systems exhibit scaling around the transi- 
tion and universality holds. 

ii) As already stated, when N is decreased slightly be- 
low N c (d), the fixed point C+ gets complex coordinates 
and looses its direct physical meaning. In particular, the 
flow no longer stops at a point, see Fig. II Obi Conse- 
quently, the correlation lengths of systems in £ do not 
diverge at T c . Strictly speaking, all systems undergo 
first order phase transitions. However two facts must 
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be noted. Firstly, all the trajectories emerging from £ 
are attracted toward a small region in coupling constant 
space, denoted by £' in Fig. II Obi Secondly, the flow in 
£' is very slow. 

From the second observation, we deduce that for all 
systems in £ the correlation lengths at the transition are 
very large — although finite — since they typically be- 
have as the exponential of the RG time spent around 
£', which is large. Therefore, the transitions are all ex- 
tremely weakly of first order for systems in £. We thus 
expect scaling behaviors with pseudo-critical exponents 
for all physical quantities, with the subtlety that this 
scaling aborts very close to T c where the true first order 
nature of the transitions shows up. 

As for the first observation — i.e. all trajectories are 
attracted toward a small region £' — , it allows to con- 
clude that all phase transitions are governed by a small 
region in coupling constant space and that, therefore, 
universality almost holds. In particular, the pseudo- 
critical exponents should be roughly the same for all sys- 
tems whose microscopic Hamiltonian corresponds to a 
point in £ . 

Let us study in greater detail the case where N is just 
below N c (d — 3). For such values of TV, it is reason- 
able to approximate £' by a point. The best approx- 
imation is clearly to choose the point in £' that mim- 
ics best a fixed point, i.e. the point where the flow is 
the slowest: the minimum of the flow |65j. To deter- 
mine this so-called "pseudo-fixed point" , Zumbach |65| 
has proposed to introduce a norm for the flow and to 
determine the point where the norm is minimum. He 
has performed this approach in the context of a NPRG 
equation (LPA of the Polchinski equation) where he has 
built the needed mathematical structures. He has shown 
that, when a minimum exists, pseudo-critical exponents 
characterizing pseudo-scaling can be associated with the 
pseudo- fixed point, in the same way that true exponents 
are associated with a true fixed point (see Appendix [D] 
for more details). 

A natural assumption to explain the pseudo-scaling be- 
haviors observed in real systems is that the minimum of 
the RG flow mimics a true fixed point even for values 
of N not very close to N c (d — 3). For the Heisenberg 
systems, this position has been advocated by Zumbach 
[65| and by the present authors [63| . 

Within our present approach we have confirmed the 
existence of a minimum of the flow, for values of N just 
below N c (d = 3), leading to pseudo-scaling and pseudo- 
universality |63| . By following this minimum we have 
confirmed that it persists down to N — 3 and have com- 
puted the associated pseudo-critical exponents, see Table 
IXII We also give in this table the exponents found by 
Zumbach within the LPA of the Polchinski equation for 
the same model [65j and recall thos e fou nd within the 
six-loop approach of Pelissetto et al. |155| . 

The values that we have obtained within our calcu- 
lation for the critical exponents are not too far from — 
some of — those found experimentally for group 2 of ma- 



Method 


Ref. 




P 


7 


V 


n 


NPRG 


f63] 


0.38 


0.29 


1.04 


0.54 


0.072 


LPA 


f65] 


0.11 


0.31 


1.26 


0.63 


0.0 


6-loop 


[155] 


0.35(9) 


0.30(2) 


1.06(5) 


0.55(3) 


0.08 



TABLE XI: The critical and pseudo-critical exponents for 
N = 3. a, (3 and 7 have been computed assuming that the 
scaling relations hold. The first line corresponds to our non- 
perturbative approach, the second to Zumbach's work. In the 
third line, we have recalled the six-loop results of Pelissetto 
et al. for comparison. 

tcrials, see Eq. Ij37(l . as well as those found numerically 
for the STA, Table IVl As usual, our truncation overes- 
timates r\ and thus, at fixed /3, underestimates v. It is 
remarkable that the values of the pseudo-critical expo- 
nents we have found at the minimum are in good agree- 
ment with those obtained within the six-loop approach. 
This strongly suggests that there is a common origin to 
these two sets of critical exponents. We shall come back 
on this point later. 



C. Scaling with or without pseudo-fixed point: the 
Heisenberg and XY cases 

Let us now argue that the preceding analysis, based 
solely on the notion of minimum, is too naive to give an 
explanation of the pseudo-critical behaviors in the phys- 
ically interesting cases. Let us also give a qualitative 
picture that supplements the concept of minimum. 

We have found that, when TV is lowered below N = 3, 
the minimum of the flow is less and less pronounced and 
that, for some value of N between 2 and 3, it completely 
disappears. Since several XY systems exhibit pseudo- 
scaling in experiments or in numerical simulations, this 
means that the concept of minimum of the flow does not 
constitute the definitive explanation of scaling in absence 
of a fixed point. One encounters here the limit of the 
concept of minimum of the flow. First, it darkens the 
important fact that the notion relevant to scaling is not 
the existence of a minimum but that of a whole region in 
coupling constant space in which the flow is slow, i.e. the 
(3 functions are small. Put it differently, the existence of a 
minimum does not guarantee that the flow is sufficiently 
slow to produce large correlation lengths. Reciprocally, 
one can encounter situations where the RG flow is slow, 
the correlation length being large so that scaling occurs 
even in absence of a minimum. The existence of a min- 
imum is thus neither necessary nor sufficient to explain 
pseudo-scaling. Second, even when the minimum exists, 
reducing the region £' to a point rules out the possibility 
of testing the violation of universality. For instance, one 
knows that for N — 3 universality is, in fact, violated, 
see Table IVl while a minimum of the RG flow is found. 
This feature cannot be reproduced by the unique set of 
exponents computed at the minimum. The opposite as- 
sumption, done first by Zumbach |65j and by the present 
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authors (6j|, was thus unjustified. 

Thus, even for very weak first order transitions, the 
beautiful simplicity of second order transitions is lost and 
the finite extend of the attractive region £' has to be 
taken into account. To be precise, one needs to define 
two subsets of £ and £': T> which is the region in £ 
leading to pseudo-scaling and TZ, the subset of £' which 
is the image of V in the RG flow, see Figs. Illal and ITTbl 
Let us now consider the characteristics of the flow when 
N is varied. 

Since for N > N c (d — 3) all the systems in £ undergo 
a second order phase transition, one expects — thanks to 
continuity arguments — that for N slightly below N c (d = 
3) , all systems in £ exhibit pseudo-scaling and thus that 
T> = £ . At the same time, £', the image of £ is almost 
point-like — see Fig llObl — and universality holds. 

As N is decreased below N c (d), two phenomena occur. 

i) While T> remains equal to £, the domain £', which 
is initially point-like, grows, see Fig. Illal This means 
that while pseudo-scaling should be generically observed, 
universality starts to be significantly violated: a whole 
spectrum of exponents should be observed, the size of £' 
providing a measure of this violation of universality. 

ii) For low values of N, the region T> leading to pseudo- 
scaling become smaller than £ , see Fig. lllb| For systems 
defined by initial conditions in T>, the correlation lengths 
are still relatively large but the pseudo-critical exponents 
can vary from system to system according to the size of 
TZ. For systems defined by initial conditions in £ but not 
in T>, the RG flow is always fast, producing small corre- 
lation lengths at T c . The corresponding systems undergo 
strong first order phase transitions. Moreover, as N de- 
creases, the flow in £' should become more and more 
rapid so that, for systems in £ , the correlation lengths 
at the phase transitions should decrease. The transitions 
are thus expected to become more strongly of first order 
for lower N. 




FIG. 11: Schematic representation of the flow a) for N below 
iV c (d) — N ~ 3 — and b) for N well below N c (d) — N ~ 2. 
For the sake of clarity, we have represented £' outside £ while 
it can be included in it. T> and TZ are represented in grey. In 
a), V = £ and TZ = £' . 

The precise values of N for which these changes of be- 
haviors occur as well as the shapes and extents of T>, 
TZ and £' can only be obtained from a detailed analy- 
sis of both the microscopic Hamiltonian and of the RG 



flow. However, independently of the details of the model 
under study, of the precise value of N c (d), etc, one ex- 
pects the following behavior: as N is decreased, a sys- 
tem that undergoes at large N a second order transition 
undergoes, for N just smaller than N c (d), a very weak 
first order transition governed by the minimum. Then, 
it should undergo a weak first order transition where the 
notion of minimum is no longer relevant and for which 
universality does not hold anymore. Finally, it should 
undergo a strong first order phase transition. In the 
spectrum of models studied numerically, it is easy to 
see that the STAR, VW, 2 and BCT models with XY and 
Heisenberg spins nicely obey this prediction. For N = 3, 
they all show scaling and the phase transitions should be 
very weakly of first order. However, their exponents are 
clearly incompatible with those of STA and with those 
associated with the minimum, see Tables Ivl and IXII The 
RG trajectories associated with these systems are thus 
expected to pass through TZ, but far from the minimum. 
One thus naturally expects that, when N is decreased 
down to N = 2, no scaling behavior is observed for these 
systems. This is indeed what is found in numerical simu- 
lations, see Table ITTT1 This strongly suggests that T> has 
shrinked between N — 3 and N — 2 and that N = 3 
corresponds to Fig. Illal and N = 2 to Fig. lllbl 



D. The integration of the RG flow for Heisenberg 
and XY systems 

In the previous section we have shown that the notion 
of minimum — or pseudo-fixed point — in the RG flow is 
neither sufficient nor necessary to explain the existence 
of scaling without a fixed point. For this reason, one 
has to resort to another method to study the physics of 
XY and Heisenberg frustrated magnets. In practice, we 
integrate numerically the RG flow around the transition 
temperature T c and determine the behavior of the phys- 
ical quantities such as the correlation length, the sus- 
ceptibility and the "magnetization" — defined as \f~k, 
see Eq. 1)1460 — as functions of the reduced temperature 
t r = (T~T C )/T C . 



1. Three difficulties 

Let us mention three difficulties encountered during 
the integration of the flow. 

First, in principle, in the absence of universality, we 
should study each system independently of the others. 
Thus, to correctly specify the initial conditions of the 
RG flow, we should also keep all the microscopic infor- 
mation relevant to the description of a given material. 
This program remains, in the most general case, a diffi- 
cult challenge since this would consist in keeping track 
of the lattice structure as well as of the infinite number 
of coupling constants involved in the microscopic Hamil- 
tonian. However this is, in principle, possible. Actually, 
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this has been done with much success for certain classes 
of magnetic systems and fluids de scribed by O(N) mod- 
els mostly within the LPA [l9l llQjj j. Our trun- 
cations — even the best one — are too restricted ap- 
proximations to reach this goal since this would at least 
require to keep the full field dependence of the poten- 
tial Uk (p, t) . We have thus used our flow equations to 
explain the generic occurrence of pseudo-scaling in frus- 
trated systems without trying to describe the behavior 
of a specific system. In practice, we have computed the 
correlation length, magnetization and susceptibility us- 
ing a simplified version of our truncation keeping only the 
potential part expanded up to order eight in the fields, 
a field-independent field renormalization and discarding 
all the current-terms involving four fields and two deriva- 
tives. We have checked that this ansatz leads to stable 
results with respect to the addition of higher powers of 
the fields and inclusion of current terms. 

Second, the truncations we have considered do not al- 
low to determine accurately the critical temperature. In- 
deed, in our approach we perform a local description of 
the potential around the nontrivial minimum Eq. <|146|) . 
For a second order phase transition this does not matter 
since the nontrivial minimum, when it exists, is always 
the true one. However, for a first order transition, the 
zero-field configuration, i.e. with <px = (f>2 = 0, plays a 
crucial role. In effect, in this case, the transition temper- 
ature precisely coincides with the temperature at which 
the energy at the nontrivial minimum and at the zero- 
field configuration are equal. Since we cannot expect that 
our truncation describes accurately the potential around 
the zero-field configuration, we are not able to compare 
the energy of this configuration with the energy of the 
nontrivial minimum and to determine the transition tem- 
perature accurately. We discuss in more details this point 
in Appendix El and show that, for a weak first order tran- 
sition, this fact should not bias significantly our analysis. 

The third difficulty encountered in the integration of 
the flow is that, in the absence of universality, the tem- 
perature dependence of the physical quantities relies on 
the precise temperature dependence of the microscopic 
coupling constants. We have used several ansdtze for the 
temperature dependence of the coupling constants and 
have observed that, although it could be important for 
the details of the results, it does not affect much the gen- 
eral conclusions. Thus, we illustrate our results with the 
simplest ansatz consisting in fixing all the couplings to 
temperature-independent values and by taking a linear 
temperature dependence for k at the lattice scale: 

n k=A =a + bT . (159) 

For each temperature, we have integrated the flow equa- 
tions and have deduced the ^-dependence of the physical 
quantities, such as the "magnetization", the correlation 
length, etc, around T c . The different coupling constants 
parametrizing the initial condition of the flow have been 
varied to test the robustness of our conclusions. This has 
allowed us to establish the following facts. 
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FIG. 12: Log-log plot of the magnetization m and of the 
correlation length £ for N — 3 as functions of the reduced 
temperature t r . The straight lines correspond to the best 
power law fit of the data. 



2. The Heisenberg case 

• For N — 3, we can find initial conditions of the flow 
such that for a wide range of reduced temperatures - 
up to four decades — the physical quantities behave as 
power laws. From an experimental viewpoint, this is all 
what is needed since scaling has been found on temper- 
ature ranges that are even smaller. The kind of pseudo- 
critical behaviors we find is illustrated on Fig. 1121 

• Varying the initial conditions of the flow, we observe 
that this phenomenon happens in a wide domain of the 
coupling constant space. This corresponds to the domain 
V previously defined, see Fig. Illal 

• Within T>, the pseudo-critical exponents vary 
smoothly: (3 varies typically between 0.27 and 0.42 and 
v between 0.56 and 0.71. These are only typical values 
since it has been impossible to explore the whole space 
of coupling constants. Since for (3 ~ 0.27 one can find 
v ~ 0.56, the exponents of group 2 are satisfactorily re- 
produced, see Tables ITVl and Ivl This shows in particular 
that there exists, in T>, a set of "microscopic" coupling 
constants that lead to the behavior observed in group 2. 

• It is easy to find initial conditions leading to pseudo- 
critical exponents in good agreement with those obtained 
in the six-loop calculation, see Table IVIII Actually, a 
whole set of initial conditions lead to exactly the same 
(pseudo-) critical exponents as those found at six-loop 
(3 = 0.30(2), v — 0.55(3). This corresponds to the region 
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of the minimum of the flow, see Table IXll 

• In contrast, we have not found initial conditions of 
the RG flow reproducing correctly the critical exponents 
of group 1, of STAR, Y^,2 and BCT as well as negative 
values for r\. This can originate i) in the overestimation 
of r\ produced by our truncation of Tk in powers of the 
derivatives at order d 2 , Eq. I|145|) . ii) in the impossibility 
to sample the whole coupling constant space, in) in the 
too simple temperature dependence of ka that we have 
considered, see Eq. H159|) . 

• For a given value of one exponent, it is possible to find 
several values for the other exponents. Thus we expect 
to find systems sharing for instance almost the same ft 
but having quite different values for v and 7. 

• At the border of T>, the temperature ranges over 
which power laws hold become smaller and smaller. In 
a log-log plot, the independence becomes less and less 
linear and the pseudo-critical exponents more and more 
sensitive to the choice of T c made for the fit. Finally, 
outside T>, no more power-law behavior is observed. 

• When we go from N = 3 to TV = 4, we have observed, 
as expected, that T> becomes far wider and that the power 
laws hold gcncrically on larger temperature ranges. This 
is consistent with our discussion of Section HX CI Recip- 
rocally, and as also expected, when going from N = 3 
to N — 2, T> becomes much smaller and the power laws 
hold generically on smaller temperature ranges. Let us 
discuss this point in greater detail now. 



3. The XY case 
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FIG. 13: Log-log plot of the of the magnetization m and of 
the correlation length £ for N = 2 as functions of the reduced 
temperature t r . The straight lines correspond to the best 
power law fit of the data. The power-law behavior observed 
far from the critical temperature breaks down for small t r . 
The behavior of the correlation length at small t r is an artefact 
of our truncation, see Appendix IeI 



• For N — 2, one observes qualitatively the same type 
of behaviors as for N = 3. However, as predicted above, 
T> is smaller and the power laws hold at best only on 
two decades of reduced temperature, which is consistent 
with what is observed experimentally. This is illustrated 
in Fig. 1131 where we have represented log-log plots of the 
magnetization and correlation length as functions of the 
reduced temperature. Note also the surprising behavior 
of the correlation length that increases at small reduced 
temperature (see Appendix IEI for an explanation of this 
phenomenon.). 

• Within T> the exponents vary on the intervals: 0.25 < 
ft < 0.38 and 0.47 < v < 0.58. 

• We find initial conditions leading to exponents close 
to those of group 2 (for Ho and Dy, see Table ITTfl : 
ft = 0.38, v = 0.58, 7 = 1.13. These results are quite 
stable with respect to changes of microscopic parameters. 
This is in agreement with the stability of (3 in group 2. 
Interestingly, these initial conditions correspond to small 
p, in our truncation Eq. i|147f) . i-s- to initial conditions 
close the 0(4)-invariant line: jl = 0, see Fig. |2 where 
the 0(4) fixed point is denoted by V [213] . Thus, during 
a large part of the flow, the trajectory remains close to 
the 0(4) fixed point before bifurcating away from this 
point. This is perhaps the reason why the value of (3 
for materials of group 2 is close to that associated with 



an 0(4) behavior — fto(4) = 0.382 — a fact that has 
been repeatedly noticed by experimentalists. Note how- 
ever that the other exponents are not close to the 0(4) 
values: ^0(4) = 0.738, 7o(4) = 1-449. 

• We also easily find initial conditions leading to 
ft = 0.25, corresponding to group 1, essentially com- 
posed of STA systems. The power laws then hold on 
smaller ranges of temperatures and the critical exponent 
ft is more sensitive to the determination of T c and to the 
initial conditions. For such values of ft, we find that v 
varies between 0.47 and 0.49, which is somewhat below 
the value found for CsMnBr3, see Table Q] 

• The two previous points suggest that both helimag- 
nets — such as Ho or Dy — and STA — such as CsMnBr3 

— can be described by the same field theory but with ex- 
ponents at the two ends of the spectrum. It is actually 
also possible that helimagnets display a different kind of 
physics because of the presence of long range inter actions 
or because of the presence of surface effects |196| . 

• As in the N — 3 case, we can easily find initial condi- 
tions leading to pseudo-critical exponents close to those 
found in the six-loop calculation, Table lvTTI For instance, 
for initial conditions leading to ft = 0.33, we find typi- 
cally: v = 0.56 and 7 = 1.07. 

• As in the Heisenberg case, we have not been able to 
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find initial conditions of the RG flow leading to negative 
values of rj. 

Let us now comment our results. 



4- Comments 

The main feature of the physics of Heisenberg and XY 
frustrated magnets — scaling behaviors without univer- 
sality — is reproduced, at least qualitatively and, to some 
extent, quantitatively. This behavior finds a natural ex- 
planation: there exists a whole domain T> in the space of 
coupling constants such that the RG trajectories starting 
in T> are "attracted" toward a region 1Z where the RG 
flow is slow so that there is pseudo-scaling. Since 1Z is 
not reduced to a point, there exists a whole spectrum of 
exponents and not a unique set. The occurence of strong 
first order phase transitions, that are observed in some 
materials and simulated systems, is explained by the RG 
trajectories starting out of V. 

Let us now stress that since universality is lost, the de- 
termination of the precise pseudo-critical exponents asso- 
ciated with a given material or system is obviously more 
difficult than the determination of the usual — universal 
— critical exponents characterizing a second order phase 
transition. As already said, computing them would in- 
deed require to know precisely the microscopic structure 
of the materials or systems studied — providing the ini- 
tial conditions of the RG flow — and to take into account 
the full field-dependence of the potential U% (p, t) . 

X. POSSIBLE TESTS OF OUR SCENARIO 

There are several tests that can be performed both 
experimentally and numerically to confirm our proposals. 
Let us start by the Heisenberg case. 

• It is not clear, up to now, whether the materials of 
group 1 — VCI2 and VBr2 — are really three-dimensional 
Heisenberg STA, at least for a temperature range wide 
enough to measure exponents. It would be very inter- 
esting to re-study these materials and to measure all ex- 
ponents for each of them. This could allow to confirm 
experimentally the absence of universality. 

• Since we predict that they can be violated, there is 
clearly a need to check the scaling relations as well as 
the negativity of rj. The experimental determination of 
the exponents 7 and v for the two groups of Heisenberg 
materials is still much too poor. It is also necessary to 
have an estimation of both the systematic and statistical 
errors to strengthen our conclusion on the negativity of 
rj. Let us however recall that the first order nature of the 
transitions in Heisenberg systems is likely to be much 
weaker than in XY systems. Thus the violations of both 
the scaling relations and the positivity of rj should be 
much more difficult to prove experimentally in this case. 

• It remains mysterious why, in CsMn(Bro.igIo. 81)31 
such strange values of the exponents 7 and v have been 



found, see Table IIVI As we have already discussed in 
point i) of Section llll G II we find unconvincing the ar- 
guments proposed in |141| to explain them. Remeasuring 
these exponents could provide accurate results for 7 and 
v from which universality and the negativity of 77 could 
be tested. 

• Most probably STAR and the Va^ model undergo 
both first order transitions since rj is found negative for 
these models. It would be extremely interesting to study 
a sequence of models that interpolate between STA and 
STAR to see how the effective exponents change with the 
deformation of the model. 

• We have already noticed that the exponents found 
in the six-loop calculation are very close to the pseudo- 
critical exponents found at the minimum of the RG flow 
in our approach. It is important to know if this is just an 
accidental coincidence or if they correspond to the same 
fixed point, real in one approach complex in the other. 

Let us now discuss the XY case. Most of the points 
discussed in the Heisenberg case can be transposed here: 
necessity to check the scaling relations and the positiv- 
ity of 77, possibility to interpolate between the STA and 
STAR. Here, however, we are in a better position to ob- 
tain conclusive results since the transition is expected to 
be more strongly of first order. 

• A better determination of v in CsMnBr3 would help 
to confirm that 77 is indeed negative. We also expect to 
have a weaker universality and thus a faster change of 
the exponents as the microscopic details of the model 
are varied. In particular, a precise determination of a in 
the different materials of group 1 could lead to incom- 
patible exponents — they are up to now only marginally 
compatible — and would give a direct proof of the lack 
of universality. 

• On the numerical side, the sequence of models inter- 
polating between STA and STAR should lead to rapidly 
varying exponents. Thus the lack of universality in this 
case should be much simpler to prove numerically than 
in the Heisenberg case. For STA, it would also be ex- 
tremely interesting to determine 77 independently by the 
two scaling relations 77 = 2(3 /v — 1 and rj — 2 — 7/1/. As 
far as we know 77 has mainly been determined using 7/f. 
According to our scenario the two determinations should 
not coincide. However, they are probably both negative. 



XI. CONSEQUENCES FOR PERTURBATION 
THEORIES 

Frustrated magnets represent a unique controversial 
example of systems for which almost all the possible 
perturbative and nonperturbative approaches have been 
used, sometimes with a very high precision. This allows 
to draw several conclusions about the relative predictive 
power of these different methods. Firstly, it appears that 
the low-order results obtained within the NLcr or GLW 
models fail to correctly describe the physics in three di- 
mensions. Indeed, we recall that the one-loop result of 
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the NLer model predicts a second order phase transition 
with a 0(4) behavior while the GLW approach leads to 
first order phase transitions for all values of N smaller 
than 21.8. Secondly, frustrated magnets probably pro- 
vides the first example where high-order perturbative re- 
sults is questionable. We shall now discuss the status of 
the various perturbative approaches at the light of our 
results. 



A. The NLer model approach 

Let us first consider the NLa model approach, focus- 
ing on the Heisenberg case since it is notorious that this 
approach does not work for XY spins. The very likely 
existence of a line N c (d) going from d = 2 to d = 4 
confirms what has been already anticipated in Section 
Ivl the predictions based on this approach are incorrect 
as for the physics in d = 3. Indeed, the shape of this 
line implies that the 0(4) fixed point predicted in the 
Heisenberg case — that likely exists at all orders of per- 
turbation theory — very probably disappears between 
two and three dimensions. Actually, following this fixed 
point, that we call C+ for an obvious reason, from d = 2 
with the simplest </> 4 -like truncation, we have found sev- 
eral interesting features. First, infinitesimally close to 
d — 2, we find that C+ is characterized by an exponent v 
of the 0(4) universality class. Second, as d is increased, 
the exponent v associated with C+ becomes more and 
more different from that characterizing an 0(4) transi- 
tion. Third, we find that an unstable fixed point, C_, 
shows up in a dimension d > 2. As the dimension is 
further increased, the fixed points C + and C_ get closer 
together and eventually coalesce in a dimension less than 
three. This phenomenon is illustrated in Fig l|14|l in the 
case of the 0(3) x 0(3) model at the lowest order in the 
field expansion [6^|. 

The collapse of the fixed points for different values of 
N generates the curve N c (d). This curves is well known 
from the perturbative expansion performed around four 
dimensions. Within our approach, this curve can be fol- 
lowed when the dimension is lowered down to d = 2. 
There, for a given — low — value of N, the curve N c (d) 
provides the value of d c (N) for which the stable fixed 
point obtained within the NLer model approach collapses 
with another — unstable — fixed point. Since this un- 
stable fixed point is not found in the low-temperature 
perturbative expansion we therefore obtain here a non- 
perturbative solution to the breakdown of the NLer model 
approach. For N = 3, one gets d c ~ 2.8. Note that ob- 
taining an accurate determination of the dimension d c 
where the fixed points collapse would require to consider 
better truncations in fields of than those we have con- 
sidered. However, as already explained in the O(N) case, 
see Section lYl D 91 the stable fixed point coalesces in this 
case with one of the multicritical points. Thus it is im- 
possible, within our truncation, to follow it smoothly for 
2-1 d ~ 2.5. With our best truncation, we are anyway 
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FIG. 14: The stable — C+ — and unstable — C_ — fixed 
points as functions of the dimension d. The fixed points are 
parametrized by the quantity k* which is proportionnal to 
the inverse temperature of the NL<r model. The fixed point 
C- appears in a dimension 8/3 and collapses with the stable 
fixed point C+ in d ~ 2.83. 



able to give an estimate of this dimension: d c ~ 2.6 — 2. 7 
which fits well with the results of Pelissetto et al. |l56j |. 
see Fig. 

Frustrated magnets thus provide a situation where 
there is a manifest breakdown of the low-temperature 
expansion of the NLer model. This is not the first oc- 
curence of such a breakdown. The case of the two- 
component nonfrustrated 0(2) system has already exem- 
plified the inadequacy of the low-temperature expansion 
to explain the existence of a phase transition for these 
systems in two dimensions, the Berezinskii-Kosterlitz- 
Thouless phase transition |52t l53j . There is however an 
important difference between the case of XY nonfrus- 
trated spins and that of Heisenberg frustrated spins. In- 
deed, in the former case, the low-temperature expansion 
performed on the corresponding 0(2) NLer model leads 
to a free theory to all orders in the temperature T in 
d > 2. This result is however known to be incorrect for 
XY spins themselves or for the systems that belong to the 
same universality class — like 4 He — that both undergo 
a phase transition in d > 2. In this case, the unability 
of the low-temperature expansion to correctly describe 
the physics makes no doubt and one i s invited to turn 
to other methods: Coulomb-gas |191| or spin-vortices 
[52l |5^| formulations in two dimensions or GLW model 
approach in three dimensions. On the contrary, in the 
case of Heisenberg frustrated spins, the low-temperature 
expansion leads to a nontrivial behavior — due to the 
nonabelian character of the 50(3) group — so that the 
inadequacy of the low-temperature perturbation theory 
is not so obvious. 

It remains to understand the very origin of this failure 
of the low-temperature perturbation theory. In the case 
of XY nonfrustrated spins, it clearly lies in the existence 
of nontrivial topological configurations, called vortices, 
that are not taken into account in a low-temperature ex- 
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pansion. In the case of Heisenberg frustrated magnets, 
the influence of nontrivial topological configurations on 
the phase transition in three dimensions has also been 
invoked (see Section V-A). It remains however to con- 
firm that these configurations indeed play a fundamental 
role and to know, for instance, if they are responsible 
for the first order character of the transitions in three 
dimensions. 

This is a delicate question. Indeed, whereas the pertur- 
bative approach to the NLer model misses topologically 
nontrivial configurations, the GLW and effective aver- 
age action approaches are very likely sensitive to such 
vortices. In effect, both approaches correctly reproduce 
the physics of three-dimensional XY nonfrustrated spin 
systems at the tra nsition which is very likely driven by 
vortices |l97t Q.98J. However, within these approaches, 
it is still not clear how the vortices are taken into ac- 
count. Therefore, disentangling vortices and spin-waves 
and understanding the respective role of each kind of ex- 
citation within the phase transition remains a theoretical 
challenge. 



B. The GLW model approach vs the NPRG 
approach 

We now discuss the relationship between the weak- 
coupling results obtained within the GLW model ap- 
proach — in particular, the six-loop computation — and 
our results. A natural question arises: how is it possi- 
ble to reconcile these results together and what does this 
imply for the different approaches? 

We have noted an important fact: the critical expo- 
nents found for N = 3 in the six-loop calculation and 
in our approach — at the minimum of the flow — are 
very close (see Sections IIX Bl and IIXD|I . We have also 
found very close exponents for N — 2 (see Section HX Dp 
with the only difference that there is no minimum in the 
flow in this case. This is a rather strong indication that 
the two sets of exponents have a common origin. This 
leads us to formulate some proposals to reconcile the two 
approaches. 

The first one is that the fixed point that appears as real 
in the six-loop calculation and complex in our approach 
is, actually, a complex one. This would mean that the 
computations of Pelissetto et al. and Calabrese et al. is, 
actually, not converged as for the nature — real or com- 
plex — of this fixed point whereas it is almost converged 
as for the exponents. We shall not speculate too much 
about the origin of this — hypothethical — failure of the 
weak-coupling approach. Let us mention again however 
that the perturbative series obtained in the case of frus- 
trated magnets appear to be rather particular since the 
critical properties deduced from them strongly depend 
on the order of the series. We recall that there is no 
nontrivial fixed point for N = 2 and TV = 3 up to three 
loops; they only appear at five loops. Also, the six-loop 
results has been obtained in a region where the pertur- 



bative expansions are not Borel summable |l55j . It is 
clear that this question deserves further investigations. 
Frustrated magnets could appear as the first exemple of 
a breakdown of a weak-coupling perturbative analysis. 

The second proposal is that, reciprocally, within the 
NPRG method, the lack of fixed point in the XY and 
Heisenberg cases is due to artefacts of the truncation 
in fields and/or derivatives. Only the recourse to other 
kinds of expansions of the effective action Tk — involving 
either the full function Uk (p, t) or the full momentum de- 
pendence — could lead to unambiguous statements. In 
this respect, we however recall that the LPA approach of 
Zumbach, that involves the full field-dependence of the 
potential, has led to no fixed point for N = 2 and N = 3. 



XII. CONCLUSION AND PROSPECTS 

On the basis of their specific s y m metry breaking 
scheme, it has been proposed |o7t l68l 17(1 llZl IB! Il30| 

that the critical physics of XY and Heisenberg frus- 
trated systems in three dimensions could be characterized 
by critical exponents associated with a new universality 
class. From this point of view, the study of frustrated 
magnets has been rather disappointing, the experimental 
and numerical contexts excluding such an hypothesis. At 
the same time, the phenomenology of frustrated magnets 
has displayed a novel kind o f crit ical behavior — generic 
scaling without universality |l!8| — requiring the use of 
new theoretical approaches. 

H igh-o rder perturbative calculations in d = 3 |l3ll 

Il32l Il55j provide an explanation to the lack of universal- 
ity in frustrated magnets: the focus character of the fixed 
point induces spiral-like R G tra jectories from which, ac- 
cording to Calabre se et al. | l32l , follows varying effective 
critical exponents |l3ll Il32j . We have, however, under- 
lined several drawbacks of this explanation. The major 
one lies in its lack of naturalness: a fine-tuning of initial 
conditions seems to be necessary to match with the phe- 
nomenology. Another drawback of the perturbative ap- 
proach is that, being restricted to investigate the physics 
in three dimensions, it cannot provide a general picture 
of what happens between two and four dimensions. In 
particular it provides no explanation to the failure of the 
NLct model approach. 

Within the framework of a NPRG approach, the 
generic and nonuniversal scaling finds a natural explana- 
tion in terms of the slowness and "geometry" of the flow. 
This method also explains the mismatch between the dif- 
ferent perturbative approaches by means of a mechanism 
of annihilation of fixed points in a dimension between two 
and three that invalidates the low-temperature perturba- 
tive approach performed from the NLct model. As said 
along this article, more work, in particular the recourse 
to ansatz involving the full field-dependence or full mo- 
mentum dependence of the effective action, is probably 
necessary to completely understand the situation. This 
includes the clarification of the relation between our ap- 
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proach and the six-loop results. However, the main fea- 
tures of frustrated magnets appear now to be well de- 
scribed. 

The main result of this article is the explanation of the 
generic character of weak first order phase transitions in 
frustrated magnets. Being given the closeness between 
these systems and others systems — see the Introduction 
— it is natural to speculate about the degree of generality 
of this phenomenon. 

Within our approach, the generic character of the weak 
first order phase transition appear to be strongly related 
to the proximity of the number of components N of the 
system under study with N c (d = 3). For frustrated sys- 
tems, it appears that this value is of the same order than 
the physically relevant values of N, N = 2 and N = 3. 
This could be a very specific property of the frustrated 
systems. We now argue that, on the contrary, this prop- 
erty is likely to be common to many other systems. 

Let us recall that the line N c (d) corresponds to the 
collapse of two fixed points, one of them governing the 
phase transition. This phenomenon cannot happen in 
theories with only one (f> 4 coupling constant (i.e. in O(N) 
models) since, in this case, there is only one fixed point 
apart from the gaussian. However, for theories with c 
coupling constants, we expect 2 C perturbative — real or 
complex — fixed points in d = 4 — e since, at one-loop, 
the /3-functions are quadratic in the coupling constants. 
When the number of components of the field is varied, 
these fixed points move in the coupling constant space 
and it is generically observed that they meet and collapse 
for some critical value N c (d). Many examples are now 
known in the literature. Let us review some of them. 

Let us first consider the generalization of the model 
studied in this article consisting in P orthonormal TV- 
component vectors. It has a P-dependent critical value 
of N given at one-loop by l6al: N c ~ 10P. For P = 3, 
one finds at two loop order |l56(: 

iV c (4-e) = 32.49 - 33. 72e . (160) 

In the Abelian Higgs model coupled to a TV-component 
scalar field, relevant to su perconduc tors , N c [d) is found 
at two- loop order to be [lIlM l200| : 

N c (d = 4 - e) = 182.9 - 242. 7e . (161) 

In a SU (2) gauge mode l coupled to bosons, it is given 
at two loop order by |20G| : 

N c {d = 4- e) = 718 - 990.8e . (162) 

In a 0(p) gauge theory coupled to N scalar fiel ds (in 
the vector representation) it is given at one-loop by |20l| : 
N c ~40p. 

In all these examples, we observe that N c (d) decreases 
very steeply when d decreases. This is in line with our 
expectation that large N and small d favour continuous 
phase transitions. In particular, as far as we know, in all 
NLcr models relevant to systems whose order parameter is 



continuous, a stable fixed point is found in d = 2+e for all 
N > 2. This is in particular the case for the NLcr model 
supposed to descri be th e physics of the abelian Higgs 
model in d = 2 0,1202 • This suggests that N c (d = 2) 
is always smaller or equal to 2. It is interesting to notice 
that this bound is probably reached in frustrated systems 
jl56j , see Section IVB3I It is thus extremely probable 
that in many systems the curve N c (d) has a similar shape 
as the one found in frustrated systems, see Fig. This 
suggests that many systems could exhibit weakly first 
order transi tions in d = 3 without any fine-tuning of 
parameters |214| . The effective average action method 
should be ideally suited to study these situations. 
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APPENDIX A: THE POSITIVITY OF THE 
ANOMALOUS DIMENSION 

In this appendix, we sketch the proof showing that 
the anomalous dimension r\ must be positive in a second 
order phase transition if the underlying theory is given by 
a usual 4 -like GLW theory. This excludes, for instance, 
theories involving gauge fields or replica field theories of 
disordered systems using the formal N — > limit. The 
argument goes as follows. On one hand, using the Kallen- 
Lhemann decomposition, it is possible to prove that the 
field renormalization Z is positive and less or equal to 
one: < Z < 1 |l!9j |. On the other hand, using the 
RG equations, it is possible to show that, around the 
fixed point describing the second order phase transition, 
Z behaves with the RG scale p as: 

Z{p) ~ M " (Al) 

when ji — ► 0, which corresponds to the long distance — 
i.e. critical — physics. By combining these two results 
we find i] > 0. 

APPENDIX B: THE INVARIANTS OF THE 
SYMMETRY GROUP 

We show, in this appendix, that all field combinations 
invariant under O(N) x 0(2) can be rewritten in terms 
of the two invariants p and r introduced in Eq. (|52|l and 
given by p = Tr( *$$) and r = ±Tr (*$$ - i p/2) 2 . This 
property is a priori nontrivial since we can easily build 
an infinite number of invariants by considering, for in- 
stance, Tr ('$.$)" for any value of n or det (*$.$). The 
result is easily obtained by using the properties of the 
characteristic polynomial of the square matrix X: 

P X (A) = det(X - Al) . (Bl) 
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In the case of a two by two matrix, the characteristic 
polynomial reads: 



P X (X) = A 2 - A TrX + det X . 



(B2) 



The Cayley-Hamilton theorem states that any matrix is 
a root of its characteristic polynomial: 



Px{X) = 



(B3) 



Applying this last result to the two by two matrix 
we get: 

(*$.$) 2 - Tr (*$.$) + det ('$.$) = . (B4) 

By taking the trace of this equation, we see that 
det (*$.<&) = p 2 /4-r. Moreover, if we multiply Eq. !|B4|) 
by *$.<!> and take the trace of this equation, we observe 
that Tr (*$.<I>) 3 can be expressed in terms of Tr ('$.<f>) 2 , 
Tr and det (*$.$) which, themselves, can be ex- 

pressed in terms of p and r. By iteration, we can show 
that all 0(N) x 0(2) invariants can be expressed in terms 
of p and t. This property can be generalized to the 
O(N) x O(P) model (with N > P), which admits P 
independent invariants. 
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with r(y) being the dimensionless cut-off: 



r(y) 



R k (yk 2 ) 

Zyk 2 



(C2) 



(C3) 



APPENDIX C: THE THRESHOLD FUNCTIONS 

We discuss in this appendix the different threshold 
functions I, m and n appearing in the flow equations, 
which encode the nonperturbative properties of the the- 
ory. We consider here a general case, where the thresh- 
old functions depend on three arguments. For particular 
truncations — for instance that discussed in the O(N) 
vectorial model — it may happen that some of these ar- 
guments are vanishing. In such case, we do not write the 
associated argument so that, for instance, i^ (w,0,0) is 
denoted by l^w) 



We recall that the tilde in dt means that only the t 
dependence of the function Rk is to be considered. As a 
consequence, we should not consider the t-dependence of 
the coupling constants to perform this derivative. There- 
fore, in the preceding equations: 



d t P = 



dRk d 



P 



dt dRk 
= -viwiv) + 2yr'(y))- 



(C4) 
(C5) 



1. Definitions 



The threshold functions arc defined as: 



Now, threshold functions can be expressed as explicit in- 
tegrals if we compute the action of dt- To this end, it 
is interesting to notice the equality: d y dtPi — dtd y Pi 1 so 
that: 



^1,^(^1)^2,0) 



1 



+ wi)"i(P 2 + w 2 ) n2 

(Cla 



mi un2 (w 1 ,W2,a) 



d t d y r(y) = - v (r(y) + yr'(y))-2y(2r'(y)+yr"(y)) (C6) 



We then get: 
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Once a regulator r(y) has been chosen, the thresh- tions. When calculating the flow equation for the cou- 
old functions can be computed numerically and, in some pling constants related to the potential part, a single 



cases, analytically. 



function I appears: 



2. Substitution rules 



We give here the rules which relate the different inte- 
grals appearing in the calculation to the threshold func- 
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For the coupling constants associated with the derivative 
terms, two more functions appear: 
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Notice that the powers of k and Z appearing in the pre- 
ceding expressions are chosen so that when the flow equa- 
tions are reexpressed in terms of dimensionless renormal- 
ized quantities, there is no explicit dependence on these 
parameters. 



3. Universal values of the threshold functions 

For particular arguments, the threshold functions take 
values independent of the choice of the regulating func- 
tion r(y). This is particularly important when we extract 
the first coefficients of the perturbative (3 functions out of 
the nonperturbative ones, since the former are universal. 
From Eq. I|C7() we can compute the value of Z 2n (0,0, a) 
which enters in the (3 function for the coupling constant 
of the GLW model around four dimensions: 



C (0,0,o) = -n / dy 



r'(y) 



a (1 + a + riy))^ 1 

= [(l + a + r(y))-X (C13) 
f 

~~ (1 + a) n ' 

The last equality follows from the asymptotic behaviors 
of r(y) that are given by Eqs. ||HTJ and 



lim r(y) = 
lim y r(y) = 1 



(C14) 



and are independant of the actual form chosen for r(y). 
Similarly, one finds: 



^ n (0,0,a) = l 



(C15) 



Also, the threshold function m2 2 (w,0, 0) takes a uni- 
versal form for large argument u>, which enters in the (3 
function of the temperature in the NLct model around 
two dimensions. Using Eq. I|C7(1 . one gets: 



lim w 2 m\ 2 {w, 0, 0) 



dy d y 



1 + r(y) + yr'{y) 



1 + r(y) 
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4. Threshold functions from the theta cut-off 

For certain regulating functions r(y), it is possible to 
compute analytically the threshold functions. Using such 
regulating functions is very helpful in practice and sim- 
plifies considerably the numerical procedures. In this sec- 
tion, we give the threshold functions associated with the 
theta cut-off, see Eq. (|86ll . One has, for a = 0: 
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APPENDIX D: THE MINIMUM OF THE RG 
FLOW 

In this appendix, we describe in more details the no- 
tions of pseudo-fixed point and of minimum of the flow. 
We then explain how these ideas have been implemented 
in practice to determine effective exponents for very 
weakly first order phase transitions. 

As described previously, the RG flow equations for 
STA with a large number of spin components (N > 
N c (d)) admit two fixed points. When N is decreased 
slightly below N c (d), the two fixed points acquire a 
small complex part and loose their direct physical rel- 
evance. Strictly speaking, there is no more attractor in 
the real coupling constant space but the flow remains 
sensitive to the presence of complex fixed points. Zum- 
bach [oil l65l log proposed that a particular point, the 
minimum of the flow, should mimic to some extent the 
behavior of an attractor. This point is defined as the lo- 
cation, in coupling constant space, where the flow is the 
slowest and the quantity: 



Mi9i}) = JEA 2 



(Dl) 



— where (3i are the /3-functions for the different coupling 
constants gi — is minimum. Let us stress on few prop- 
erties of the minimum of the flow: 



where, again, we have used the asymptotic behaviors of 

r(y), Eq. ifHT4) . 



• when a true fixed point exists, A({g*}) — and, in 
this case, the minimum is a fixed point. 
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• When two fixed points annihilate, we are left with 
a single minimum of the flow sitting right at the 
position where the fixed points have collapsed. 

• For trajectories getting close to the minimum, the 
RG time spent in its vicinity is large and so is the 
correlation length. 

We therefore see that a minimum shares some features 
with a true fixed point. 

One easily obtains the equation characterizing a mini- 
mum: 



8A 
dgi 



E 

3 



with 



M, 



1,3 



Mijfc = 



dgi 



(D2) 



(D3) 



Under the assumption that the minimum of the flow 
mimics correctly the attractor of the trajectories, it is 
natural to compute the critical pseudo-critical exponents 
in the standard way. The anomalous dimension is ob- 
tained by evaluating i]({g.i}) at the minimum of the flow 
and v by diagonalizing the matrix Mij at this point. It 
is important to notice that the pseudo-critical exponents 
thus obtained are invariant under reparametrization of 
coupling constants, as it should be, since Eq. (|D2|I trans- 
forms as components of a vector. 



APPENDIX E: THE DISCONTINUOUS 
CHARACTER OF THE PHASE TRANSITION 

In this appendix, we discuss in more details the prob- 
lems that we encounter in our description of the first 
order phase transition that occurs in frustrated magnets. 
We also explain the surprising increase of the correla- 
tion length observed at small reduced temperatures (see 

Fig. ni)). 

To this end, let us discuss the following toy model of 
first order phase transition. We consider a scalar Z2- 
invariant model characterized by one field <f> and by a 
"0 4 - 6 " potential: 



U(6) = r— - — + — 
KV> 2 2 6 



(El) 



As usual, we assume that r varies linearly with the tem- 
perature. 

For low reduced temperatures — small r — the poten- 
tial has a local minimum for cj> = and a global minimum 
for: 



/Min 



(r)) 2 = l + VT^ 



(E2) 



so that the system exhibits a spontaneous magnetization, 
see Fig. 1151 When the temperature — r — is increased, 
the energy difference between the configurations = 




FIG. 15: Shape of the potential Eq. iEll for different tem- 
peratures, i.e. different values of the parameter r. The plain 
lines correspond to r = r c and r = r* , while the dotted lines 
correspond to different generic values or r. 



and 4> = <fi Mln (r) decreases and eventually vanishes for 
r — r c = 3/4 which defines the critical temperature. For 
r larger than r c , the ground state of the system is given 
by the configuration <f> — so that the system has no 
more spontaneous magnetization. Therefore, when one 
crosses the critical temperature, ones observes a jump 
of the magnetization from </> Mln (r c ) to 0, which is the 
consequence of the competition between two minima of 
the potential, see Fig. 

For r > r c , the field configuration tfi Mln (r) which is no 
longer the ground state becomes a metastable state. One 
sees from Eq. (IE2|I that, for r > r* = 1, this metastable 
state disappears and we are left with <j) = as the only 
physically relevant state, see Fig. ^] Finally, it must be 
noted that, for r = r* , the curvature of the potential at 
the configuration Min (r*) vanishes: U"((j) Min (r*)) = 0. 
This means that the susceptibility in the metastable state 
diverges at r = r* . Similarly, one can show that the 
correlation length in the metastable state also diverges. 

Let us now come back to the NPRG method. In the 
truncation of the effective average action that we use — 
an expansion in powers of the fields of the form Eq. 1147|) 

- we retain only local informations on the potential 
around its nontrivial minimum — which is equivalent 
to the configuration </> Mln (r) discussed above. In par- 
ticular we do not accurately describe the physics around 
the zero-field configuration <fri — <f>2 — 0. We are thus 
unable to compare the energies of different local minima 
and to determine the temperature of transition at which 
the energies of the two minima are equal. Also, in a small 
domain of temperatures — equivalent here to r c < r < r* 

— the configuration that we probe corresponds actually 
to the metastable state and not to the true equilibrium 
state. However, these phenomena should not induce a 
large bias in our analysis as long as the transition is 
weakly of first order since, in this case, the temperature 
range where metastable states exist is very small. This 
means that the error induced on the determination of the 
critical temperature is very small too. 

Moreover in our study, as in the toy model above, we 
should observe when r reaches r* — the temperature at 
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which the metastable state must disappear — the as- 
sociated divergence of the correlation length discussed 
previously. This is precisely what we found in frustrated 
magnets for small reduced temperature — see Fig. 113b . 

Note that this increase of the correlation length as well 
as the error associated with our determination of the crit- 



ical temperature both rely on the truncation that we have 
considered. These problems can be cured by considering 
truncations of the form Eq. 1145|l which retains the full 
field-dependence of the potential (see for instance [45] for 
a treatment of a first order phase transition in a NPRG 
approach) . 
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